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FOREWORD 


The  rapid  advances  made  in  helicopter  and  V/STOL  aircraft 
development  in  the  past  few  years  have  spotlighted  areas  in  which 
significant  aerodynamic  problems  have  been  encountered,  and  in  some 
cases  the  problems  still  exist.  Fortunately,  a  technological  maturity 
has  now  been  attained  in  the  industry,  making  possible  investigations 
through  knowledgeable  scientific  approaches  that  have  produced  an 
enlightened  understanding  of  the  problems  and,  in  many  cases,  have 
led  to  practical  solutions.  The  next  generation  of  flight  vehicles, 
currently  in  the  design  and  development  stages,  are  offering  challenges 
to  the  aerodynamicist  and  engineer,  and  it  is  evident  that  future  vehicle 
developments  will  demand  an  ever-increasing  rate  of  technological 
advance  in  the  knowledge  and  understanding  of  aerodynamic  phenomena. 

Several  years  have  passed  since  a  technical  specialists*  meeting 
was  held  to  direct  attention  specifically  to  the  low-speed  aerodynamic 
problem  areas  of  helicopters  and  V/STOL  vehicles.  Therefore,  in  the 
interest  of  disseminating  up-to-date  information,  the  cosponsors  of  this 
symposium,  the  U.S.  Army  Aviation  Materiel  Laboratories  (USAAVLABS) 
and  Cornell  Aeronautical  Laboratory,  Inc.  (CAL),  believe  such  a  meeting 
among  technical  specialists  in  the  field  would  be  timely.  It  is  hoped  that 
this  symposium  will,  through  the  presentation  of  selected  technical 
papers,  establish  the  state  of  the  art  of  aerodynamic  analysis  in  the 
basic  problem  areas  and  spotlight  those  critical  areas  where  research 
is  urgently  needed.  The  ultimate  objective  is  to  identify  those  advances 
required  in  the  state  of  the  art  that  can  assure  the  availability  of  the 
analytical  tools  needed  to  develop  and  analyze  the  next  generation  of 
helicopters  and  V/STOL  aircraft. 

In  keeping  with  these  objectives,  five  technical  sessions,  each 
dealing  with  specific  basic  areas  of  aerodynamic  analysis  associated 
with  V/STOL  aircraft,  have  been  formed.  In  addition,  a  panel  session 
has  also  been  scheduled  in  which  outstanding  members  of  industry  and 
government  from  three  countries  will  present  their  recommendations 
for  areas  of  research  that  need  to  be  pursued  if  the  state  of  the  art  is 
to  advance  at  the  required  rate. 

It  is  believed  that  the  formal  presentation  of  th«  selected  papers 
and  the  panel  presentations  and  recommendations  will  stimulate  con¬ 
structive  discussions  among  the  specialists  in  attendance.  While  a 
significant  amount  of  time  has  been  allotted  during  the  sessions  for 
this  to  be  accommodated,  it  is  hoped  the  sessions  will,  in  addition, 
stimulate  discussions  and  serious  thought  between  the  attendees  and 
the  technical  members  of  the  various  organizations  who  were  unable 
to  attend.  In  order  to  foster  this  latter  objective,  the  informal  discussions 
(questions  and  answers,  etc.)  will  be  recorded,  printed,  and  mailed  to 
each  attendee  for  inclusion  in  his  copy  of  the  proceedings. 


The  proceedings  of  the  symposium  have  been  bound  in  four 
volumes  ~  a  separate  volume  for  the  technical  sessions  of  each  day,  and 
one  volume  for  the  panel  sessions.  We  are  indebted  to  the  authors  for 
preparing  their  manuscripts  in  a  form  that  could  be  directly  reproduced. 
This  material  was  published  as  provided  by  the  authors  and  was  neither 
checked  nor  edited  by  CAL  or  USAAVLABS. 

The  cosponsors  of  the  symposium  are  grateful  to  the  many  people 
who  contributed  to  its  success.  In  particular,  our  thanks  go  to  Colonel 
Harry  L.  Bush,  Commanding  Officer  of  the  U.S.  Army  Aviation  Materiel 
Laboratories,  and  Mr.  Ira  G.  Ross,  President  of  Cornell  Aeronautical 
Laboratory,  Inc.,  who  opened  the  sessions;  to  Mr.  Charles  W.  Harper, 
our  keynote  speaker;  to  Major  General  William  Bunker,  Deputy 
Commanding  General,  U.S.  Army  Materiel  Command,  for  his  address 
at  the  symposium  banquet;  to  the  five  session  chairmen  -- 
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INTRODUCTION 

A  theoretical  aerodynamic  analysis  of  a  ducted  propeller 
at  angle  of  attack  is  presented  in  References  1  and  2.  This  po¬ 
tential  flow  analysis  can  be  used  to  predict  the  force  and 
moment  on  the  duct  in  terms  of  the  propeller  thrust.  The  useful¬ 
ness  of  any  theory  is  dependent  on  evaluation  by  comparison  with 
valid  data.  Much  of  the  previous  data  (summarized  in  Ref.  1)  are 
restricted  to  hovering  or  axial  flight  and/or  low  Reynolds  number 
with  flow  separation  from  the  duct.  Full-scale  test  data  for  two 
ducted-propeller  configurations  at  angle  of  attack,  presented 
herein,  allow  a  correlation  study  and  a  good  evaluation  of  the 
potential  flow  theory.  The  duct  force  and  moment  are  predicted 
by  use  of  Reference  2  and  compared  with  experiment.  The  analysis 
is  extended  to  predic  the  duct  pressure  distribution  and  boundary- 
layer  separation  which  are  also  compared  with  experiment. 

The  analysis  is  based  on  a  combination  of  two  previous 
potential  flow  theories  for  a  ring-wing  at  a  small  angle  of 
attack,  Reference  3,  and  for  a  lightly  loaded  ducted  propeller  in 
axial  flow,  Reference  4.  In  Reference  2,  the  duct  thrust,  normal 
force,  and  pitching  moment  were  predicted  by  representing  the 
propeller  as  a  uniformly  loaded  actuator  disk.  Estimates  of  the 
normal  force  and  pitching  moment  on  propeller  blades  operating  in 
a  duct  were  small  compared  with  those  on  the  duct.  Small  amounts 


of  duct  taper,  camber,  and  thickness  were  estimated  to  have 
little  effect  on  the  duct  force  and  moment;  hence,  the  duct  is 
represented  herein  as  a  thin  cylinder  to  predict  its  force  and 
moment. 

For  comparison  with  test  data,  the  theory  is  usually 
used  to  predict  the  thrust,  normal  force,  and  pitching  moment  on 
the  cylinder  in  terms  of  the  measured  total  thrust  (on  both  the 
disk  and  cylinder)  .  The  predicted  duct  thrust  can  then  be  com¬ 
pared  with  the  value  deduced  from  the  measured  duct  pressure  dis¬ 
tribution.  The  predicted  duct  normal  force  and  pitching  moment 
are  compared  with  the  values  measured  for  the  ducted  propeller 
unit. 

The  work  described  herein  was  sponsored  by  the  Ames 
Research  Center  of  NASA  under  Contract  No.  NAS 2-2647. 


EXPERIMENTAL  PROGRAM 

Two  full-scale  ducted  propellers  have  been  tested  at 
angle  of  attack  in  the  NASA  Ames  Research  Center  40-  by  80-foot 
wind  tunnel.  The  first  of  these  was  the  "Doak"  ducted  propeller 
mounted  on  the  tip  of  a  semispan  wing  (Ref.  5) .  This  unit,  which 
is  referred  to  herein  as  the  4- foot  model,  consists  of  a  duct  with 
a  chord-to-diameter  ratio  of  0.608  and  a  profile  thickness-to- 
chord  ratio  of  0.158;  an  eight-bladed  propeller  with  fixed  blade 

Q 

pitch  of  15  at  the  tip;  a  set  of  seven  inlet  guide  vanes  of  vari¬ 
able  pitch;  and  a  set  of  nine  stators  aft  of  the  propeller.  A 
photograph  and  a  sketch  of  the  cross  section  of  this  ducted  pro¬ 
peller  are  shown  in  Figures  1(a)  and  2(a).  Complete  details  of 
the  4-foot  model  and  the  test  methods  are  described  in  Reference 
5.  The  second  unit  tested  was  the  "Bell  X22-A"  ducted  propeller 
mounted  on  the  tip  of  a  stub  wing.  This  unit,  referred  to  as  the 
7-foot  model,  consists  of  a  duct  with  a  chord-to-diameter  ratio 
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of  0.525  and  a  profile  thickness-to-chord  ratio  of  0.172;  a 
three-bladed  variable  pitch  propeller;  and  six  streamlined  support 
struts  aft  of  the  propeller.  For  some  of  the  tests,  ar  elevon 
extended  across  the  duct  exit  plane  with  a  thickness  of  0.44 
feet  and  a  chord  of  approximately  5  feet.  A  photograph  and  a 
sketch  of  the  cross  section  of  this  model  are  shown  in  Figures 
1(b)  and  2(b)  .  The  test  setup  and  techniques  were  similar  to 
those  used  with  the  4-foot  model. 

For  the  4-foot  model  the  data  included  direct  measure¬ 
ment  of  lift,  drag,  and  pitching  moment  reported  in  Reference  5 
and  duct  surface  pressure  at  27  orifice  stations  around  the  duct 
profile  and  at  10  azimuth  locations.  Total  pressure  was  measured 
along  the  duct  radii  directly  ahead  of  the  inlet  vanes,  directly 
behind  the  propeller,  and  across  the  duct  exit  plane  with  pres¬ 
sure  rakes  at  10  azimuth  locations. 

Data  for  the  7-foot  model  included  direct  measurements 
of  forces  and  moments  and  duct  surface-pressure  measurements  at 
ID  orifice  stations  around  the  duct  profile  for  3  azimuth  loca¬ 
tions.  Since  the  duct  pressure  data  were  taken  at  only  3  azi¬ 
muth  locations,  pressure  integrations  were  made  to  obtain  duct 
thrust  only  at  zero  angle  of  attack. 

In  addition  to  the  data  described  above,  stalling  of  the 
4-foot  model  duct  was  deduced  from  observation  of  tufts,  sound 
level,  and  pitching  moment  as  described  in  Reference  5. 

METHOD  OF  PREDICTION 
Potential  Flow  Model 

For  a  ducted  propeller  at  angle  of  attack,  the  analysis 
of  Reference  2  predicts  the  thrust,  normal  force,  and  pitching 
moment  on  the  duct  in  terms  of  the  propeller  thrust.  The  analysis 
represents  the  propeller  and  stator  blading  and  the  centerbody 
as  a  uniformly  loaded  actuator  disk.  There  is  no  swirl  in  the 
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slipstream,  and  the  circulation  about  each  propeller  and  stator 
blade  is  assumed  to  be  invariant  with  radius  and  time.  The 
boundary  layer  and  frictional  drag  of  the  duct  are  neglected. 

The  ducted  propeller  is  assumed  isolated  in  the  free  stream. 

Calculative  examples  have  been  run  to  investigate  the 
effects  of  duct  camber  and  radial  variation  of  disk  loading  in 
order  to  arrive  at  a  simple,  yet  realistic  analytical  flow  model. 
As  a  result  of  this  work,  the  following  assumptions  will  be  em¬ 
ployed: 

(a)  The  actuator  disk  loading  is  applied  only  to  the 
annular  area  swept  out  by  the  propeller  blades  and  not  to  the  aft 
part  of  the  centerbody  or  inner  duct  surface. 

(b)  The  duct  thrust  is  taken  as  that  due  to  leading- 
edge  suction  on  a  thin  cylinder. 

(c)  The  effects  of  duct  camber  and  thickness  are  neg¬ 
lected  except  in  computing  the  duct  pressure  distribution  and 
boundary  layer.  (It  is  shown  in  Ref.  2  that  this  is  a  good 
approximation  for  the  duct  for^^  and  moment.) 

(d)  The  propeller  ard  stators  are  represented  as  a 
uniformly  loaded  actuator  disk.  The  centerbody,  support  wing, 
and  elevon  are  neglected. 

Duct  Force  and  Moment 

To  predict  the  duct  force  and  moment,  the  duct  is  repre¬ 
sented  as  a  thin  cylinder  through  the  actual  duct  trailing  edge 


Sketch  A.-  Flow  model  for  predicting 
duct  force  and  moment. 
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All  the  trailing  vorticity  is  placed  on  the  extension 
of  the  cylinder.  This  approximation  requires  the  free-stream 
crossflow  to  be  small  compared  with  the  slipstream  velocity 
(V  sin  a  «  Vj)  .  The  pressure  jump  across  the  actuator  disk 
Ap  is  equal  to  the  increased  total  pressure  in  the  slipstream 
as  given  by 

Ad  V1  .  f V  cos  a  +  v"\2 

f’~ r  ' 1  ■  L — v — )  - 1 

=  2  ?r  cos  a  +  -  sin2  a  (1) 

v  V2 

The  propeller  thrust  is  taken  as 

TP(D)  =  Ap  Ap 


where  Ap  is  the  annular  area  swept  out  by  the  propeller  blades. 
Hence,  the  propeller  thrust  coefficient  is 


TP(P)  _ 
Aq  A 


The  coefficients  for  the  duct  force  components  and  pitching 
moment  are  given  by  Equations  (20)  through  (23)  of  Reference  2  as 

^DIP)  '  fl  Bi"  °  C“S  “  +  0  (4) 

SD(p)  *  £3  Bi"2  °  +  f4  CO"  (5) 

=  4fg  sin  a  cos  a  t  (fgf6  +  f?)  ^  sin  a  (6) 


where  the  f  coefficients  are  given  versus  the  duct  chord-to 
n 

diameter  ratio  c/D  in  Table  I,  Reference  2. 
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Use  of  blade  element  theory  in  Reference  2  indicated 
that  the  normal  force  and  pitching  moment  on  the  propeller  are 
generally  small  compared  with  those  on  the  duct.  Therefore,  we 
shall  take  the  coefficients  for  the  ducted  propeller  unit  as 


(7) 


D(P) 


(8) 


c  =  c 

"■dp  ®D(P) 


(9) 


The  strength  of  the  vortex  cylinder  surrounding  the 
slipstream  is  found  by  substituting  Equations  (1) ,  (3) ,  and  (5) 
into  (7)  and  solving  for 


T 

V 


The  general  procedure  for  prediction  will  be  to  deter¬ 
mine  y/V  from  Equation  (10)  using  a  measured  value  of 
and  then  to  evaluate  the  preceding  expressions.  This  value  of 
/V  is  also  used  to  predict  the  duct  surface-pressure  distribu- 
:  Lon  and  boundary  layer. 


Duct  Pressure  Distributions 

The  duct  pressure  distribution  is  obtained  by  first  pre¬ 
dicting  the  bound  vorticity  and  velocity  distribution  over  the 
thin  duct  camberline,  then  including  the  effect  of  duct  thickness 
to  obtain  the  surface  velocity  distribution,  and  finally  using 
Bernoulli's  law  to  get  the  duct  surface  pressure. 
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The  flow  model  used  to  predict  the  duct  surface  pres¬ 
sure  is  shown  in  Sketch  B. 

Numbers  1  through  4  designate 
the  camberline  or  duct  surfaces 
indicated. 


o 


Sketch  B.-  Flow  model  for  predicting  duct  pressure  distribution. 


Duct -bound  vortici^ 


All  the  bound  and  trailing  vorticity  is  placed  on  a  ref¬ 
erence  cylinder  through  the  duct  trailing  edge.  The  distribution 
of  duct-bound  vorticity  is  composed  of  elementary  yd  vortex 
rings  which  are  axially  symmetric  and  elementary  y  rings  with 
strength  proportional  to  cos  $.  The  axially  symmetric  part  of 
the  flow  field,  composed  of  V  cos  a  plus  that  induced  by  the 
7d  rings  and  the  y  cylinder,  is  set  tangent  to  the  duct 
camberline.  This  boundary  condition  is  actually  imposed  at  the 
duct  reference  cylinder  and  it  determines  the  strength  of  the 
Yd  rings. 

The  ya  rings,  together  with  the  axial  vortex  filaments 
which  trail  from  them,  are  required  to  cancel  the  velocity  across 
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the  duct  reference  cylinder  due  to  the  free-stream  crossflow 
V  sin  a. 

The  slope  of  the  duct  camber line  is  expressed  as  a  four- 
term  cosine  series 


or  _ 


cos  nd 


where 


cos  6  =  - 


Integration  gives  the  duct  camberline  shape  as 


rs  -R 


*o  (^f^)  ♦  R,  C1  ~  4°' 2  0 

+  Rg  cos  9  -  -j  cos3  9^ 

+  R3  Q-  %  +  |  cos2  6  -  ■ |  cos4  9^) 


The  Rn  coefficients  are  determined  such  that  Equation  (12)  fits 
the  geometric  camber line  at  x/c  =  0,  1/4,  1/2,  3/4,  and  1.  The 
actual  duct  profiles  and  the  approximate  analytical  profiles  are 
compared  in  Figure  3.  An  effective  camber  due  to  thickness  €t 
as  subtracted  from  Equation  (11)  is  described  in  Reference  6. 

This  effect  is  due  to  the  radial  velocity  induced  by  the  source 
rings  used  to  generate  duct  thickness.1  We  shall  take  €t  for  a 

■’’Evaluation  of  the  axial  velocity  induced  by  the  source  rings 
from  Reference  6  indicates  that  it  is  nearly  the  same  as  for  a 
straight,  2-dimensional  wing  with  the  same  thickness  as  the 
duct.  This  approximation  will  be  employed  in  the  present 
analysis. 


NACA  0018  thickness  distribution  using  Equation  (44) ,  Reference 
6. 2  Thus,  the  slope  of  the  effective  cantoerline  is  specified  as 

3 

ee  =  e  -  et  =  ^  Rn*  008  n6  (13) 

o 


To  set  the  flow  tangent  to  the  camberline,  we  specify 


that 


VT  +  VYd  €e 


(v  COS  a  +  y 


The  u  and  v  terms  are  the  axial  and  radial  velocity  compo¬ 
nents  induced  along  the  duct  reference  cylinder  by  the  y  cyl¬ 
inder  extending  from  the  duct  trailing  edge  and  the  rings 

bound  to  the  duct  reference  cylinder.  These  terms  are  expressed 
as  a  six- term  cosine  series  as  follows.  The  components  induced 
by  the  vortex  cylinder  can  be  found  from  Reference  7  in  terms  of 
elliptic  integrals.  Fourier  analysis  by  machine  then  give 


B  *  cos  n0 
n 


(15) 


cos  n0 


(16) 


2The  duct  profiles  for  the  two  ducts  do  not  correspond  to  any 
standard  airfoil  section.  However,  it  was  found  that  the  NACA 
0018  is  a  good  approximation  to  the  actual  thickness  distribu¬ 
tion  of  both  configurations,  particularly  over  the  forward 
portion  of  the  airfoil  where  the  pressure  distributions  are  of 
greatest  interest. 
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If  the  duct -bound  vorticity  is  put  in  terms  of  a  Glauert  series 

as3 


sin  nS 


(17) 


then  the  corresponding  induced  axial  velocity  component  is  from 
Equation  (18) ,  Reference  2 

C3  ‘  C1  C4  *  C2 

+  - ^ -  CC3  20  +  - 2 -  cos  3® 

4  6 

C  -  C 

+  -2—q — -  cos  40  (18) 


The  radial  velocity  component  can  be  found  from  Equation  (14)  of 
Reference  4  as 


Tables  2.1  through  2.4  of  Reference  4. 

The  six  unknown  CR  coefficients  appearing  for  yd  are 


determined  by:  (1)  substituting  the  foregoing  cosine  series  for 
the  terms  in  Equation  (14),  (2)  expanding  the  right-hand  side  as 


Between  the  actuator  disk  and  duct  trailing  edge  Yd  includes 
the  vorticity  shed  from  the  disk  onto  the  inner  duct  surface. 
The  Kutta  condition  of  no  flow  across  the  duct  trailing  edge  is 
satisfied  because  Yn  Is  nonsingular  at  6  -  -w  in  Equation 
(17).  D 
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a  cosine  series,  and  (3)  equating  each  of  the  six  harmonics  of 
Equation  (14) .  The  resulting  six  linear  algebraic  equations  are 
then  solved  by  machine  for  the  six  CR.  The  remaining  part  of 
the  duct  vorticity  y  is  the  same  as  determined  by  Weissinger 
for  a  ring-wing  at  angle  of  attack,  and  is  given  by  Equation  (8) 
of  Reference  2  as 


V  sin  a  cos  <t> 


c  cot 
o 


!  +  Is 


sin  n© 


Weissinger1 s  values  for  the  cr  coefficients  are  given  in 
Table  1,  Reference  2.  The  axial  velocity  induced  by  ya  along 
the  duct  reference  cylinder  is  given  by  Equation  (19)  of  Refer¬ 
ence  2  as 


cV  sin  a  cos  <p 


^¥-Hc°+®+^Cc°+%) 

+  [co  •  i4D6£  (s  +t)  +  t]  cos  6 


C3  '  C1  C4  "  C2 

+  - - cos  26  +  — -  cos  3 6 

4  6 


CS  '  C3 

+  - g -  cos  40 


Finally,  the  velocity  distribution  over  the  four  duct  camberline 
surfaces  shown  in  Sketch  B  can  be  evaluated  by  substitution  of 
preceding  expressions  into 


s +  s 

V  *  - - 2 - 


+  V  cos  a  +  u  +  u  ,  +  u 

a/  f\f  ft 
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where  the  plus  sign  refers  to  the  inner  surfaces  2  and  3  and  the 
minus  sign  the  outer  surfaces  1  and  4.  At  the  sides  of  the  duct, 
7  induces  a  tangential  velocity  component.  We  shall  not  con¬ 
sider  the  component,  since  it  varnishes  at  the  vertical  plane  of 
symmetry  to  which  we  shall  restrict  the  pressure  distribution 
predictions. 


Duct  surface  velocity 

To  obtain  the  duct  surface  velocity  from  the  velocity 
distribution  over  the  camberline  (Eq.  (22)),  we  apply  corrections 
to  the  discontinuous  (vorticity)  and  continuous  terms  in  Equa¬ 
tion  (22) .  These  correction  factors  are  deduced  as  follows  from 
theoretical  results  for  two-dimensional  (c/D  -  0)  thick  airfoils 
in  Reference  8.  For  a  thin  flat  plate  at  angle  of  attack,  the 
bound  vorticity  is 

^  =  2a  cot  |  (23)^ 

At  a  =  1/(2tt)  the  lift  coefficient  is  unity  and 


(24) 


This  chordwise  distribution  is  listed  in  Table  I  for  comparison 
with  values  for  a  NACA  0018  airfoil  from  Reference  8.  It  is 
evident  that  thickness  has  little  effect  on  the  bound  vorticity 
except  near  the  airfoil  nose  where  x  <  0.1c.  Hence,  we  shall 
correct  the  singular  (cot(6/2))  terms  in  y D  and  y by  re¬ 
placing  cot (6/2)  by  the  corresponding  value  for  a  NACA  0018  air¬ 
foil  when  x  <  0.1c.  For  example,  at  the  leading  edge 
(1/2t r)  cot(6/2)  is  replaced  by  1.342.  The  surface  velocity  dis¬ 
tribution  for  a  two-dimensional  NACA  0018  airfoil  at  a  =  0  is 
shown  in  the  last  column  of  Table  I.  In  this  case,  the  continu¬ 
ous  part  of  the  camberline  velocity  distribution  is  V  and  the 
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continuous  part  of  the  surface  velocity  distribution  is  VF(x) , 
where  F(x)  is  given  by  the  last  column  in  Table  I.  By  analogy, 
we  shall  obtain  the  continuous  part  of  the  duct  surface  velocity 
by  multiplying  the  continuous  part  of  camber line  velocity  by 
F(x) .  Thus,  we  obtain  the  duct  surface-velocity  distribution 
from  the  duct  camber line  distribution  (Eq.  (22))  as 


u 


=  + 


F(X) 


(25) 


The  first  parenthesized  term  is  corrected  only  when  x  <  0.1c  as 
described  previously. 


Duct  pressure  distribution 

The  duct  surface-pressure  coefficient  is  obtained  from 
Equation  (25)  by  Bernoulli's  law  as 


(26) 


However,  on  the  inner  duct  surface  downstream  of  the  actuator 
disk,  the  total  pressure  is  increased  by  Ap  and  Cp  is  in¬ 
creased  by  Ap/q  due  to  the  pressure  jump  across  the  disk. 


Duct  Boundary-Layer  Model  and  Assumptions 

Our  purpose  in  computing  the  duct  boundary  layer  is  to 
predict  those  conditions  under  which  separation  will  occur  on  the 
windward  inside  surface  before  the  propeller.  We  shall  use  the 
axisymmetric  Truckenbrodt  method  to  obtain  the  momentum  thickness 
(6)  and  shape  factor  (H)  of  the  boundary  layer  and  then  we 
shall  use  an  empirical  rule  (H  =  1.8)  to  predict  turbulent 
separation. 
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The  Truckenbrodt  method  (Refs.  9  or  10)  has  been  com¬ 
pared  to  known  results  in  both  laminar  and  turbulent  flow  and 
found  to  give  excellent  agreement.  We  have  computed  the  laminar 
boundary- layer  parameters  for  Howarth  flow  on  a  flat  plate  and 
they  agree  within  5  percent  of  the  exact  Howarth  solution.  We 
have  also  computed  turbulent  boundary- layer  parameters  for  a 
NACA  65(216) -222  profile  at  10.1°  angle  of  attack  using  the  mea¬ 
sured  velocity  distribution.  The  computed  momentum  thickness, 
shape  factor,  and  separation  point  all  agree  well  with  measured 
results. 

We  shall  make  use  of  the  following  assumptions  in  calcu¬ 
lating  the  duct  boundary  layers: 

1 3.)  The  boundary- layer  momentum  thickness  is  small  in 
comparison  with  the  duct  radius  so  that  the  presence  of  the 
boundary  layer  does  not  appreciably  alter  the  potential  flow 
around  the  duct. 

(2)  The  predicted  surface  velocity  (Eq.  (25))  is  taken 
as  the  boundary- layer-edge  velocity  distribution. 

(3)  The  velocity  profiles  in  the  boundary  layer  are 
approximated  by  a  one-parameter  family  of  curves. 

As  a  transition  criterion,  it  is  assumed  that  the  bound¬ 
ary  layer  changes  from  laminar  to  turbulent  at  a  Reynolds  number 
of  about  400  based  on  momentum  thickness  and  local  surface 
velocity.4  In  cases  with  high  adverse  pressure  gradients  near 
the  leading  edge,  the  laminar  boundary  layer  separated  (H  =  4.04) 
before  Re^  =  400.  Assuming  that  the  separated  region  acts  as  a 
turbulence  trip,  the  point  of  transition  is  located  at  the  first 
indication  of  laminar  separation.  In  Reference  11,  it  is  reported 

4 

The  transition  criterion  is  based  on  a  measured  critical  Reynolds 
number  on  a  flat  plate,  Vx/v  =  3xl05  (Ref.  10) .  This  corresponds 
to  V6/V  »  365  in  a  laminar  boundary  layer  on  a  flat  plate. 
Since  the  location  of  turbulent  separation  is  insensitive  to 
small  variations  in  transition  location,  transition  was  assumed 
to  occur  when  350  <  us6/v  <  425. 


that  small  laminar  separation  bubbles,  which  did  not  show  up  in 
the  pressure  distribution,  did  serve  as  turbulence  trips. 

The  shape  factor  (H)  is  used  to  estimate  the  location 
of  boundary- layer  separation  from  the  duct.  The  usual  criteria 
for  boundary- layer  separation  is  that  turbulent  separation  occurs 
when  H  has  attained  a  value  between  1.8  and  2.4,  and  laminar 
separation  occurs  when  H  =  4.04. 

The  Truckenbrodt  solution  to  the  boundary  layer  is  based 
on  integration  of  the  energy  equation.  The  energy  equation  is 
written  in  the  form 


(27) 


where  6**  is  the  energy  thickness  of  the  boundary  layer,  d  is 
the  dissipation  term,  and  t  is  the  energy  of  turbulent  motion 
term.  The  right  side  of  Equation  (27)  is  the  shear  stress  work 
in  the  boundary  layer.  Truckenbrodt ' s  final  expression  for  mo¬ 
mentum  thickness  as  a  function  of  boundary- layer-edge  conditions 
is 


(28, 

where 
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laminar  flat-plate  drag  coefficient  from  the 
Blasius  relation  based  on  q  and  s. 
turbulent  flat-plate  drag  coefficient  from  the 
Schultz-Grunow  relation  based  on  q  and  s. 

laminar  flow:  n  =  1 
turbulent  flow:  n  =  1/6 

The  actual  calculation  of  momentum  thickness  is  carried 
out  in  a  surface  coordinate  system  where  z  is  measured  from  the 
leading  edge  along  the  surface  of  the  duct  and  s  is  the  surface 
length  of  the  duct.  The  point  of  transition  of  the  boundary 
layer  from  laminar  to  turbulent  is  denoted  by  z .  After  the 
boundary  layer  is  computed  in  the  surface  coordinate  system,  the 
results  are  transferred  into  a  coordinate  system  based  on  axial 
distance  x  and  chord  length  c. 

Boundary-layer  parameters  were  computed  on  the  inside 
surface  of  both  ducts  using  measured  and  predicted  velocity  dis¬ 
tributions.  Computation  of  these  parameters  using  the  measured 
velocity  distribution  was  started  at  the  leading  edge.  When  the 
predicted  velocity  distribution  was  used  to  compute  the  boundary- 
layer  parameters,  the  computation  was  started  either  at  the  lead¬ 
ing  edge  or  at  the  stagnation  point  on  the  outside  surface. 

RESULTS 

The  predicted  force  and  moment  coefficients  are  given 
by  Equations  (3)  through  (10)  in  terms  of  the  parameters  Ap/A, 
c/D,  and  fn>  Values  for  these  parameters  are  listed  in  Table  II 
for  both  ducts.  The  f  are  found  by  interpolation  from  Table  1, 
Reference  2.  The  total  measured  force  and  moment  coefficients 
are  deduced  from  force  balance  data.  Individual  duct  and  pro¬ 
peller  thrust  coefficients  are  determined  from  integrations  of 
measured  pressure  distributions  over  the  duct  surface  ar.d  total 
pressure  rake  data  just  aft  of  the  propeller. 
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Division  of  Thrust  in  Axial  Flow 


Four- foot  model 

For  a  series  of  axial  flow  runs  with  the  4-foot  model 

duct,  the  measured  thrust  coefficients  are  listed  in  Table  III, 

The  sum  of  the  individual  duct  and  propeller  measured  thrusts 

C_  and  C_  is  equal  to  the  total  measured  thrust  C_ 

rD(P)  tP(D)  tDP 

within  7  percent.  The  values  of  y/V  and  GpD^pj  predicted 

from  the  measured  C_  ,  are  listed  next.  Finally,  the  values 

AP(D) 

of  y/v  and  C_  predicted  from  the  measured  C_  ,  .  are 

tD(P)  ^PjD) 

shown.  The  latter  values  are  obtained  by  use  of  Equations  (3) , 

(1)  ,  and  (5)  in  succession  rather  than  by  the  usual  method  (Eqs. 

(10)  and  (5))  . 

Seven-foot  model 


Measured  values  of  C,_  and  C_  are  shown  in 

TDP  ^D(P) 

Table  III  together  with  the  predicted  value  of  ctD(P) *  The 
effect  of  increased  propeller  blade  pitch  is  indicated  by  the 
last  three  runs. 


Duct  Force  and  Moment  at  Angle  of  Attack 

Four-fcot  model 

Normal  force  and  pitching  moment  coefficients  predicted 
versus  the  total  thrust  coefficient  are  shown  in  Figures  4  and 
5,  together  with  measured  values.  The  data  generally  lie  well 
below  the  predicted  curves,  particularly  the  normal-force  coef¬ 
ficient.  The  measured  duct  normal  force  (deduced  from  the  pres¬ 
sure  distribution)  is  considerably  smaller  than  the  measured 
total  normal  force  as  shown  in  Table  IV. 

Seven-foot  model 

Similar  predicted  and  measured  coefficients  for  the 
larger  duct  are  shown  in  Figures  6  and  7.  The  data  are  all  for 
the  same  propeller  pitch  and  with  the  elevon  either  off  or 
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aligned  with  the  duct  axis.  The  elevon  appears  to  have  most 
effect  on  the  measured  normal  force  at  the  highest  angles  of 
attack. 

Duct  Pressure  Distributions  in  Axial  Flow 

Following  the  procedure  described  previously,  the  first 
step  in  predicting  the  duct  pressure  distribution  is  to  fit  Equa¬ 
tion  (12)  to  the  duct  camber line  and  solve  for  the  four  camber- 
line  coefficients  Rn<  The  computed  values  for  Rn  are  given 
in  Table  II  for  both  ducts.  The  thickness  distributions  for  both 
ducts  were  approximated  by  a  NACA  0018  profile.  The  resulting 
duct  profiles  (Fig.  3)  fit  the  nose  radius  of  curvature  for  both 
ducts.  However,  the  hook  in  the  camber line  near  the  leading 
edge  of  the  larger  duct  is  not  well  represented.  Attempts  to  fit 
the  hook  with  the  four-term  series  caused  a  very  poor  fit  over 
most  of  the  camber line. 

The  effective  camber  due  to  thickness  (et  in  Eq.  (13)) 
is  given  by  Reference  6  for  a  NACA  0008  profile  and  c/D  =  1/2. 
Multiplication  by  18/8  gives,  for  an  0018  profile  and  c/D  =  1/2, 


=  0.001  +  0.040  cos  6  +  0.013  cos  2 G  -  0.001  cos  36 

This  distribution  was  used  for  both  ducts.  The  resulting  Rn* 
coefficients  for  the  effective  camber line  (Eq.  (13))  are  given 
in  Table  II. 

Fourier  analysis  by  machine  gives  the  series  coeffi¬ 
cients  Bn  and  Bn*  which  appear  in  Equations  (15)  and  (16)  for 
the  velocity  components  induced  by  the  vortex  cylinder  trailing 
from  the  duct.  Machine  computed  values  are  given  in  Table  V  for 
the  two  ducts. 

The  coefficients  which  appear  in  Equation  (19)  are 

found  by  interpolation  from  Tables  2.3  and  2.4  of  Reference  4. 
Values  for  the  two  ducts  are  given  in  Tables  VI  and  VII. 
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The  Cn  coefficients  in  Equation  (17)  for  the  duct- 
bound  vorticity  yd  must  be  computed  for  each  case  since  they 
depend  upon  y/V.  However,  for  no  duct  camber  (R  *  =  0)  the  C 
depend  only  upon  c/D,  and  these  values  are  given  in  Table  V. 

The  cR  coefficients  in  Equation  (20)  for  the  duct- 
bound  vorticity  ya  can  be  found  by  interpolation  from  Table  1 
of  Reference  2.  Values  for  the  two  ducts  are  given  in  Table  V. 

Four- foot  model 

Duct  pressure  distributions  computed  for  the  smaller 
duct  in  axial  flow  are  shown  in  Figure  8  at  three  thrust  levels. 
Measured  data  are  shown  by  the  vertical  lines  to  indicate  the 
variation  with  azimuth.  The  discontinuous  change  in  predicted 
pressure  on  the  inner  duct  surface  is  caused  by  the  pressure 
jump  across  the  assumed  uniformly  loaded  actuator  disk. 

Seven-foot  model 

For  a  series  of  axial  flow  runs  with  the  larger  model, 
the  measured  and  predicted  duct  pressure  distributions  are  shown 
in  Figure  9.  The  first  three  runs  are  for  propeller  blade  pitch 
8  =  19°,  and  the  last  is  for  3  =  29°. 

Duct  Pressure  Distributions  at  Angle  of  Attack 

Four -foot  model 

Pressure  distributions  measured  and  predicted  for  the 
smaller  duct  at  angle  of  attack  are  shown  in  Figure  10.  The 
corresponding  pressure  distribution  for  a  =  0  is  shown  in  Fig¬ 
ure  8(b)  . 

Total  pressure  distributions  measured  vertically  across 

the  duct  exit  plane  are  shown  in  Figure  11.  These  data  show 

that  the  propeller  loading  was  not  uniform  and  that  the  flow  was 

separated  from  the  lower  inner  duct  surface  2  even  at  low  angles 
of  attack. 
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Seven-foot  model 


Pressure  distributions  measured  and  predicted  for  the 
7-foot  model  duct  at  angle  of  attack  are  shown  in  Figures  12  and 
13.  The  effect  of  angle  of  attack  at  a  high  advance  ratio  is 
shown  in  Figures  12  and  9(b) ;  and  at  a  low  advance  ratio  in  Fig¬ 
ure  13  and  9(c) .  All  the  data  at  angle  of  attack  were  taken  13° 
off  the  vertical  plane  of  symmetry  (0  =  0  in  Sketch  B) . 

Duct  Boundary  Layers 

With  the  foregoing  duct  surface  geometry  and  both  mea¬ 
sured  and  predicted  surface-pressure  distributions,  duct  boundary 
layers  were  computed  using  the  method  previously  described. 
Bernoulli's  equation  was  used  to  relate  surface  velocity  and 
pressure  distributions.  The  boundary  layers  were  computed  only 
on  the  lower  inner  surface3  of  the  computed  duct  profile. 

Four-foot  model 

For  the  smaller  duct  in  axial  flow  at  various  advance 
ratios,  the  bound ary- layer  momentum  thickness,  shape  factor,  and 
surface-velocity  distributions  are  shown  in  Figure  14.  The  mo¬ 
mentum  thickness  was  assumed  zero  at  the  leading  edge  in  each 
case.  The  turbulent  boundary  layers  computed  from  the  predicted 
velocity  distributions  did  not  separate  from  the  duct  ahead  of 
the  propeller.  However,  when  the  measured  velocity  distribution 
was  used  for  the  hovering  flight  condition  (J  =  0)  and  a  second 
high  thrust  condition  (J  =  0.178),  the  boundary  layer  separated 
slightly  ahead  of  the  propeller. 

The  boundary  layer  results  for  the  duct  at  a  =  80°  and 
J  =  0.166  are  shown  in  Figure  15,  where  the  measured  and  predic¬ 
ted  surface-velocity  distributions  (from  Fig.  10(d))  are  also 

5The  computed  duct  shape  is  the  superposition  of  a  NACA  0018 
thickness  distribution  on  the  duct  camberline  shape  from  Equa¬ 
tion  ( 12)  . 


20 


compared.  The  momentum  thickness  and  shape  factor  were  computed 
from  'he  leading  edge  using  both  predicted  and  measured  velocity 
distributions.  The  boundary  layer  computed  from  the  predicted 
velocity  distribution  separates  before  that  from  the  measured 
velocity  distributions.  For  the  predicted  velocity,  the  boundary 
layer  was  also  computed  from  the  stagnation  point  at  x/c  =  0.60 
with  almost  no  change  in  the  results. 

For  the  larger  duct  at  high  advance  ratio  (J  “  0.52)  and 
various  angles  of  attack,  the  bound ary- layer  momentum  thickness, 
shape  factor,  and  surface  velocity  distribution  are  shown  in  Fig- 
tire  16.  Using  predicted  velocity  distributions,  the  boundary 
layers  were  computed  from  the  leading  edge  at  a  =  0°  and  20°, 
and  from  x/c  =  0.05  near  the  outside  stagnation  point  at 
a  =  40°.  The  measured  velocity  distribution  was  also  used  for 
a  =  0°  and  40°.  Only  the  last  of  these  computed  boundary  layers 
separated.  This  result  agrees  with  the  measured  pressure  dis¬ 
tributions  where  separation  is  evident  at  a  =  40°  (Fig.  12(b)). 

For  the  measured  pressure  distributions  at  lower  advance 
ratio  (Fig.  13),  separation  was  computed  at  x/c  =  0.21  for 
u  ==  20°  and  x/c  =  0.18  for  a  =  40°.  The  effect  of  separation 
was  not  evident  in  the  measured  pressure  distribution  except  at 
higher  angles  of  attack. 

Duct  Stall  Boundary 

The  location  of  boundary- layer  separation  on  the  smaller 
duct,  as  computed  from  measured  and  predicted  velocity  distribu¬ 
tions,  is  shown  by  the  points  in  Figure  17(a) .  The  upper-dashed 
curve  is  the  boundary  corresponding  to  separation  midway  between 
the  leading  edge  and  the  propeller  (at  x/c  =  0.15).  This 

boundary  is  estimated  from  the  points  computed  with  measured 
velocity  distributions.  The  other  dashed  curves  correspond  to 
separation  at  other  values  of  x/c  as  estimated  from  the  com¬ 
puted  points.  Also  shown  in  this  figure  is  the  stall  boundary 
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and 


deduced  frow  the  measured  pitching  moment,  sound  level, 
visualization  of  tufts  on  the  inner  duct  lip  (Ref.  5)  . 

For  both  ducts  in  axial  flow,  the  location  of  boundary- 
layer  separation  computed  from  measured  velocity  distributions 
is  shown  in  Figure  17(b)  . 

All  of  the  computed  separation  points  on  both  ducts  are 
shown  in  Figure  18.  The  points  are  plotted  versus  a  Reynolds 
number  based  on  a  length  £  obtained  from  a  straigh  line  through 
the  separation  point  on  a  u(x)  plot  as  indicated  in  Sketch  C 
and  Figure  18. 


Sketch  C.-  Linear  approximation  to  duct  velocity  distributions. 

The  curved  line  in  Figure  18  shows  that  the  separation  points 

for  a  series  of  straight  line  velocity  profiles,  computed  by  the 

Truckenbrodt  method  as  before,  are  all  at  nearly  a  constant 

velocity  ratio,  u_„/u  *  0.65.  The  points  which  lie  below 

scp  max 

the  curve  are  for  velocity  profiles  which  are  concave  downward, 
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for  example,  curve  (a)  in  Sketch  C,  or  the  dashed  curve  in  Fig¬ 
ure  15.  On  the  other  hand,  the  points  which  lie  above  the  curve 
sure  for  convex  upward  profiles  such  as  curve  (b)  in  Sketch  C,  or 
the  lower  dashed  cruve  in  Figure  14. 

DISCUSSION 

Division  of  Thrust  in  Axial  Flow 

In  social  flow,  the  4-foot  model  duct  thrust  predicted 
both  from  the  measured  total  thrust  and  the  measured  propeller 
thrust  (Table  III)  agrees  with  the  measured  value  within  +10/-12 
percent.  This  agreement  seems  satisfactory  in  view  of  the  fact 
that  the  theory  neglected  the  observed  non-uniformity  of  propel¬ 
ler  blade  loading  and  the  prevalence  of  flow  separation  from  the 
duct  diffuser  in  addition  to  many  other  real  effects. 

For  the  7-foot  model  duct  in  axial  flow,  the  measured 
duct  thrust  generally  exceeds  the  predicted  value  (Table  III) . 

The  flow  separated  from  the  inner  duct  surface  upstream  of  the 

propeller  when  C_  >  2  according  to  the  bound ary- layer  pre- 

ADP 

dictions  using  the  measured  pressure  distributions  (Fig.  17) . 

The  predictions  did  not  account  for  the  hook  in  the  duct  camber¬ 
line  near  the  leading  edge  or  the  duct  frictional  drag.  However, 
calculations  indicate  that  both  of  these  effects  are  appreciable 
only  at  low  t'nrust  levels.  Calculations  also  indicate  that  the 
thrust  on  the  7-foot  model  (but  not  the  4-foot  model)  duct  is 
predicted  more  accurately  when  the  pressure  rise  across  the 
actuator  disk  is  assumed  to  act  on  the  aft  part  of  the  inner  duct 
surface. 

Duct  Force  and  Moment  at  Angle  of  Attack 

The  main  difference  between  the  measured  and  predicted 
performance  of  the  two  ducted  propellers  is  the  normal  force  at 
angle  of  attack.  Thu  measured  values  for  the  4-foot  model  duct 
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(Fig.  4)  are  generally  much  smaller  than  for  the  7-foot  model 
duct  (Fig.  6) ,  particularly  at  high  thrust  levels.  The  predicted 
normal  force  is  nearly  the  same  for  the  two  ducts  and  generally 
in  good  agreement  with  the  data  for  the  larger  duct  without  the 
elevon  (Fig.  6(a)).  The  data  and  predictions  both  indicate  that 
there  was  considerable  flow  separation  from  the  4-foot  model  duct 
at  angle  of  attack.  However,  it  appears  that  the  very  low  normal 
force  measured  for  the  smaller  duct  at  high  thrust  level  was 
caused  largely  by  the  concentration  of  propeller  thrust  loading 
toward  the  hub.  Except  for  this  effect,  one  might  expect  the  nor¬ 
mal  force  on  the  smaller  unit  to  be  greater  than  predicted  for  the 
duct  alone  because  of  the  large  centerbody  and  large  number  of 
stator  and  propeller  blades.  Theoretical  estimates  indicate 
that  concentration  of  the  propeller  loading  toward  the  hub  does 
cause  the  duct  normal  force  to  decrease  significantly  at  high 
thrust  level  but  has  little  effect  or.  the  duct  thrust.  The  pitch¬ 
ing  moment  appears  to  be  well  predicted  for  both  ducts  (Figs.  5 
and  7)  . 

Duct  Pressure  Distributions  in  Axial  Flow 

The  pressure  distributions  for  the  4-foot  model  duct 
(Fig.  8)  appear  to  be  reasonably  well  predicted.  At  the  higher 
thrust  levels,  and  particularly  for  the  hovering  flight  condition, 
the  pressure  on  the  inner  duct  surface  aft  of  the  propeller  is 
lower  than  predicted.  This  is  probably  due  to  the  fact  that  the 
propeller  loading  is  concentrated  near  the  hub  so  that  the  full 
disk  pressure  jump  is  not  applied  suddenly  to  the  duct  surface. 

The  pressure  distributions  for  the  larger  duct  (Fig.  9) 
are  accurately  predicted  on  the  outer  surface.  However,  on  the 
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more  critical  inner  surface,  the  pressure  is  predicted  less  accu¬ 
rately  than  for  the  smaller  duct.  This  is  particularly  true  when 
the  blade  pitch  and  thrust  coefficient  are  high  (Fig.  9(c)). 
Although  the  presence  of  the  elevon  appears  to  have  little  effect 
on  the  duct  pressure  distribution,  increased  blade  pitch  at  a 
fixed  thrust  lowers  the  measured  pressure  over  the  inner  duct 
surface.  The  effect  is  complicated  by  the  fact  that  separation 
of  the  boundary  layer  slightly  ahead  of  the  propeller  was  com¬ 
puted  for  all  of  the  pressure  distributions  shown  at  the  higher 
blade  pitch.  Such  separation  and  possible  reattachment  after  the 
propeller  would  probably  be  affected  by  the  loading  of  the  pro¬ 
peller  blade  tips.7 

The  discontinuity  in  duct  pressure  predicted  for  a  uni¬ 
form  actuator  disk  loading  is  not  evident  for  either  ducted  pro¬ 
peller.  This  can  be  attributed  to  the  low  loading  of  the  propel¬ 
ler  blade  tips  for  the  smaller  unit  as  mentioned  earlier.  For 
the  larger  unit  the  blade  tips  were  apparently  more  highly  loaded, 
but  the  computed  duct  boundary  layer  separated  near  the  propeller 
except  when  the  blade  loading  was  low.  The  separation  of  the 
duct  boundary  layer  may  have  caused  the  pressure  discontinuity  on 
the  inner  duct  wall  to  be  smoothed  out  similar  to  a  shock  wave¬ 
boundary  layer  interaction.  The  predicted  pressure  gradient  on 
the  inner  duct  surface  can  be  improved  by  smoothing  out  the  dis¬ 
continuity  across  the  actuator  disk. 

7 Recently,  rake  data  were  obtained  approximately  10  inches  up¬ 
stream  of  the  duct  exit  for  a  =  0,  3  =  29°,  J  =  0,  0.22,  0.45, 
0.62.  The  measured  total  pressure  was  very  uniform  for  each 
one  of  these  runs  except  within  the  duct  boundary  layer  (about 
2  in.  thick)  and  toward  the  centerline  (r  18  in.  where 
Ap  *  0) .  The  low  dynamic  head  in  the  central  part  of  the  slip¬ 
stream  is  believed  to  have:  (1)  caused  the  velocity  over  the 
inner  duct  surface  to  be  higher  than  predicted,  and  (2)  reduced 
the  effectiveness  of  the  elevon. 
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Duct  Pressure  Distributions  at  Angle  of  Attack 

The  predicted  pressure  distributions  at  angle  of  attack 
are  in  overall  qualitative  agreement  for  both  ducts  (Figs.  10, 

12,  and  13) .  However,  the  accuracy  is  poorest  on  the  most  criti¬ 
cal  surface  2.  For  surface  2  the  pressure  is  underpredicted  more 
for  the  smaller  than  for  the  larger  duct.  This  is  in  agreement 
with  the  fact  that  the  normal  force  on  the  smaller  duct  was  over¬ 
predicted  whereas  the  normal  force  on  the  larger  model  was  not. 

For  the  smaller  duct  the  inaccuracy  of  prediction  is  ascribed  to 
the  prevelence  of  flow  separation  from  the  duct  and  the  nonuni¬ 
formity  of  the  propeller  loading. 

For  the  larger  duct  at  low  tirust  coefficient  (Fig.  12(a)), 
the  pressure  gradient  on  the  uppermost  surface  4  is  more  adverse 
than  predicted.  This  is  believed  to  be  r  ised  by  the  hook  near 
the  leading  edge  of  the  actual  camberline  (Fig.  3)  which  cannot 
be  accurately  represented  by  the  four-term  series  expression  for 
the  camberline.  The  corresponding  suction  peak  is  much  smaller 
for  the  smaller  duct.  As  the  angle  of  attack  is  increased  to 
40°  in  Figure  12(b)  ,  the  flow  separates  from  surface  4.  The  flow 
is  also  separated  from  surface  2  in  Figure  19(b)  as  shown  by  the 
measured  pressure  distribution  and  the  boundary- layer  calcula¬ 
tions. 

At  higher  thrust  coefficient  (Fig.  13) ,  the  stagnation 
point  moves  onto  surface  4  and  the  pressure  gradient  becomes 
favorable.  The  computed  separation  on  surface  2  was  near  the 
propeller  and  the  separation  bubble  did  not  appear  in  the  data 
until  a  exceeded  the  values  shown  here.  Hence,  the  predicted 
pressure  distributions  are  relative  accurate. 

Duct  Boundary-Layer  Calculations 

In  Figure  14,  the  difference  in  pressure  gradient  up¬ 
stream  of  the  propeller  causes  the  boundary  layer  corresponding 
to  the  measured  pressure  to  separate  ahead  of  the  propeller  while 
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the  predicted  pressure  boundary -layer  does  not.  However,  when 
the  boundary  layer  separates  farther  forward  similar  results  are 
obtained  using  measured  and  predicted  pressure  distributions 
(Fig.  15) . 

Starting  point 

In  Figure  15,  the  boundary  layer  computed  from  the  lead¬ 
ing  edge  is  compared  with  that  computed  from  the  outside  surface 
stagnation  point.  Up  to  the  leading  edge  the  pressure  gradient 
is  highly  favorable  and  the  momentum  thickness  at  the  leading 
edge  is  quite  small.  There  is  only  a  small  difference  in  momen¬ 
tum  thickness  on  the  inside  surface  and  negligible  difference 
(<  0.01c)  between  the  computed  separation  points.  Therefore,  the 
additional  effort  required  to  compute  the  boundary  layer  from 
the  stagnation  point  does  not  seem  required. 

Seven-foot  model 

The  measured  pressure  distribution  on  the  duct  at 
a  =  40°  (Fig.  12(b))  indicates  a  separated  region  at  x/c  *  0.05. 
The  boundary  layer  computed  from  the  leading  edge  using  the  mea¬ 
sured  pressure  distribution  separates  at  x/c  =  0.032.  Since 
this  pressure  distribution  clearly  indicates  a  separation  region, 
the  computed  result  is  encouraging  in  terms  of  the  method  used  to 
predict  the  boundary  layer.  No  separation  was  computed  using  the 
predicted  pressure  distribution,  but  this  distribution  was  in¬ 
accurate  because  of  the  flow  separation. 

Four- foot  model 

A  summary  of  all  the  computed  separation  points  for  the 
smaller  duct  is  shown  in  Figure  17(a).  The  dashed  curves  indicate 
the  location  of  separation  on  the  windward  inner  duct  surface  2 
as  estimated  from  the  computed  points.  These  curves  are  obviously 
rough  estimates  because  of  the  small  number  of  points.  The 
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upper  dashed  curve,  which  can  be  r ,t\ mated  more  accurately  than 
the  others,  gives  the  estimated  ilio’-t  conditions  (angle  of 
attack  versus  advance  ratio)  for  separation  about  midway  between 
the  duct  lip  and  propeller  (x/c  =  0.15).  The  stall  boundary 
deduced  in  Reference  5  from  sound,  tuft,  and  pitching  moment 
data  lies  somewhat  above  the  upper  dashed  curve,  hence,  thir 
stall  boundary  apparently  corresponded  to  separation  very  near 
the  leading  edge  of  the  duct.  The  lowest  dashed  curve  indicates 
that  separation  is  expected  to  occur  before  the  propeller  except 
for  low  angle  of  attack  and  high  advance  ratio.  The  dashed  curve 
for  separation  at  x/c  =  0.24  indicates  that  this  is  the  calcu¬ 
lated  location  of  separation  for  hovering  flight  (any  a  at 
J  =  0)  and  that  there  are  values  of  a  for  J  >  0  which  fix 
the  separation  point  at  x/c  =  0.24.  A  "duct  stall  boundary" 
could  be  defined  to  correspond  with  a  fixed  location  of  separa¬ 
tion  before  the  propeller.  This  boundary  would  then  correspond 
with  any  one  of  the  curves  indicated  in  Figure  12(a)  depending  on 
the  chosen  location  of  separation. 

Boundary-Layer  Separation  in  Axial  Flow 

The  computed  separation  points  for  both  ducts  in  axial 
flow  are  compared  in  Figure  17(b) .  When  the  total  thrust  coeffi¬ 
cient  is  near  20  it  can  be  seen  that  the  boundary  layer  on  the 
larger  duct  separates  slightly  farther  forward,  however,  the 
larger  duct  carries  a  considerably  larger  fraction  of  the  thrust 
(Table  III) .  The  computed  boundary  layers  on  both  ducts  in  axial 
flow  are  separated  ahead  of  the  propeller  when  >  2.  With¬ 

out  unsteady  effects  due  to  finite  propeller  blades,  the  boundary 
layer  was  computed  beyond  the  propeller  for  both  ducts  in  axial 
flow  using  the  measured  velocity  distribution.  Separation  oc¬ 
curred  on  the  smaller  duct  at  x/c  =  0.65  for  Gj*Dp  =  0*684  and 
on  the  larger  duct  at  x/c  =  0.65  for  OTdp  =  0.890.  Thus,  the 
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computed  boundary  layers  on  the  inner  surface  of  both  ducts  sepa¬ 
rate  even  for  thrust  coefficients  smaller  than  unity. 

Boundary-Layer  Separation  Summary 

The  computed  separation  results  for  both  ducts  are  com¬ 
pared  in  Figure  18.  The  solid  line  gives  the  velocity  ratio  for 
separation  as  predicted  for  a  constant  adverse  velocity  gradient 
indicated  in  Figure  18.  The  points  which  lie  above  the  curve  are 
for  increasingly  adverse  velocity  gradients  such  that  u(x)  is 
convex  upward.  The  points  below  the  curve  are  for  u(x)  curves 
which  are  sharply  peaked  near  the  leading  edge  and  convex  down¬ 
ward.  It  can  be  seen  that  all  of  the  computed  separation  poin.,3 

useD 

lie  within  a  rather  narrow  band  of  velocity  ratio,  0.6  <- — e<0.8. 

umax 

By  use  of  Figure  18,  one  can  predict  the  location  of  separation 
for  a  typical  duct  velocity  profile  rather  closely. 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  ability  of  previously  developed  theory  for  a  ducted 
propeller  at  angle  of  attack  to  predict  the  duct-to-propeller 
thrust  ratio,  the  normal  force,  and  the  pitching  moment  has  been 
evaluated  by  use  of  wind  tunnel  data  for  two  large-scale  ducted 
propellers.  The  theoretical  predictions  were  extended  to  include 
the  pressure  distribution,  boundary  layer,  and  stall  boundary  for 
the  duct. 

The  experimental  data  presented  herein  show  that  the 
flow  over  a  ducted  propeller  at  angle  of  attack  is  generally  very 
complex  with  much  free  vorticity  generated  by  non-uniform  blade 
loading  and  separation  of  the  flow  from  the  duct  and  centerbody. 
Specifically,  the  rake  data  for  the  4-foot  model  show  that  the 
blade  loading  was  concentrated  near  the  Y  )  and  that  the  flow 
was  generally  separated  from  the  inner  duct  surface. 

The  theoretical  flow  model  concentrates  the  free  vor¬ 
ticity  into  a  thin  duct  boundary  layer  and  a  single  vortex 
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cylinder  trailing  from  the  duct.  Nevertheless,  this  simple  model 
succeeds  in  predicting  at  least  qualitatively  the  force  and  mo¬ 
ment,  the  pressure  distribution,  and  the  separation  of  the  bound¬ 
ary  layer  over  the  entire  operating  range  of  propeller  thrust  and 
free-stream  angle  of  attack. 

The  following  specific  conclusions  are  drawn  from  the 
experimental  data  and  theoretical  calculations  reported  herein. 

(1)  The  duct  thrust  force  and  pitching  moment,  as  pre¬ 
dicted  for  a  thin  cylinder  surrounding  am  actuator  disk,  corre¬ 
spond  reasonably  well  with  the  measured  data.  In  hovering  flight 
the  4- foot  model  duct  carries  about  50  percent  of  the  total 
thrust  and  the  7-foot  model  duct  about  60  percent.  For  the 
larger  model,  the  duct  thrust  is  predicted  more  accurately  with 
the  assumption  that  the  pressure  rise  across  the  actuator  disk 
acts  on  the  inner  duct  surface  aft  of  the  propeller.  This  is  not 
the  case  for  the  smaller  duct,  apparently  because  the  propeller 
loading  is  low  near  the  blade  tips. 

{2)  The  duct  normal  force  at  angle  of  attack  is  well 
predicted  for  the  7-foot  model  and  considerably  overpredicted  for 
the  4-foot  model.  This  difference  is  also  believed  to  be  caused 
by  the  concentration  of  propeller  loading  nearer  the  hub  for  the 
smaller  model,  since  theoretical  estimates  indicate  that  the  duct 
normal  force  decreases  with  reduced  loading  of  the  propeller 
blade  tips. 

(3)  The  duct  pressure  distributions  predicted  for  axial 
flow  correspond  reasonably  well  with  the  measured  data.  However, 
when  the  thrust  level  and  blade  pitch  are  both  high,  the  pressure 
on  the  inner  surface  of  the  7-foot  duct  is  lower  than  predicted, 
and  the  predicted  discontinuity  in  pressure  across  the  propeller 
is  not  evident  in  the  data.  These  differences  are  believed  to  be 
caused  by  low  loading  of  the  central  part  of  the  propeller  and 
separation  of  the  duct  boundary  layer  as  predicted  before  the 
propeller  except  at  low  thrust. 
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(4)  The  predicted  duct  pressure  distributions  at  angle 
of  attack  are  in  qualitative  agreement  with  the  data,  but  they 
are  not  always  quantitatively  accurate,  particularly  in  the  crit¬ 
ical  region  Tor  flow  separation,  inside  the  windward  duct  lip. 

(5)  For  both  ducts  in  axial  flow,  separation  of  the 
boundary  layer  on  the  inner  surface  is  predicted  to  occur  before 
the  propeller  except  when  the  advance  ratio  is  high  and  the 
thrust  coefficient  is  low.  Separation  is  predicted  slightly  far¬ 
ther  forward  on  the  7-foot  model  duct  than  on  the  4-foot  mode’ 
duct  for  the  same  thrust  coefficient  in  axial  flow,  however,  t  e 
duct-to-propeller  thrust  ratio  is  higher  for  the  larger  duct. 

The  prediction  of  flow  separation  from  the  inner  duct  lip  of  the 
4-foot  model  duct  at  angle  of  attack  corresponds  well  with  the 
stall  boundary  deduced  experimentally  from  tufts,  sound  level, 
and  pitching  moment. 

The  general  conclusion  of  this  study  is  that  the  differ¬ 
ence  between  the  measured  performance  of  the  two  ducted  propellers 
and  the  theoretical  predictions  is  caused  mainly  by  diff  rences 
in  blade  loading  distribution  and  the  prevalence  of  flow  separa¬ 
tion  from  the  inner  duct  surface.  More  data  and  theory  are  needed 
to  define  and  predict  the  distribution  of  blade  loading,  the  duct 
boundary  layer,  and  the  interaction  between  the  propeller  blade 
tips  and  a  region  of  flow  separation. 

It  is  recommended  that  total  pressure  rake  data  be  ob¬ 
tained  to  determine  the  radial  distribution  of  blade  loading 
together  with  data  to  locate  transition  and  separation  of  the 
duct  boundary  layer.  Such  data  could  be  obtained  for  hovering 
flight  quite  readily  without  the  use  of  a  wind  tunnel. 
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EFFECT  OF  PROFILE  THICKNESS  ON 
SURFACE  VELOCITY  DISTRIBUTION 


x/c 

td  i  e 

2V  2tt  COt  2 

t  =  0,  C£  =  1 
(Eq.  (24)) 

'Yd 

2V 

NACA  0018,  c g  —  1 
(p.  325,  Ref.  8) 

us 

■f  -  p(x> 

NACA  0018,  a  =  0 

(p.  325,  Ref.  8) 

0 

00 

1.  342 

0 

0.005 

2.25 

1.178 

0.682 

.025 

.994 

.861 

1.103 

.05 

.695 

.662 

1.228 

.  3. 

.478 

.479 

1.276 

.2 

.  319 

.320 

1.275 

.4 

.  195 

.  184 

1.205 

.  6 

.  130 

.113 

1. 116 

.8 

.079 

.063 

1.025 

1.0 

0 

0 

0 

TABLE  II 


SHAPE  PARAMETERS  FOR  THE  TWO  DUCTED  PROPELLERS 


Parameter 

4-foot 

model 

7-foot 

model 

Dp  (in.) 

48.0 

84.7 

Dh  (in.) 

D^rs.) 

16.0 

17.5 

54.3 

93.3 

c  (in.) 

33.  0 

49.0 

t  (in.) 

5.2 

8.4 

max' 

<**»«« 

0.158 

0.172 

Xp/C 

.293 

.286 

Ap/A 

.70 

.79 

c/D 

.608 

.525 

3.  30 

3.10 

f 

0.  54 

0.53 

2 

f 

1.95 

1.90 

3 

f 

0.93 

0.92 

4 

f5 

.  29 

.22 

f 

1.49 

1.50 

6 

f 

1.92 

1.87 

7 

f  f  +f 

2.  35 

2.20 

5  6  7 

Ro 

-0.007 

-0.040 

R 

-  .007 

-  .068 

l 

R2 

-  .040 

-  .058 

R 

.039 

-  .013 

3 

Ro* 

-  .008 

-  .041 

R  * 

-  .047 

-  .108 

l 

V 

-  .053 

-  .071 

R  * 

.040 

-  .012 

3 

38 


TABLE  III.-  CONCLUDED 


Parameter 
(for  Rn*  =  0) 

4-foot  model 
duct 

1-  - 

C 

o 

0.4922 

c 

l 

.7402 

c 

2 

-  .3142 

C3 

.2107 

C4 

-  .1589 

c 

.  1277 

c 

o 

.955 

c 

1 

-  .280 

c 

2 

-  .027 

c  ,c  .c 

3  4  3 

0 

-0.2487 

.2910 

-  .1625 
.1065 

-  .0796 
.0640 
.1509 

-  .0774 
.0147 

-  .0035 
.0014 

-  .0007 


7-foot  model 
c 


0.5259 

.7270 

-  .3154 
.2107 

-  .1591 
.1280 
.995 

-  .247 

-  .020 

0 


COEFFICIENTS  FOR  4-FOOT  MODEL  DUCT  (EQ.  (19)) 


; i 

3» 

i 

i 


vO 

0.00000 

0 

.000025 

0 

-  .00134 

0 

.00151 

irt 

0 

0. 00008 

0 

-  .00226 

0 

. 00222 

0 

•vj- 

-0. 00002 

0 

-  .00466 

0 

. 00359 

0 

-  .00088 

n 

0 

-0. 01403 

0 

.00690 

0 

-  .00136 

0 

CM 

0.03874 

0 

. 01860 

0 

-  .00234 

0 

.00001 

rH 

0 

0.09140 

0 

-  .00468 

0 

.00002 

0 

o 

0.07745 

. 15490 

.02791 

-  .00007 

-  .00016 

-  .000005 

.00000 

/ 

/  x 

0 

1 

2 

3 

4 

5 

6 

f 


Figure  1 


(a)  4-foot  model. 

. -  Ducted  propellers  mounted  in  the  Ames 
40-  by  80-foot  wind  tunnel. 
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Experiment 


Figure  4.-  Normal  force  versus  thrust  for  the  4-foot  model  at 


Experimental 


Figure  6.-  Normal  force  versus  thrust  for  the  7-foot  model 


Experiment 


Figure  7.-  Pitching  moment  versus  thrust  for  the  7-foot  model  at 
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Figure  9.-  Continued, 
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Figure  9.-  Continued 


Figure  10.-  Continued 
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Figure  15.-  Computed  boundary  layers  for  the  4-foot  model  duct 
at  a  =  80  ,  J  =  0.166. 


(a)  Observed  stall  boundary  and  computed  separation 
points  for  the  4-foot  model  duct. 


(b)  Computed  separation  points  for  both  ducts  at  a  =  0. 


Figure  17.-  Flow  separation  from  both  ducts. 
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AEROTHERMODYNAMIC  PERFORMANCE  OF  A  HIGH 
BYPASS  TIP  TURBINE  CRUISE  FAN  SYSTEM 
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Lift  Fan  Systems  Operation 

Advanced  Technology  and  Demonstrator  Programs  Department 
General  Electric  Company 
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ABSTRACT 

One  of  the  promising  propulsion  systems  for  use  in  V/STOL  aircraft 
systems  is  the  high  bypass  lift/cruise  fan.  This  system  incorporates  the 
turbomachinery  of  the  tip  turbine  driven  lift  fan  presently  powering  the 
XV-5A  research  aircraft.  These  high  bypass  systems  require  inlet  and  exhaust 
system  geometry  outside  the  range  of  most  existing  test  results. 

A  test  program  was  initiated  to  provide  large  scale  model  aero¬ 
dynamic  and  thermodynamic  performance  characteristics  for  this  lift/cruise 
system  installation.  The  model  incorporated  a  propulsion  system  consisting 
of  the  X-376  pitch  fan  from  the  XV-5A  aircraft,  with  a  bypass  ratio  of  ten, 
and  a  T58  core  engine.  Model  geometry  was  obtained  by  changing  the  fan  inlet 
and  exhaust  system.  Testing  was  performed  in  the  NASA- Lang ley  16  foot  wind 
tunnel  covering  a  range  of  zero  to  0.85  Mach. 

This  paper  presents  the  results  of  this  test  program  in  terms  of: 
o  Estimated  remote  fan  and  axisymmetric  nacelle  performance 

o  Performance  of  fan  and  close-coupled  engine-fan  nacelles  including 

interference  problem  caused  by  nacelle  intersections 


o  Possible  solutions  to  problems  of  close-coupled  system  as  applicable 
to  future  test  programs  and  designs  of  full  scale  high  bypass  lift/ 
cruise  fan  systems. 
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inlet  leading  edge  area 

stream  tube  area  for  inlet  mass  flow 

nacelle  area  based  on  maximum  diameter 


nacelle  drag  coefficient, 


N 


A 

40  TT 


pressure  coefficient 
nacelle  external  drag 

nacelle  leading  edge  diameter 
inlet  throat  diameter 
nacelle  maximum  diameter 

inlet  loss  correction  factor  for  mass  flow  ratio 
inlet  length  from  maximum  diameter  to  leading  edge 

nacelle  length 

total  installation  length 

Mach  number  for  local  critical  flow 

free- 8 t ream  Mach  number 

free-stream  dynamic  pressure 

average  velocity  at  inlet  throat  based  on  total  pressure  and 
flow 
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i  -  fan  flow  coefficient,  ratio  of  inlet  velocity  to  fan  tip  speed 

Y  -  fan  pressure  coefficient,  ratio  of  fan  pressure  rise  to  dynamic 

pressure  based  on  fan  tip  speed 

‘tyg  -  fan  stage  efficiency,  ratio  of  fan  horsepower  developed  to  turbine 
power  extracted 

ui  -  inlet  loss  coefficient  based  on  throat  velocity  head 

INTRODUCTION 

Recent  incrr.’.sed  emphasis  on  V/STOL  aircraft  systems  has  shown  the 
need  of  several  new  families  of  propulsion  systems.  One  of  the  promising 
systems  in  the  high  bypass  category  is  the  lift/cruise  fan.  This  propulsion 
system  incorporates  the  principle  of  the  tip  turbine  driven  fan.  The  main 
propulsion  system  component  is  the  fan  rotor,  consisting  of  a  single  stage 
fan  with  the  turbine  blades  attached  around  the  periphery  of  the  fan  blades. 
This  component  is  an  outgrowth  of  the  technology  of  the  fan  systems  presently 
used  in  the  XV-5A  research  aircraft.  Figure  1  shows  the  present  wing  fan 
rotor  system  used  in  the  XV-5A  aircraft. 

In  addition  to  the  rotor,  the  complete  propulsion  package  consists 
of  the  inlet  and  exhaust  systems  for  the  fan,  and  the  manifolds  required  for 
ducting  the  core  engine  exhaust  gases  to  the  fan  tip  turbine.  Figure  2  shows 
a  typical  lift/cruise  fan  system.  The  fan  inlet  consists  of  the  bulletnose 
and  nacelle  front  cowl.  The  cowl  is  designed  to  meet  both  static  takeoff  and 
cruise  requirements  as  well  as  providing  an  enclosure  for  the  hot  gas  ducting. 
The  installation  as  shown  in  a  close-coupled  arrangement  where  the  gas 
generator  or  core  engine  nacelle  is  close  to  or  part  of  the  fan  inlet 
system.  Another  typical  installation  is  the  remote  arrangement  that  per¬ 
mits  locating  the  core  engine  in  any  convenient  position  and  ducting  of  the 
hot  gas  flow  to  the  fan  nacelle  through  a  pylon  mount. 

The  exhaust  system  for  the  fan  forms  the  nacelle  afterbody  and 
provides  a  means  of  area  control  for  the  fan  system.  The  typical  arrangement 
as  shown  uses  a  convergent  plug  nozzle  with  an  inflatable  rubberized  fabric 
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boot  on  the  plug,  providing  a  two  position  nozzle.  Both  the  hot  turbine  and 
cold  £an  flows  are  discharged  through  the  common  nozzle  annulus;  the  cold 
flow  being  adjacen*.  to  the  nozzle  plug  system. 

Design  studies  have  shown  that  a  typical  optimum  lift/cruise  fan 
system  will  have  a  bypass  ratio  of  between  5  and  8  and  a  fan  pressure  ratio 
of  1.3  to  1.5.  At  present  there  is  no  operating  full  scale  fan  hardware  in 
this  range  of  pressure  and  bypass  ratios,  one  of  the  closest  systems  is 
being  the  pitch  control  fan  of  the  XV-5A.  This  fan  system,  designated 
X-376,  has  a  bypass  ratio  of  about  10,  a  pressure  ratio  of  1.1,  and  a  fan 
tip  diameter  of  36  inches.  In  order  to  obtain  much  needed  installation 
aerodynamics  applicable  to  the  lift/cruise  fan  systems,  a  test  program 
utilizing  the  X-376  fan  system  was  initiated.  The  tests  were  performed 
in  the  NASA- Lang ley  16  foot  transonic  wind  tunnel  under  the  sponsorship 
of  the  U.S.  Army  Aviation  Materials  Laboratory. 

This  paper  presents  some  of  the  most  important  performance  criteria 
obtained  during  conduction  of  the  test  program.  The  results  of  the  test 
program  are  reported  fully  in  Reference  1. 

THE  TEST  MODEL 

The  model  (Figure  3)  used  in  the  test  program  consisted  of  the 
following  major  components  arranged  as  a  close-coupled  lift/cruise  fan 
system: 

o  The  X-376  fan  system 
o  A  T-58  core  engine 
o  The  engine  inlet  nacelle 
o  Three  interchangeable  fan  inlet  cowls 
o  Three  interchangeable  afterbody  systems 
o  Three  interchangeable  nozzle  plug  systems. 

The  combination  of  interchangeable  components  provided  for  twenty- 
seven  possible  tes"  configurations.  Eight  of  the  most  realistic  configura¬ 
tions  were  tested.  Pertinent  geometric  data  for  the  eight  test  configurations 


4 


are  given  in  Figure  4.  The  test  models  covered  a  range  of  design  cruise 
Mach  numbers  from  0.55  to  0.8,  with  variations  of  both  nozzle  area  and 
afterbody  length. 

The  test  model  was  sting  mounted  in  the  wind  tunnel  (Figure  5), 
with  the  engine  inlet  located  below  the  fan  inlet  centerline.  This  arrange¬ 
ment  is  inverted  from  the  normal  expected  installation  as  an  aircraft 
propulsion  system,  but  was  required  to  facilitate  installation  in  the 
wind  tunnel. 

Instrumentation  was  provided  for  measurement  of  pressure  distri¬ 
butions  on  both  the  internal  and  external  surfaces  of  the  model.  In  addi¬ 
tion,  the  flow  conditions  were  measured  at  the  inlet  and  discharge  planes 
of  the  core  engine  and  fan  system.  This  instrumentation  provided  data  for 
evaluating  propulsion  system  internal  and  nacelle  external  performance. 

Total  model  forces  were  measured  by  a  three  component  strain  gage  balance 
located  at  the  attachment  of  the  model  to  the  sting  mount. 

The  test  program  included  investigation  of  system  performance  for 
a  range  of  Mach  numbers  up  to  0.85  at  various  power  settings  of  the  propul¬ 
sion  system.  Performance  under  conditions  of  variable  angle  of  attack  were 
investigated  on  one  of  the  test  models,  but  will  not  be  discussed  in  this 
presentation. 

PRESENTATION  OF  RESULTS 

Since  the  model  used  in  this  test  program  included  a  complete 
propulsion  system  in  addition  to  an  aerodynamic  model,  the  test  results  will 
be  presented  as  two  different  categories.  The  first  category  will  involve 
the  aerothermodynamic  performance  of  the  fan  system  under  conditions  for 
forward  speed.  Normal  propulsion  system  parameters  will  be  used  such  as 
thermodynamic  efficiency  and  fan  pressure  ratio.  The  second  category  and 
the  main  body  of  the  results  will  involve  the  aerodynamic  performance  of  the 
fan  inlet  and  afterbody  systems.  The  performance  in  this  case  will  be 
presented  as  local  surface  velocity  distributions  using  pressure  coefficient 
notation.  Particular  attention  will  be  directed  to  the  levels  of  the  maximum 
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surface  velocities  and  the  conditions  required  for  establishment  of  local 
velocities  greater  than  critical  (Mach  number  equal  to  1).  The  conditions 
required  for  local  critical  flow  will  be  discussed  and  compared  with  similar 
existing  test  results. 

Test  Results  -  Propulsion  System 

The  X-376  fan  system  is  presently  installed  as  the  pitch  control  fan 
of  the  XV-5A  aircraft.  In  this  installation, the  fan  is  installed  with  its 
axis  normal  to  the  free-stream  velocity.  Fan  performance  has  been  investi¬ 
gated  for  a  range  of  flight  speeds  up  to  100  knots.  The  Langley  installation 
as  tested,  provided  the  first  source  of  test  data  for  the  fan  installed  with 
its  axis  aligned  with  the  flow  direction  at  speeds  up  to  a  Mach  number  of 
0.8S.  A  second  installation  difference  existed  in  that  the  fan  exhaust 
system  incorporated  a  confluent  or  mixed  flow  nozzle  system.  This  type  of 
nozzle  is  characterized  by  the  merging  of  both  the  fan  turbine  and  fan 
discharge  streams  prior  to  accelerating  the  mixed  flow  in  the  nozzle  throat. 
In  this  system,  there  is  an  interaction  of  the  two  flow  streams  that  is  of 
interest  in  terms  of  effects  on  propulsion  system  performance. 

Figure  6  shows  a  typical  fan  operating  map  obtained  while  testing 
one  of  the  eight  test  models.  Fan  operation  in  this  case  is  depicted  as  the 
variation  of  stage  efficiency  and  pressure  coefficient  with  fan  flow  coeffi¬ 
cient.  These  are  conventional  parameters  presently  used  in  defining  perfor¬ 
mance  of  a  single  stage  fan  system  with  relatively  low  pressure  ratios.  The 
data  shown  in  the  figure  represents  the  complete  test  range  from  static  to 
a  Mach  number  of  0.8.  The  complete  range  of  core  engine  power  settings  from 
"windmill"  to  "maximum"  is  also  included  in  the  data.  Windmill  is  the  con¬ 
dition  with  the  core  engine  inoperative  and  fan  rotational  speed  provided  by 
the  ram  drag  forces.  The  significance  of  this  single  performance  charac¬ 
teristic  is  the  convergence  of  all  data  into  a  single  characteristic,  that 
also  agrees  quite  well  with  the  predicted  performance.  This  effect  lends  to 
confidence  in  predicting  cruise  fan  performance  based  on  static  or  predicted 
performance.  An  almost  identical  characteristic  was  obtained  for  the  seven 
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other  test  models.  Indicating  that  reasonable  changes  in  nacelle  geometry 
do  not  produce  significant  changes  in  propulsion  system  performance. 

The  second  important  propulsion  system  performance  parameter  that 
was  verified  by  this  program,  was  the  performance  of  the  exhaust  nozzle 
system.  In  analysis  of  mixed  flow  cruise  fan  system,  it  has  been  normal  to 
assume  the  static  pressure  of  the  fan  and  fan  turbine  flows  to  be  equal  at 
the  plane  where  mixing  is  initiated.  Measurements  of  the  pressures  in  this 
plane  of  the  model,  verified  this  assumption  within  reasonable  accuracy. 

All  eight  test  models,  covering  a  range  of  exhaust  nozzle  geometry  showed 
that  these  pressure  levels  were  equal. 

Utilizing  the  fan  map  presented  above,  the  verification  of  equal 
static  pressure  at  the  discharge  plane,  and  an  assumption  of  an  ideal  inlet 
and  exhaust  system,  the  performance  of  the  X-376  fan  system  as  a  lift/ 
cruise  fan  can  be  obtained  readily  and  is  shown  in  Figure  7.  Here,  fan 
net  thrust,  power  requirements  and  optimum  nozzle  area  are  presented 
versus  flight  Mach  number.  Two  levels  of  nacelle  drag  coefficient  are 
shown.  As  can  be  seen,  nacelle  drag  has  significant  effects  on  system 
thrust  for  a  low  pressure  ratio  fan  system.  Effects  of  nacelle  drag  are 
less  for  a  higher  pressure  ratio  system. 

Test  Results  -  Inlet  System 

Some  of  the  most  significant  results  obtained  during  the  test  pro¬ 
gram  involved  the  aerodynamics  of  the  numerous  nacelle  geometries  in  a 
typical  close-coupled  lift/cruise  fan  installation.  Representative  remote 
lift/cruise  fan  installation  performance,  may  be  assumed  to  exist  on  the 
axisymmetric  part  of  the  nacelle  opposite  the  core  engine  nacelle.  The 
following  is  a  discussion  of  some  of  the  inlet  performance  obtained  by 
analysis  of  the  test  results. 

Fan  inlet  internal  performance  may  be  presented  in  two  ways;  first, 
the  effects  of  inlet  geometry  on  inlet  total  pressure  losses  and,  second,  the 
relationships  of  local  maximum  surface  velocities  on  the  inlet  loading  edge. 
The  three  inlet  geometries  tested  are  shown  in  Figure  8  for  reference. 
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Figure  9  presents  typical  inlet  loss  coefficient  variation  with  mass 
flow  ratio  for  the  three  inlet  systems.  The  fairing  of  the  actual  test 
results  is  shown  as  a  solid  line.  The  dashed  curves  represent  an  empirical 
correlation  of  the  test  data  using  the  following  equation: 

_  ^ 
ui  =  0.017  ——  (K) 

UTH 

The  first  part  of  the  equation  is  an  empirical  factor  relating  inlet 
loss  with  inlet  length  and  diameter  at  static  conditions.  The  second  term  is 
an  empirical  constant  that  is  a  function  of  mass  flow  ratio  as  shown  in 
Figure  10.  The  agreement  of  the  empirical  equation  and  the  actual  test  data 
is  very  good.  This  correlation  can  be  a  useful  tool  in  predicting  perfor¬ 
mance  of  lift/cruise  fan  system  where  inlet  losses  have  significant  effects 
on  total  installed  performance. 

The  second  indicator  of  inlet  performance  is  the  maximum  surface 
velocity  near  the  inlet  leading  edge.  Figure  11  shows  a  typical  family  of 
inlet  internal  pressure  distributions  for  one  of  the  test  inlets  covering 
a  range  of  test  Mach  numbers  and  mass  flow  ratios.  Correlation  of  the  peak 
pressure  coefficient  with  mass  flow  ratio  is  shown  in  Figure  12.  The  two 
longer  inlet  systems,  inlet  1  and  2,  exhibit  similar  characteristics.  Inlet 
3  shows  significantly  higher  velocities  at  the  higher  speed  condition,  high 
values  of  ^1/Aq.  This  difference  may  be  attributed  to  the  extension  of  the 
bulletnose  foreward  of  the  inlet  leading  edge.  The  data  shows  that  an  inlet 
system  of  the  type  tested,  all  having  an  inlet  leading  edge  diameter  of  1.13 
times  the  throat  diameter,  will  have  a  peak  local  velocity  of  about  1.6  times 
the  calculated  average  velocity.  Using  local  sonic  flow  as  a  criteria,  this 
limits  the  inlet  mass  flow  at  static  condition  to  84  percent  of  the  theoret¬ 
ical  choking  value. 

Tc;.,ts  under  conditions  of  variable  angle  of  attack  showed  that  the 
nacelle  angle  of  attack  must  be  greater  than  12  degrees  before  the  peak  inlet 
internal  velocity  exceeds  the  levels  as  occur  at  static  conditions.  This 
tends  to  indicate  that  the  design  conditions  for  the  inlet  throat  size  is 
established  by  the  static  operating  condition. 
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Nacelle  external  performance  may  be  Investigated  by  observing  the 
levels  of  maximum  velocity  on  the  inlet  external  surface  under  conditions  of 
forward  speed  and  varying  mass  flow.  Figure  13  shows  a  typical  set  of 
nacelle  pressure  distributions  taken  at  low  Mach  numbers  and  at  different 
mass  flow  ratios.  The  pressure  distributions  are  typical  of  a  low  speed 
inlet  of  this  type  and  exhibit  highest  local  velocities  under  conditions 
of  minimum  mass  flow  ratio.  It  should  also  be  pointed  out  that  the  levels 
of  pressure  coefficient  around  the  nacelle  system  vary  considerably;  minimum 
pressure  coefficients  occur  in  the  region  of  the  fairing  of  the  fan  inlet 
to  the  engine  nacelle. 

Cross-ploting  the  numerous  pressure  distributions  for  the  eight 
test  models  for  a  range  of  mass  flow  ratios  and  ftee-stream  Mach  numbers 
yielded  the  following  conclusions: 

o  The  peak  velocities  on  the  fan  inlet  external  surface  are  Independent 
of  afterbody  geometry  for  the  range  of  models  tested. 

o  Peak  pressure  coefficient  levels  are  generally  higher  in  the  region 
of  the  fairing  and  approach  critical  flow  condition  at  Mach  numbers 
lower  than  expected  based  on  axisymmetric  inlet  design  criteria. 

o  Characteristic  increases  in  peak  pressure  coefficient  were  observed 
with  conditions  of  decreasing  mass  flow  ratio  and  increasing  Mach 
number. 

Faired  test  data,  showing  effects  of  free-stream  Mach  number  and 
inlet  circumferential  position,  are  presented  in  Figure  14,  for  the  three 
inlet  geometries.  The  rise  in  pressure  coefficient  with  Mach  number  agrees 
quite  well  with  typical  theories  of  compressibility.  The  variation  of 
critical  Mach  number  with  circumferential  position  is  also  apparent,  showing 
interaction  effects  at  the  fairing  section  of  the  inlet. 

Effect  of  mass  flow  on  local  peak  pressure  coefficient  are  also 
shown  in  Figure  15. 
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Low  Mach  number  conditions  are  shown  in  these  figures  because 
restriction  in  propulsion  system  operation  did  not  permit  large  excursions 
of  mass  flow  at  the  higher  Mach  numbers.  This  restriction  existed  because 
at  high  speeds  fen  rotational  speed  at  windmill  conditions  and  maximum  power 
conditions  were  equal,  and  consequently  mass  flow  did  not  vary. 

Making  use  of  this  data  for  relatively  short  inlet  systems  and  the 
previous  correlations  performed  in  Reference  2,  a  set  of  design  character is- 
tics  for  selection  of  inlet  geometry  based  on  local  critical  flow  were 
obtained.  These  characteristics  are  shown  in  Figure  16  and  17.  Using  the 
same  criteria  as  in  Reference  2,  Figure  16  shows  the  variation  of  critical 
Mach  number  with  both  inlet  length  and  inlet  leading  edge  diameter.  In 
general,  increases  in  inlet  length  and  increases  in  inlet  diameter  result 
in  higher  Mach  number  capability  for  the  inlet.  However,  mass  flow  ratio 
must  be  greater  than  that  shown  in  Figure  17  in  order  for  the  critical  Mach 
number  characteristics  to  apply.  This  figure  shows  that  increases  in  inlet 
diameter  cause  additional  restrictions  in  mass  flow  ratios  at  which  the 
inlet  can  operate.  An  iterative  procedure  of  designing  an  inlet  system, 
taking  into  consideration  inlet  mass  flow  ratio,  design  Mach  number  and 
geometric  limitations,  is  required. 

FAIRING  RECONTOUR  BASED  ON  TEST  RESULTS 

As  previously  pointed  out  in  the  discussion,  the  region  in  the 
vicinity  of  the  fairing  exhibits  large  negative  pressure  coefficients  and 
consequently  low  critical  Mach  numbers.  Making  use  of  the  test  data  of  Figure 
15, 16  and  17,  it  is  possible  to  develop  a  criteria  for  redesign  of  the  fairing 
region  for  nearly  uniform  circumferential  velocity  distributions.  Comparison 
of  conditions  for  equal  pressure  coefficients,  tends  to  indicate  that  the 
fairing  region  is  operating  with  a  mass  flow  of  about  1.4  times  the  equivalent 
axisymmetic  condition.  This  infers  that  the  inlet  leading  edge  in  the  fairing 
is  spilling  nit  only  the  normal  amount  of  flow,  but  an  additional  40  percent 
that  is  not  spilled  over  the  inlet  in  the  vicinity  of  the  engine  nacelle. 

This  level  seems  quite  reasonable.  Making  use  oi  this  additional  spillage  or 
reduced  mass  flow  ratio,  it  is  possible  to  determine  a  set  of  characteristics 


applicable  to  redesign  of  the  inlet  contour  in  the  engine  to  fan  nacelle 
fairing  region.  Figure  18  shows  the  proposed  design  parameters  applicable 
to  the  fairing  design.  The  factors  presented  are  an  increase  in  inlet  length 
and  a  reduction  in  inlet  diameter  as  required  for  simultaneous  critical  flow 
over  the  complete  nacelle  system.  This  set  of  characteristics  were  obtained 
by  simply  recalculating  the  design  criteria  as  presented  for  axisymmetic 
inlet  systems,  (Figure  16  and  17)  at  a  40  percent  reduced  mass  flow  ratio. 
Comparison  of  required  inlet  length  and  inlet  leading  edge  diameter  then 
resulted  in  the  characteristics  as  presented. 

A  comparison  of  a  typical  cruise  fan  nacelle  system  in  the  original 
test  and  the  proposed  modified  configurations  is  shown  in  Figure  19.  It  is 
apparent  from  the  sketch  that  the  redesign  produces  an  inlet  that  has  both 
a  swept  back  effect  and  a  general  thickening  of  the  inlet  lip  in  the  vicinity 
of  the  fairing.  The  sweep  back  angle  would  be  about  10  degrees  for  an  inlet 
system  designed  for  a  critical  Mach  number  of  about  0.8. 

An  interesting  verification  of  this  redesign  criteria  was  obtained 
while  testing  one  of  the  models  under  conditions  of  variable  angle  of  attack. 
Tests  on  inlet  1,  showed  that  general  decreases  in  local  pressure  coefficient 
and  consequently  local  Mach  numbers  occurred  in  the  fairing  region  under 
conditions  of  negative  angle  of  attack.  Negative  angle  of  attack  exists 
when  the  engine  inlet  moves  down  in  front  of  the  fan  inlet.  A  negative 
angle  of  attack  of  about  6  degrees  was  required  to  produce  pressure  distri- 
butions  equivalent  to  design  Mach  number  in  the  fairing  region.  This  level 
of  6  degrees  compares  quite  favorably  with  the  10  degrees  discussed  a' eve. 

NACLIIE  EXTERNAL  DRAG 

A  region  of  lift/cruise  fan  performance  lacking  either  experimental 
and  analytical  data  is  total  nacelle  external  drag.  This  test  program  pro¬ 
vided  a  source  of  this  type  of  experimental  test  data.  Utilizing  total  model 
forces,  as  measured  by  the  balance  system,  and  fan  system  net  thrust,  based 
on  measured  flow  conditions,  the  nacelle  external  drag  was  evaluated.  Only 
one  questionable  parameter  is  required  in  determining  the  nacelle  drag  and 
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this  is  the  exhzust  nozzle  thrust  coefficient.  For  evaluetion  of  the  test 
date,  the  nozzle  thrust  coefficient  was  estimated  using  the  following  two 
asswptions: 

o  The  nozzle  velocity  coefficient,  because  of  flow  angularity,  will 
be  0.99. 

o  All  other  nozzle  losses  will  be  due  to  surface  skin  friction  that 
was  calculated  based  on  measured  surface  velocities  and  an  average 
skin  friction  coefficient. 

These  two  assumptions  resulted  in  nozzle  thrust  coefficients  between 
0.95  and  0.97.  The  low  levels  of  thrust  coefficient  are  due  to  the  scrubbing 
drag  on  the  large  areas  of  the  nozzle  plug  system.  It  should  be  noted  that 
the  large  plug  geometries  are  not  typical  of  the  1.3  to  1.5  pressure  ratio 
fan  system  but  are  required  because  the  X-376  is  only  a  1.1  pressure  ratio 
and  was  used  in  a  typical  higher  pressure  ratio  nacelle  system. 

Making  use  of  these  estimated  thrust  coefficients  and  the  measured 
balance  forces,  the  nacelle  total  low-speed  drag  as  measured  is  shown  in 
Figure  20.  This  nacelle  drag  coefficient  included  the  following: 

o  Fan  inlet  external  friction  drag 
o  Afterbody  friction  drag 

o  Complete  nacelle  form  drag  including  the  fairing  region 
o  Core  engine  inlet  friction  and  form  drag. 

Not  included  in  the  drag  coefficient  are: 
o  Inlet  internal  losses 

o  Exhaust  nozzle  losses  including  plug  friction  drag. 

Measured  system  drag  coefficients  varied  between  0.09  and  0.04  depend¬ 
ing  on  mass  flow  ratio  and  inlet  geometry.  Estimated  engine  nacelle  friction 
drag  coefficient  was  estimated  to  be  0.016  and  is  shown  on  the  figure.  Since 
engine  nacelle  form  drag  is  probably  a  small  fraction  of  the  total  drag  as 
shown,  the  drag  numbers  as  measured  and  reduced  by  0.016,  represent  a  typical 
total  nacelle  drag  coefficient  level.  The  top  band  of  the  data  represents  the 
longer  nacelle  systems  and  the  lower  region,  the  shorter  configurations. 
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The  dreg  coefficients  as  previously  discussed  represent  the  low- 
speed  values  only  in  the  range  of  free-stream  Mach  numbers  between  0.2  and 
0.4.  A  premature  drag  rise  was  also  indicated  by  the  balance  measurements, 
similar  in  nature  to  the  premature  critical  flow  problems  in  the  fairing 
region.  Using  faired  data  for  the  eight  test  models,  it  was  observed  that 
the  drag  rise  with  increasing  Mach  number  was  a  function  of  the  nacelle  inlet 
configuration  only.  The  ratio  of  drag  to  low  speed  drag  for  the  three  inlets 
is  shown  in  Figure  21  as  a  function  of  free-stream  Mach  number.  Note  that 
drag  rise  occurred  at  Mach  numbers  of  0.4,  0.5  and  0.6  for  inlets  3,  2  and 
1,  respectively.  Design  Mach  numbers  for  the  inlet  were  0.55,  0.7  and  0.8. 
Although  the  exact  cause  of  the  4rag  rise  is  not  readily  apparent,  it  is 
probable  that  the  problems  of  the  fairing  design  and  associated  interactions 
are  the  cause  of  these  trends.  This  is  a  region  of  performance  that  will 
remain  unanswered  until  a  test  program  of  the  redesign  system  is  accom- 
polished . 

CONCLUSIONS 

Analysis  of  high  speed  tests  of  a  family  of  typical  lift/cruise  fan 
propulsion  systems  indicated  the  following: 

1.  The  lift/cruise  fan  system  is  not  adversely  affected  by  simulated  high 
speed  operation. 

2.  Nacelle  low  speed  drag  coefficients,  not  including  core  engine  nacelle 
drag  ranged  between  about  0.035  at  high  mass  flow  ratios  and  0.065  at 
near  design  mass  flow  ratios. 

3.  The  models  as  tested  exhibited  high  spillage  and  consequently  low 
critical  Mach  numbers  in  the  region  of  the  core  engine  to  fan  nacelle 
fairing.  A  proposed  redesign  of  the  fairing  was  established  using  the 
test  data. 
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Figure  3.  Cross  Section  of  Test  Model. 
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Figure  11.  Typical  Inlet  Internal  Pressure  Distributions 
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Figure  14.  Effects  of  Mach  Number  on  Inlet  External 
Pressure  Coefficient. 
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Figure  15-  Effects  of  Mass  Flow  Ratio  on  Inlet  External 
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Figure  21.  Nacelle  Drag  Rise. 
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THRUST  DEFLECTION  NOZZLES  FOR  VTOL  AIRCRAFT 

by 

VINCENT  J.  DISABATO 
Pratt  and  Whitney  Aircraft 
East  Hartford,  Connecticut,  United  States 


INTRODUCTION 

Airplane  studies  reported  in  the  literature  have  indicated  that  an  aircraft 
using  both  lift  cruise  engines  and  lift  engines  provides  a  good  propulsion  system  for 
a  VTOL  aircraft.  The  lift  cruise  engines  are  sized  for  horizontal  flight.  For  VTOL 
thrust  from  the  lift  cruise  engine  is  deflected  downward,  and  the  resulting  thrust  is 
supplemented  by  separate  lift  engines  to  provide  the  necessary  total  thrust.  This 
paper  discusses  various  means  of  deflecting  the  thrust  from  lift  cruise  engines  in  an 
efficient  manner.  The  thrust  deflector  should  have  a  vertical  thrust  to  deflector 
weight  higher  than  the  thrust  to  weight  ratio  of  the  lift  engines. 

The  thrust  deflection  device  should: 

1.  Be  lighter  than  a  corresponding  lift  engine. 

2.  Deflect  the  thrust  as  efficiently  as  possible. 

3.  Provide  horizontal  performance  equal  to  or  nearly  equal  to  an  undeflected 
nozzle. 

Data  from  nozzle  model  tests  and  thrust  deflection  devices  conducted  in  the 
early  1960's  is  presented.  This  data  illustrates  the  performance  trends  that  can  be 
expected  from  various  thrust  deflection  devices  and  engine  nozzles. 

GENERAL  THRUST  DEFLECTOR  NOZZLE  DESIGN 

The  wide  range  of  flight  capability  required  of  future  VTOL  aircraft  demands 
nozzles  that  will  be  extremely  versatile  to  meet  critical  design  points  efficiently  and 
still  maintain  a 'deflection  capability.  Technology  in  the  area  of  nozzles  for  horizontal 


subsonic  and  supersonic  flight  has  been  developing  steadily  with  the  aid  of  scale  model 
and  flight  testing.  Thrust  deflection,  however,  is  relatively  new  and  adds  more  com¬ 
plexity. 

The  function  of  the  thrust  deflection  system  is  to  turn  the  engine  airflow  from 
an  axial  direction  to  a  vectored  or  vertical  direction  and  thus  provide  lift.  The 
vertical  thrust  produced  by  the  deflector  depends  upon  the  static  pressures  acting  on 
the  walls  of  the  nozzle  and  the  wall  static  pressure  must  be  higher  than  the  static 
pressure  surrounding  the  nozzle  to  produce  a  desirable  thrust  force. 

Figure  1  shows  the  regions  of  pressure  change  on  the  inside  and  outside  wall, 
(from  Ref.  3,  Figure  73).  The  static  pressure  distributions  in  a  curved  channel,  up¬ 
stream  of  a  nozzle,  differs  from  that  of  a  straight  duct  due  to  the  centrifugal  force 
associated  with  the  turning.  This  radial  force  produces  an  increase  in  pressure  on  the 
outside  wall.  On  the  inside  wall,  the  pressure  decreases  around  the  bend  and  then 
increases  to  the  exit. 

The  regions  of  adverse  pressure  gradients  can  cause  the  flow  to  detach  or 
separate  with  a  loss  in  total  energy.  The  flow  separation  is  most  severe  on  the 
inner  wall  where  the  highest  adverse  pressure  gradient  occurs.  This  is  shown 
graphically  in  Figure  2,  (from  Ref.  1,  Figure  11).  The  velocity  distribution  is  shown, 
at  the  exit  of  the  ber.d,  as  a  function  of  the  distance  from  the  inner  wait  for  inlet 
Mach  numbers  of  0.32  and  0.6b.  It  can  be  seen  from  the  figure  that  the  flow  is 
separated,  and  no  flow  is  encountered  for  7  to  12  percent  of  the  duct. 

The  totai  pressure  loss  distribution  for  an  inlet  Mach  number  of  0.2  is  shown 
in  Figure  3,  (from  Ref.  2,  Figure  4).  The  inner  wall  shows  the  largest  contribution 
to  total  pressure  loss;  this  is  true  for  the  inlet  Mach  numbers  of  interest  from  0.10 
to  0.60. 

dmv2 

Since  the  centrifugal  force  is  a  function  of  the  turning  radius,  F  - - 

the  parameter  r/d  is  important,  where  "r"  is  the  radius  of  the  centerline  of  the 
constant  area  duct  and  "d"  is  the  diameter.  The  turning  radius  ratio,  pressure 
loss  variation  is  shown  in  Figure  4,  (from  Ref.  1,  Figure  10),  for  a  few  inlet  Mach 
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numbers  The  total  pressure  loss  through  the  constant  area  duct  is  now  defined. 

It  can  be  assumed  that  the  addition  of  a  nozzle  to  the  duct  will  not  appreciably 
change  the  pressure  loss  through  the  duct  for  an  engine  operating  pressure  ratio 
(Pt/Pa)  of  approximately  two.  The  total  pressure  loss  through  the  duct  results  in  a 
loss  in  gross  thrust  coefficient.  The  total  pressure  loss  is  converted  to  gross  thrust 
coefficient  loss  by  the  following  equation. 
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The  gross  thrust  coefficient  loss  is  shown  in  Figure  5  as  a  function  of  the 
duct  turning  radius  for  a  few  inlet  Mach  numbers. 

For  a  constant  area  duct  turn,  the  following  compromises  are  evident  from 
a  study  of  Figure  5: 

1.  Low  inlet  Mach  number  is  necessary  for  a  high  gross  thrust  coefficient. 

2.  A  turning  radius  of  approximately  2. 5  duct  diameters  is  optimum  for 
high  gross  thrust  coefficient. 

GENERAL  THRUST  DEFLECTOR  REQUIREMENTS 

The  cycle  used  in  the  lift-cruise  engine  has  a  strong  influence  on  the  design 
of  the  thrust  deflector.  The  number  of  deflectors  for  each  engine  can  vary  from  one 
or  two  for  a  simple  jet  engine  to  as  many  as  four  for  a  fan  engine  (two  for  the  fan 
and  two  for  the  primary  stream).  The  burner  length  requirements  for  the  augmenting 
system  can  affect  the  location  of  the  deflectors .  Augmentation  may  also  cause  rejection 
of  otherwise  attractive  deflection  schemes  because  of: 

1.  Excessive  cooling  requirements. 

2.  Variable  area  nozzle  requirements. 

The  engine  installation  in  the  aircraft  affects  the  design  of  the  deflector 
system  in  the  following  ways: 
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VTOL  mode 


The  basic  stability  of  the  aircraft  during  take-off  and  landing  requires  that 
the  sum  of  the  moments  produced  by  the  thrust  of  the  engines  around  the  aircraft 
c.g.  be  in  static  balance.  Using  lift  engines  in  addition  to  the  deflected  lift  cruise 
engine  allows  greater  freedom  in  the  choice  of  engine  and/or  engine  deflection 
locations,  i.e.  the  engines  need  not  be  concentrated  at  the  aircraft  c.g. 

Transition  mode 

The  deflection  device  must  be  continuously  controllable  during  transition  to 
provide  an  axial  component  to  accelerate  the  aircraft  while  maintaining  enough  lift. 

No  unwanted  roll  components  should  exist.  The  deflector  should  respond  quickly  and 
accurately  during  transition. 

Forward  flight  mode 

The  nozzle  must  have  high  performance  in  the  horizontal  mode  where  all  of 
the  usual  nozzle  thrust  and  drag  problems  with  external  flow  will  be  present. 

DISCUSSION  OF  DEFLECTOR  TYPES 

Deflection  devices  which  have  been  studied  in  the  past  few  year?  include: 

1.  Swivel  nozzle 

2.  Ventral  nozzle 

3.  Ventral  nozzle  with  rotating  cascade 

4.  Aft-hood  deflector 

Swivel  Nozzle 

The  swivel  nozzle  uses  the  same  duct  for  the  horizontal  and  vertical  flight 
modes.  The  nozzle  is  rotated  on  a  bearing  to  change  the  thrust  vector  from  horizontal 
to  vertical  and  is  shown  in  Figure  6.  The  gas  path  for  the  horizontal  mode  requires  a 
degree  of  flow  turning.  T1  :c  inherently  causes  a  slight  performance  loss  in  the 
horizontal  mode.  As  the  nozzle  is  rotated  for  transition,  the  Tow  turning  is  increased 
by  a  small  amount  and,  therefore,  the  transition  performance  ard  the  vertical  per¬ 
formance  can  be  good.  During  transition,  however,  the  swivel  nozzle  develops  a 
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side  force  as  it  rotates  around  the  bearing.  This  side  force  must  be  counterbalanced 
if  only  one  lift  cruise  engine  is  planned  for  the  aircraft,  usually  requiring  the  use  of 
two  swivel  nozzles. 

The  swivel  nozzle  is  simple  and  can  have  a  simple  and  reliable  actuation 
system.  The  thrust  to  weight  of  this  type  of  system  is  in  the  range  of  35  to  50 
according  to  our  studies.  The  swivel  nozzle,  however,  is  not  flexible  enough  to  be 
used  or  modified  effectively  for  efficient  supersonic  horizontal  flight.  The  swivel 
nozzle  causes  an  increase  in  cross  sectional  area  of  the  aircraft.  Minimizing  this 
increase  of  frontal  area  results  in  a  much  increased  weight  of  the  swivel  nozzle, 
hence  a  poorer  thrust  to  weight  ratio. 

Ventral  Nozzle 

The  ventral  nozzle  is  located  in  the  afterburner  section  of  the  engine.  A 
flow  blocker,  possibly  a  clamshell,  is  used  to  make  the  transition  from  the  horizontal 
to  vertical  mode.  The  ventral  nozzle  device  shown  in  Figure  7  has  a  clamshell  to 
close  off  part  of  the  horizontal  flow  and  open  up  part  of  the  ventral  nozzle  so  a  constant 
total  amount  of  flow  comes  from  the  engine  at  any  clamshell  position.  The  clamshell 
is  closed  until  all  of  the  flow  is  coming  out  the  ventral  nozzle  to  obtain  vertical 
thrust. 

During  horizontal  flight,  the  gas  path  is  unhampered  by  the  deflection  device. 
This  allows  the  maximum  augmentation  for  horizontal  operation  and  any  type  of  super¬ 
sonic  nozzle  can  be  used  with  the  ventral  nozzle  deflector. 

The  thrust  to  weight  of  the  ventral  nozzle  would  be  in  the  range  of  18  to  25. 

This  is  somewhat  heavier  than  the  swivel  nozzle,  but  produces  no  increase  in  aircraft 
frontal  area  due  to  the  deflector. 

Ventral  Nozzle  with  Rotating  Cascade 

The  ventral  nozzle  with  cascades  is  very  similar  to  that  without  cascades. 

The  flow  blocker,  however,  in  the  nozzle  with  cascades  completely  blocks  the  horizontal 
flow  at  the  same  time  uncovering  the  ventral  opening.  The  transition  from  the 


5 


horizontal  to  vertical  mode  is  then  made  by  the  rotation  of  the  cascade.  This  scheme 
is  shown  in  Figure  8. 

The  thrust  to  weight  of  this  deflector  system  would  also  be  in  the  range  of 

18  to  25. 

Aft-Hood  Deflector 

The  aft-hood  deflector  deflects  the  flow  by  a  curved  blocker  located  in  a 
supersonic  nozzle  downstream  of  the  afterburner.  Augmentation  is  available  for 
deflected  thrust  operation.  The  aft-hood  deflector  is  shown  in  Figure  9. 

The  transition  from  the  horizontal  to  vertical  mode  is  made  by  opening  doors 
in  the  bottom  of  the  ejector  shroud  and  rotating  the  hood  to  the  desired  position  for 
the  vertical  mode.  The  primary  or  afterburner  nozzle  is  then  opened  to  the  max- 
afterburning  position,  thus  transferring  the  minimum  jet  area  or  control  area  to  the 
deflector  exit.  This  is  done  so  that  the  flow  can  be  turned  at  a  low  Mach  number 
reducing  the  turning  losses  through  the  system. 

The  transition  from  the  vertical  to  horizontal  mode  is  made  by  rotating  the 
hood  to  the  stowed  position.  When  a  thrust  vector  of  approximately  45  degrees  is 
obtained,  the  control  area  is  switched  back  to  the  primary  nozzle  by  closing  down  the 
primary  nozzle.  The  hood  is  quickly  rotated  into  the  stowed  position,  and  the  ejector 
shroud  is  closed. 

The  gas  path  is  unhampered  by  the  deflection  system  and  max-augmentation 
is  possible  when  the  hood  is  stowed.  This  system  has  very  good  thrust  vectoring 
capability  and  supersonic  performance. 

The  thrust  to  w-eight  of  this  system  is  in  the  range  of  29  to  30,  or  somewhere 
between  the  ventral  nozzle  and  the  swivel  nozzle. 

EXPERIMENTAL  TEST  FACILITY 

Thrust  deflector  nozzle  performance  data  presented  in  this  paper  was  obtained 
from  scale  model  tests  in  the  United  Aircraft  Research  Laboratory  static  cold  flow 
facility.  The  facility  was  originally  designed  to  conduct  static  tests  on  airbreathing 
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engine  exhaust  models.  It  has  been  adapted  for  thrust  deflector  testing  by  the  addition 
of  a  rotating  plenum.  The  balance  and  rotating  plenum  are  shown  in  Figure  1C. 

The  rotating  plenum  is  used  to  rotate  the  force  vector  produced  by  the  thrust 
deflector  to  a  horizontal  position.  The  angle  of  the  rotating  plenum  or  the  angle  cf 
the  thrust  vector  is  recorded  while  the  force  magnitude  is  measured  by  the  balance. 
Forces  are  measured  with  a  thoroughly  developed  null  type  balance.  A  hollow  tube, 
supported  by  a  pair  of  flexures  to  the  balance  housing,  is  constrained  to  move  only 
in  the  thrust  or  drag  direction.  Balance  displacements  are  indicated  by  a  Schaevitz 
coil.  Nulling  is  achieved  by  adjustment  of  the  balance  base  pressure.  Once  *he 
balance  is  nulled,  the  pressures  are  recorded;  these  measurements  are  transformed 
into  nozzle  thrust  coefficients.  Airflow  measurements  are  made  in  accordance  with 
standard  ASME  techniques  using  a  bellmouth  flow  meter. 

TYPICAL  DEFLECTED  JET  PERFORMANCE 

Swivel  Nozzle 

A  swivel  nozzle  scale  model  of  a  nonaugmented  engine  system  is  shown  in 
Figure  11.  The  Mach  number  at  the  inlet  to  the  model  is  0. 35  and  the  model  is 
approximately  one  tenth  full  size.  The  performance  of  this  nozzle  at  various  nozzle 
pressure  ratios  and  rotation  angles  is  shown  in  Figure  12.  The  horizontal  performance 
is  2  to  3  percent  lower  than  a  convergent  nozzle  would  obtain  due  to  a  bend  in  the  duct 
upstream  of  the  swivel  bearing.  The  performance  falls  smoothly  during  transition 
by  another  2  to  3  percent  as  the  nozzle  is  rotated  due  to  worsening  of  the  internal 
geometry . 

The  flow  coefficient  is  nearly  constant  during  the  transition.  The  thrust  angle 
lags  the  nozzle  geometric  angle  for  60  and  90  degrees,  but  this  can  be  compensated  for 
by  overturning  the  nozzle  to  obtain  the  desired  thrust  angle. 

Ventral  Nozzle 

A  ventral  nozzle  is  shown  in  Figure  13.  The  model  simulates  a  possible  full 
scale  installation  having  a  ventral  nozzle  located  downstream  of  the  turbine  exit  plane 
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in  the  afterburner.  Afterburning  while  in  the  deflected  mode  with  the  ventral  nozzle 
in  this  position  would  be  difficult;  therefore,  the  model  jet  area  simulated  a  nonafter¬ 
burning  full  scale  engine  condition  for  the  vertical  mode.  The  inlet  Mach  number  to 
the  ventral  nozzle  tested  was  0. 20. 

The  performance  of  the  ventral  nozzle  is  shown  in  Figure  14.  The  thrust 
coefficient  is  4  to  6  percent  lower  than  a  convergent  nozzle.  The  results  were  obtained 
for  a  simple  design  and  it  is  possible  to  improve  the  performance  with  a  more  refined 
design. 

The  flow  coefficient  at  a  pressure  ratio  of  2. 0  is  about  10  percent  below  a 
conventional  convergent  nozzle.  The  hole  can  be  sized  on  this  basis  to  pass  the  eagine 
flow,  or  the  design  can  be  modified  to  improve  the  flow  coefficient.  The  transition 
from  the  horizontal  mode  of  this  deflection  system  is  made  by  closing  down  the 
horizontal  nozzle  and  opening  the  ventral  nozzle.  This  type  of  transition  is  shown  in 
Figure  15.  The  transition  performance  is  compromised  with  this  type  of  operation 
because  the  flow  is  divided  between  the  horizontal  and  ventral  nozzle. 

Ventral  Nozzle  with  Rotating  Cascade 

A  ventral  rotating  cascade  nozzle  model  is  shown  in  Figure  16.  This  device 
is  similar  to  the  ventral  nozzle  previously  described.  Performance  for  the  transition 
from  the  horizontal  to  vertical  mode  is  shown  for  the  case  of  blocking  completely  the 
horizontal  passage  and  opening  the  ventral  opening  while  the  cascade  is  rotated  for 
thrust  vectoring  the  model  simulated  nonafterburning  engine  operation.  The  inlet 
Mach  number  to  the  ventral  nozzle  was  0.  20.  This  model  is  identical  to  the  ventral 
nozzle  shown  above  except  for  the  rotating  cascade. 

The  performance  of  the  ventral  nozzle  with  cascade  is  shown  in  Figure  17  for 
90  degrees  thrust  deflection.  The  addition  of  the  cascade  to  the  ventral  nozzle  decreased 
the  thrust  by  another  4  percent  for  a  total  loss  of  almost  10  percenc  from  a  convergent 
nozzle.  Further  work  on  the  design  of  the  rotating  cascade  and  the  approach  to  the 
cascade  could  improve  the  performance;  however,  we  have  found  that  cascades  do  have 
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an  inherent  loss  associated  with  their  design. 

The  transition  performance  is  shown  in  Figure  18.  The  performance 
decreases  rapidly  as  the  cascade  is  rotated  due  to  the  vector  angle  of  the  flow  from 
the  cascade  which  is  not  measured.  This  side  force  would  be  canceled  with  a  two 
engine  installation.  Two  nozzles  would  be  necessary  to  cancel  the  side  force  in  a 
single  engine  installation.  It  is  obvious  that  this  deflector  is  not  desirable  below 
75  degrees  or  80  degrees  deflection. 

Aft-Hood  Deflector 

An  aft-hood  deflector  model  is  shown  in  Figure  19.  Since  this  deflection 
device  is  capable  of  partial  augmentation,  a  partial  augmentation  engine  condition 
was  simulated  by  the  model.  The  inlet  Mach  number  to  the  nozzle  was  0.315.  The 
performance  is  shown  in  Figure  20.  The  thrust  coefficient  of  this  system  is  good, 
within  2  or  3  percent  of  a  convergent  nozzle  at  full  deflection.  The  flow  coefficient 
is  lower  than  what  is  expected  from  a  convergent  nozzle;  however,  the  effective 
area,  AjQ,  is  sufficient  for  the  engine  condition.  If  an  unaugmented  engine  condition 
were  simulated  the  turning  losses  would  further  decrease.  The  transition  performance 
is  shown  in  Figure  21  for  a  pressure  ratio  of  two.  It  can  be  seen  from  the  figure  that 
the  aft-hood  transition  performance  is  also  good,  losing  only  2  to  3  percent  additional 
thrust  at  the  worst  condition. 

GROUND  EFFECTS 

The  use  of  deflected  jet  exhaust  for  vertical  thrust  raises  the  problem  of 
ground  proximity  effects  on  the  engine  exhaust  nozzle  performance.  The  deflected 
nozzle,  when  in  the  vertical  mode,  would  be  close  to  the  ground  causing  large 
changes  in  nozzle  performance.  A  series  of  tests  have  been  performed  to  determine 
the  effect  of  ground  proximity  on  nozzle  performance. 

The  performance  characteristics  presented  are  for  an  uninstalled  circular 
convergent  nozzle.  There  are  no  surfaces  except  the  nozzle  itself  and  the  ground 
upon  which  the  flow  acts.  The  performance  is  shown  in  Figure  22.  It  is  seen  that  the 


nozzle  performance  decreases  rapidly  with  decreasing  distance  to  the  ground.  This 
means  that  serious  engine  flow  suppression  as  well  as  a  loss  in  thrust  will  occur  if 
the  nozzle  exit  is  positioned  too  close  to  the  ground. 

TRANSITION  TIME 

A  simple  study  was  made  to  illustrate  that  the  transition  time  from  hovering 
to  axial  flight  of  a  given  VTOL  aircraft  can  be  affected  by  deflector  characteristics. 
Two  lift  cruise  engine  deflection  systems  are  compared,  the  ventral  nozzle  and  the 
ventral  nozzle  with  rotating  cascades.  These  two  systems  were  used  because  they 
would  require  no  changes  in  the  lift  drag  relationships  of  the  aircraft. 

A  typical  lift  cruise  airplane  having  60, 000  pounds  gross  weight  was  assumed 
for  this  analysis,  powered  by  two  lift  cruise  engines  anf  our  lift  engines.  The  lift 
cruise  engines  have  a  take-off  thrust  of  24,000  pounds  each  and  the  lift  engine  6,000 
pounds  each.  The  total  available  vertical  thrust  is  72,000  poundf. 

The  transition  time  is  shown  in  Figure  23.  The  results  show  that  the  greater 
horizontal  thrust  of  the  cascade  initially  increases  the  velocity  of  the  airplane  at  a 
faster  rate  than  the  ventral  nozzle  without  the  cascade.  After  about  ten  seconds,  the 
opening  of  the  axial  nozzle  results  in  a  faster  rate  of  acceleration  for  the  ventral 
nozzle  without  the  cascade.  The  ventral  nozzle  with  the  rotating  cascade  completes 
the  transition  more  slowly  than  the  ventral  nozzle  without  the  cascade  but  is  at  a 
higher  altitude. 

This  comparison  reaches  no  conclusions  other  than  the  fact  that  the  choice 
of  deflector  has  a  significant  effect  on  transition  time.  Obviously,  a  short  transition 
time  is  desirable,  but  the  time  must  be  evaluated  against  the  flight  path  (altitude  vs 
distance)  and  stability. 

CONCLUSIONS 

It  was  found  that  the  inner  wall  of  a  bend  shows  the  largest  contribution  to  the 
total  pressure  loss  of  the  bend.  Good  thrust  performance  will  be  obtained  when  the 
inlet  Mach  number  to  the  bend  is  low  and  the  turning  radius  ratio  is  approximately  2.  5 
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duct  diameters. 

The  conclusions  reached  in  the  paper  for  each  of  the  deflection  devices  are 
shown  in  the  following  comparison. 


90°  Deflec¬ 
tion  Thrust 
Coef.  Cvr 

Vectoring 

Capability 

Supersonic 

Capability 

Effect  on 

Horizontal 

Performance 

Deflector 

Thrust/ 

Weight* 

Swivel 

Nozzle 

0.94 

Good 

Limited 

1  to  3% 
Decrease 

35  to  50 

Ventral 

Nozzle 

0.93 

Limited 

Good 

No  Effect 

18  to  25 

Ventral 
with  Rotat¬ 
ing  Cascade 

0.89 

Limited 

Good 

No  Effect 

18  to  25 

Aft- Hood 

0.96 

Good 

Good 

No  Effect 

20  to  30 

Deflector 

♦Vertical  Thrust  Divided  by  Weight  added  by  Deflection  Device 

The  performance  levels  shown  are  for  rather  crude  test  models.  In  each 
case  one  can  expect  to  refine  these  models  and  airflow  characteristics  and  reach  higher 
performance  levels.  Improvements  in  thrust  to  weight  ratios  are  also  possible. 

On  the  other  hand,  approximately  the  same  degree  of  effort  went  into  each  of 
these  designs  and  it  is  possible  that  the  deflector  that  best  satisfies  the  requirements 
of  a  good  thrust  deflector  is  the  Aft-Hood  Deflector. 

REFERENCES 

1)  Higginbotham,  Wood,  Valentine:  "A  Study  of  the  High-Speed  Performance 
Characteristics  of  90  Degree  Bends  in  Circular  Ducts",  NACA  TN  3696, 

June  1956. 

2)  Wilbur,  Stafford  W. :  "An  Investigation  of  Flow  in  Circular  and  Annular 

90  Degree  Bends  with  a  Transition  of  Cross  Section",  NACA  TN  3995,  June 
1957. 


11 


3)  Kaufmann,  Walter:  "Fluid  Mechanics",  McGraw  Hill  1963. 

4)  Horgan,  Migdal:  "Thrust  Nozzles  for  Supersonic  Transport  Aircraft",  ASME 
63-AHGT-73,  January  1963. 


NOMENCLATURE 


PT 

P 

D 

r/D 

M 

Fgi 

Q 

A 

6 

Cvr 

q 


Total  pressure 
Static  pressure 
Diameter 

Turning  radius  ratio 

Mach  number 

Ideal  gross  thrust  =  AjQPt 


Flow  coefficient  equal  to  the  measured-to-ideal  weight  flow 
Physical  area 

Angle  between  the  engine  centerline  and  the  deflected  gross  thrust 

vedtor  in  the  vertical  plane 

Deflected  gross  thrust/ideal  gross  thrust 

Dynamic  pressure 


SUBSCRIPTS 

1 

2 

J 

a 


Region  before  flow  turn 
Region  after  flow  turn 
Nozzle  exit 
ambient 


12 


LOCAL  PRESSURE  LOSS  COEFFICIENT 


-OUTER  WALL 


Figure  3  Total  Pressure  Loss  Distribution 


15 


0.10*- 

0 


■  _ ■  1 _  ■  I 

1.0  2.0  3.0  4.0  5.0 

TURNING  RADIUS  RATIO -r/D 


otal  Pressure  Loss  Variation  with  Turning  Radius 


16 


Figure  0  Swivel  Nozzle  Operation  Axial  Thixist  Position 
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Figure  7  Ventral  Nozzle  Installation  VTOL  Position 
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Figure  15  Ventral  Nozzle  Transition  Performance 
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Figure  19  Aft- Hood  Deflector  iVtodel 
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Figure  21  Aft-Hood  Deflector  Transition  Performance 
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Figure  22  Ground  Proximity  Effects  on  Nozzle  Performance 
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Figure  23  Transition  Time 
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SHROUDED  PROPELLER  RESEARCH  AT 
MISSISSIPPI  STATE  UNIVERSITY  LEADING  TO 
APPLICATION  ON  THE  UNITED  STATES  ARMY  XV-11A 

by 

JOE  F.  THOMPSON  and  S.  C.  ROBERTS 

Aerophysics  Department,  Mississippi  State  University 
State  College,  Mississippi,  U.S.A. 


INTRODUCTION 

The  failure  of  the  conventional  propeller  to  provide  power  for  climb 
at  speeds  below  36  m.p.h.  on  an  Army  L-21  Cub  with  suction  boundary  layer 
control  provided  the  impulse  for  shrouded  propeller  research  at  Mississippi 
State  University  some  ten  years  ago.  As  is  well  known,  the  static  perfor¬ 
mance  of  an  open  propeller  designed  for  maximum  efficiency  at  cruising  speed 
leaves  much  to  be  desired.  However,  high  static  and  low  speed  thrust  is  just 
as  much  a  requirement  as  rs  high  lift  for  the  achievement  of  STOL  operation. 
Although  shrouded  propellers  have  been  used  for  some  time  to  increase  the 
static  thrust  of  tug  boats,  only  recently  has  the  interest  in  STOL  operation 
prompted  their  employment  on  aircraft. 

The  circulation  about  the  airfoil  section  of  the  shroud  induces  a 
velocity  increment  at  the  propeller  plane,  especially  near  the  propeller 
tips.  This  induced  velocity  compensates,  to  some  degree,  for  the  wide  change 
of  angle  of  attack  near  the  tips  with  flight  speed,  and  thereby  allows  the 
propeller  to  operate  at  near  optimum  conditions  over  a  wider  range  of  flight 
speeds.  With  a  properly  designed  shroud  a  higher  efficiency  can  be  obtained 
for  the  same  thrust,  or  the  thrust  can  be  increased  without  a  loss  in  effi¬ 


ciency. 


This  circulation  about  the  shroud  also  produces  a  "lift"  vector  on 
the  shroud,  which  when  angled  forward  results  in  a  thrust  component  on  the 
shroud  which  can  be  in  excess  of  the  propeller  thrust.  Finally,  when  the 
clearance  between  the  propeller  and  the  shroud  is  made  small,  the  shroud 
effectively  end-plates  the  propeller,  reducing  the  losses  due  to  vortex 
shedding. 

The  combination  of  these  effects  results  in  a  considerable  increase 
in  thrust,  especially  in  the  static  condition  and  in  low  speed  flight.  On 
the  minus  side,  the  shroud  is  an  additional  drag  producing  element.  However, 
with  careful  design  the  break-even  velocity  (flight  velocity  at  which  shroud 
drag  equals  shroud  thrust)  can  be  made  quite  large. 
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SHROUD  LEADING  EDGE  SEPARATION 


Early  Model  Studies 

Early  studies  on  model  shrouded  propellers  at  Mississippi  State 
University  revealed  the  necessity  of  preventing  leading  edge  separation  on 
the  shroud  if  its  full  potential  for  thrust  augmentation  is  to  be  realized. 
This  separation  occurs  in  precisely  the  case  when  the  increased  thrust  from 
the  shroud  is  most  important  -  the  static  and  low  speed  conditions.  At  zero 
or  low  flight  speed  the  flow  must  turn  around  the  shroud  leading  edge  from 
behind,  and  is  apt  to  separate  just  downstream  of  the  leading  edge  on  the 
inside  of  the  shroud  (Figure  1).  The  possibility  of  this  separation  is 
heightened  by  the  fact  that  a  constant  pitch  propeller  operates  inefficiently 
at  the  tips  due  to  the  angle  of  attack  change  with  radius.  Once  separation 
occurs,  the  portion  of  the  propeller  within  the  separated  region  is  ineffec¬ 
tive  or  even  detrimental . 

Full-Scale  Studies 

Full-scale  tests  were  conducted  on  a  shroud  mounted  around  the 
pusher  propeller  of  an  Anderson-Greenwood  AG-14  aircraft  (Figure  2).  Leading 
edge  separation  occurred  on  the  shroud  with  an  untwisted  propeller  (Figure  3). 
However,  when  the  inflow  velocity  profile  across  the  propeller  plane  was 
measured  and  the  propeller  was  twisted  to  match  this  angle  of  attack  distri¬ 
bution,  separation  was  suppressed  (Figure  3),  and  thrust  increased  signifi¬ 
cantly.  With  this  configuration,  shroud  thrusts  even  greater  than  propeller 
thrusts  were  achieved. 

Separation  Suppression  by  Boundary  Layer  Control 
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The  possibility  of  suppressing  shroud  leading  edge  separation  by 
suction  boundary  layer  control  in  the  shroud  was  investigated,  again  on 
the  shroud  on  the  AG-14.  The  severity  and  extent  of  the  separation  was 
considerably  reduced  from  the  case  of  no  boundary  layer  control.  In  this 
particular  case  the  application  of  boundary  layer  control  resulted  in  an 
increase  in  thrust  from  350  pounds  to  490  pounds  with  an  expenditure  of  5 
horsepower  of  suction  power.  As  a  point  of  interest,  the  original  thrust 
of  this  aircraft  with  the  open  propeller  was  280  pounds.  This  investigation 
also  showed  that  headwinds  can  cause  separated  flow  to  attach,  and  tailwinds 
can  lead  to  separation.  The  flow  is  also  very  sensitive  to  crosswinds,  with 
separation  occurring  on  the  windward  side  and  attachment  on  the  leeward. 

THE  SHROUDED  PROPELLER  OF  THE  XAZ-1  MARVELETTE4 

Design 

The  shrouded  propeller  on  the  AG-14  was  subjected  to  extensive 
testing,  both  in  the  static  condition  and  in  flight  (Figure  4).  This  full- 
scale  testing  provided  the  experience  and  design  principles  necessary  for 
the  design  of  the  shrouded  propeller  to  be  used  on  the  XAZ-1  MARVELETTE 
research  aircraft.  The  XAZ-1  MARVELETTE  (Figure  5)  was  built  as  a  test  bed 
for  the  high  lift,  high  thrust,  low  drag  concepts  to  be  employed  in  the 
design  of  the  U.  S.  Army  XV-11A  STOL  research  aircraft. 

There  are  several  reasons  that  dictate  the  use  of  a  shrouded  propel¬ 
ler  in  a  pusher  configuration,  rathe:  then  a  '•.ractor  configuration.  The 
shroud  is  a  major  stability  component  and  must  be  located  aft  of  the  center 
of  gravity  of  Che  aircraft,  or  else  a  tremendous  fin  area  will  be  required 
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to  overcome  the  adverse  yaw  effect  of  the  shroud.  Also,  it  had  been  found 
in  research  on  a  nigh-lift  L-19  that  the  annular  region  of  low  dynamic  pres¬ 
sure  surrounding  an  open  propeller  caused  early  separation  upon  impingement 
on  the  wings,  thus  negating  the  boundary  layer  control  system  on  the  wings 
and  limiting  the  maximum  lift  attainable.  Therefore,  the  shroud  should  be 
located  behind  the  wings  on  an  STOL  aircraft  in  order  to  assure  smooth, 
undisturbed  flow  over  the  wings.  A  third  reason  for  an  aft  location  of  the 
shroud  is,  of  course,  visibility. 

The  low-drag  airfoil  section  used  for  the  shroud  was  mounted  at 
zero  angle  of  attack.  The  leading  edge  radius  was  increased  to  suppress 
possible  leading  edge  separation.  The  shroud  diffusion  angle  was  approxi¬ 
mately  5.3  degrees,  which  is  Pmall  enough  to  prevent  possible  diffuser  sepa¬ 
ration  problems.  The  three-bladed  MARVELETTE  propeller  was  of  variable  pitch, 
with  blades  of  constant  chord,  and  in  the  original  configuration  the  blades 
were  untwisted.  This  allowed  twist  modification  to  be  easily  incorporated 
from  inflow  measurements  from  ground  running  and  flight  tests  for  optimum 
propeller  performance.  The  propeller  was  driven  by  a  fiber  glass  shaft 
attached  to  a  Continental  C-90  engine  which  developed  90  horsepower  at  2,500 
r . p.m. 

The  shroud  was  made  of  fiber  glass.  The  longitudinal  and  directioral 
control  surfaces  were  incorporated  in  the  shroud  i.t  order  to  eliminate  the 
conventional  tail  sections  and  reduce  the  drag  at  higher  speeds  (Figure  6). 
These  controls  were  constructed  so  that  very  little  gap  existed  between  the 
control  surfaces  at  the  trailing  edge  of  the  shroud,  irrespective  of  the  con¬ 
trol  surface  position.  Two  segments  on  either  side  of  the  shroud  acted  as 
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rudders,  and  four  segments,  top  and  bottom,  served  as  elevators.  Massive 
fiber  glass  struts  maintained  the  close  alignment  of  the  shroud  and  carried 
the  shroud  thrust  loads. 

Ground  Test  Results 

Since  it  was  necessary  to  determine  whether  any  serious  separation 
or  irregularities  occurred  in  the  flow  around  the  shrouded  propeller,  a 
considerable  number  of  flow  visualization  studies  were  performed  using  tuft 
and  napthalene  sublimation  techniques.  Intermittent  leading  edge  separation 
on  the  shroud  is  evident  in  the  tuft  pictures  in  Figure  7.  However,  the 
separation  region  is  small  and  the  flow  is  generally  quite  smooth.  Sublima¬ 
tion  studies  revealed  considerable  separation  from  the  spinner  and  the  blade 
root  cuffs.  This  was  reduced  by  modification  of  the  cuffs.  Sublimation  also 
showed  the  flow  over  the  propeller  blades  to  be  laminar  to  a  large  extent, 
and  the  tip  clearance  to  be  sufficiently  small  to  reduce  or  to  nearly  elimi¬ 
nate  tip  vortex  effects. 

The  inflow  velocity  was  examined  at  several  angular  positions  and  was 
quite  symmetrical,  with  only  slight  deviations  occurring  at  the  bottom  near 
the  fuselage  (Figure  8).  There  was  considerable  variation  in  inflow  velocity 
with  radius.  Changes  in  camber  of  the  variable  camber  wing,  or  equivalent 
flap  angle,  up  to  30  degrees  had  very  little  effect  on  the  flow  into  the 
shrouded  propeller. 

Changes  in  aircraft  airspeed  had  a  very  detrimental  effect  on  propel¬ 
ler  blade  effectiveness,  as  can  be  clearly  seen  in  Figure  9,  where  the  propel¬ 
ler  blade  was  in  the  original  untwisted  condition.  At  an  airspeed  of  70  m.p.h., 
the  outer  two-thirds  of  the  blade  was  operating  with  only  a  small  radial 
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portion  at  an  optimum  angle  ct  attack.  At  110  m.p.h.,  only  the  outer  one- 
third  was  operating,  and  the  remainder  of  the  propeller  was  windmilling. 

This  was  probably  the  reason  for  the  nonuniform  inflow  velocity  distributions, 
and  was  undoubtably  the  reason  for  the  poor  climb-out  performance  of  the  air¬ 
craft  experienced  at  90  m.p.h.  Twisting  the  propeller  blade  to  match  the  in¬ 
flow  velocity  distribution  helped  to  relieve  this  problem,  as  can  be  seen  from 
the  results  using  the  twisted  blade  in  Figure  10.  At  110  m.p.h.,  ovr,  half 
of  the  blade  was  operating  at  a  reasonably  constant  angle  of  attack.  This 
modification  considerably  improved  the  climb-out  characteristics  of  the 
MARVELETTF  aircraft. 

The  static  thrust  generated  by  the  shrouded  propeller  of  the  MARVELETTE 
was  found  from  tethered  tests.  Although  the  tests  were  performed  early  in  the 
morning  when  the  wind  was  generally  calm,  the  aircraft  was  located  inside  the 
hangar  to  eliminate  any  crosswinds  which  could  possibly  promote  flow  separa¬ 
tion  on  the  shroud.  Readings  were  also  taken  with  the  moveable  control  sur¬ 
faces  in  various  positions,  but  it  was  found  that  control  positions  had  little 
or  no  effect  on  static  thrust.  The  static  thrust  of  the  shrouded  propeller 
was  approximately  100  percent  greater  than  the  static  thrust  of  the  equivalent 
open  propeller  driven  by  the  same  horsepower  engine. 

Some  dynamic  thrust  measurements  were  taken  by  towing  a  ground  vehi¬ 
cle,  but  this  method  was  dangerous  and  impractical  at  speeds  greater  than  30 
m.p.h.  In  order  to  obtain  thrust  measurements  at  higher  speeds,  an  instru¬ 
ment  was  developed  which  would  give  an  accurate  and  repeatable  plot  of  air¬ 
craft  displacement  from  a  fixed  point  as  a  function  of  time.  This  technique 
of  measuring  the  variations  of  thrust  with  forward  velocity  was  very  effective 


up  to  speeds  below  the  tzkeoff  airspeed  of  the  aircraft.  With  careful  pilot¬ 
ing  technique,  it  could  be  used  at  ground  speeds  in  excess  of  the  minimum 
flying  speed  of  the  aircraft. 

From  Figure  11,  the  theoretical  break-even  velocity  was  152  f.p.s., 
assuming  that  the  ring  tall  is  additional  to  the  conventional  tailplane; 
however,  in  the  MARVELETTE  case,  where  the  conventional  tailplane  was  elimi¬ 
nated,  the  break-even  velocity  was  162  f.p.s.  The  experimental  curve  of  total 
thrust  is  considerably  below  the  theoretical  curve  from  simple  momentum 
theory  (Figure  12).  This  is  to  be  expected  because  in  the  theoretical  case 
it  was  assumed  that  the  complete  propeller  was  operating  at  an  optimum  blade 
angle  of  attack  which,  of  course,  was  not  the  case.  Both  curves  have  essen¬ 
tially  the  same  shape,  and  the  experimental  curve  is  approximately  75  percent 
of  the  theoretical  curve.  The  differences  between  theoretical  and  experimental 
results  are  due  primarily  to  interference  effects  and  to  insufficient  or 
improper  blade  twist,  especially  around  the  blade  root.  Shroud  leading  edge 
laminar  separation  bubbles  and  flow  separation  around  the  spinner  also  contrib¬ 
ute  to  the  difference  in  results. 

The  simple  momentum  theory^  involves  a  number  of  simplifying  assump¬ 
tions  such  as  no  rotation  of  flow  leaving  the  shroud,  constant  inflow  velocity, 
wake  area  equal  to  shroud  exit  area,  and  static  pressure  at  shroud  exit  equal 

to  ambient  pressure.  The  method  of  singularities,  as  has  recently  been  very 

12  3 

elegantly  developed  by  Therm  Advanced  Research,  Inc.  *  *  ,  is  far  superior, 
provided  that  its  requirements  are  met.  However,  due  to  the  fact  that  shrouded 
propellers  usually  operate  on  the  rear  of  the  vehicles  where  the  inflow  veloc¬ 
ity  consists  of  a  combination  of  free-stream  velocity  and  body  wake  velocity, 
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it  is  difficult  to  design  an  optimum  shrouded  propeller  without  knowing  the 
influence  of  the  vehicle  on  the  inflow  velocity  field  of  the  propeller. 

Flight  Test  Results 

As  a  result  of  flight  testing  of  the  MARVELETTE  (Figure  13),  it  was 
necessary  to  Install  additional  jury  strut."  on  the  shroud  to  maintain  the 
close  propeller  tip  clearance.  Early  in  these  flights  it  became  clear  that 
the  large  moment  of  inertia  about  the  lateral  and  vertical  axes  resulted  in 
a  tendency  toward  an  oscillation,  particularly  in  pitch.  The  aircraft  was 
completely  manageable,  but  increased  damping  in  pitch  appeared  desirable. 
Consequently,  small  fins  were  added  to  the  four  spaces  between  the  rudders 
ad  elevators.  These  fins  were  expected  to  Increase  the  tall  volume,  make 
the  outside  of  the  duct  more  effective,  and  end-plate  the  movable  surfaces  to 
improve  their  effectiveness. 

The  added  fins  apparently  succeeded  in  improving  both  steadiness  and 
controllability,  although  some  pitching  oscillation  was  still  present.  It  is 
probable  that  most  of  the  oscillation  resulted  from  lost  motion  in  the  eleva¬ 
tor  control  linkage.  The  pitch  oscillation  and  a  slight  tendency  to  wallow 
was  most  pronounced  with  reduced  power  settings  used  for  landing  approach. 
This  unsteadiness  was  probably  the  result  of  reduced  flow  through  the  duct 
at  low  power  settings.  Stability  was  generally  good,  however,  and  there  was 
no  difficulty  in  making  coordinated  turns.  So  much  directional  stability  was 
present  that  little  use  of  the  rudder  was  required.  The  need  for  pitch  con¬ 
trol  appears  to  be  considerably  greater  with  a  ducted  propeller  than  with  an 
open  one. 

The  most  important  conclusion  which  can  be  drawn  from  the  results 
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of  the  flight  testing  program  is  that,  although  the  MARVELETTE  was  an  uncon¬ 
ventional  aircraft,  the  general  handling  characteristics  were  very  similar 
to  those  of  a  conventional  airplane.  Tha  ducted  propeller  performed  satis¬ 
factorily,  with  sufficient  thrust  available  to  accelerate  the  airplane  from 
0  to  70  m.p.h.  in  17  seconds,  even  though  the  blade  twist  angles  were  not 
yet  optimized  for  maximum  thrust  at  low  forward  velocities. 

The  major  effect  of  the  shroud  was  thrust  augmentation, although  the 
directional  stability  contribution  was  sufficient  to  overshadow  the  rudder 
effect  in  normal  flying.  The  flight  results  showed  that  the  assumptions  used 
in  designing  the  shroud  and  control  surfaces  were  adequate  to  insure  adequate 
control  of  the  airplane. 

THE  SHROUDED  PROPELLER  OF  THE  XV-11A 

Design 

The  results  of  static  and  flight  tests  on  the  shrouded  propeller  of 
the  MARVELETTE  directed  the  design  of  the  U.  S.  Army  XV-11A  (Figure  14),  an 
aerodynamic  research  aircraft  specifically  designed  to  explore  the  problem 
areas  inherent  in  wide  speed  range,  fixed-wing  airplanes  of  interest  in  Army 
operations.  These  problem  areas  include  high  lift,  low  drag,  high  thrust, 
and  the  compatibility  of  these  characteristics  with  the  requirements  for  good 
visibility,  large  usable  volume,  and  rough  field  or  pantabase  operation. 

The  shroud  is  fiber  glass  and  of  the  same  airfoil  section  (NACA  4416) 
used  on  the  MARVELETTE,  with  the  leading  edge  radius  again  increased  to  sup¬ 
press  separation.  The  shroud  inside  diameter  is  66^  inches,  and  the  chord 
is  30  inches.  The  propeller  is  of  7  inch  chord  and  is  a  standard  Aero  Products 
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propeller  which  was  designed  for  the  T-63  engine,  except  that  it  had  the 
outer  eleven  (11)  inches  removed  to  fit  into  the  shroud.  The  clearance 
between  propeller  and  shroud  is  3/32  inch  when  the  propeller  is  still.  The 
clearance  decreases  somewhat  when  the  propeller  is  operating,  because  of 
blade  expansion  under  centrifugal  force.  The  propeller  is  powered  by  an 
Allison  T-63  gas  turbine  engine,  developing  250  horsepower  at  2,000  r.p.m. 
Conventional  tail  control  surfaces  are  used  (Figure  15)  rather  than  those 
incorporated  in  the  shroud  as  on  the  MARVELETTE.  The  previously  mentioned 
pitch  control  problems  encountered  with  the  MARVELETTE,  plus  the  fact  that 
the  heavy  control  linkage  system  inside  the  MARVELETTE  shroud  caused  some 
deformation  resulting  in  some  hits  by  the  propeller,  prompted  this  change. 

Tests  Results 

The  XV- 1 1A  began  its  flight  testing  (Figure  16)  in  December  1965, 
and  no  difficulties  have  yet  been  encountered  with  the  propulsion  system. 

The  shrouded  propeller  develops  1,400  pounds  of  thrust,  giving  a  thrust  per 
horsepower  of  5.6  lb/HP.  Since  the  aircraft  gross  weight  is  only  2,400 
pounds,  it  is  well  capable  of  STOL  operation.  In  fact,  the  thrust  at  the 
engine  flight  idling  position  is  so  great  that  the  propeller  must  be  put 
in  reverse  pitch  for  landing  in  short  distances.  Flight  speeds  of  150  m-D.h. 
have  been  attained  at  43  percent  power  levels,  and  the  break-even  speed, 
where  -oud  drag  exceeds  shroud  thrust,  is  estimated  to  be  160  m.p.h. 

SHROUD  FLOW  PATTERNS  IN  REARWARD  AND  FORWARD  MOTION 

Although  many  investigations  have  been  directed  toward  determining 
the  thrust  and  overall  performance  characteristics  of  shrouded  propellers, 
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only  little  consideration  has  been  given  to  determining  the  flow  pattern 
abcut  the  shroud. 

Shrou i  thrust  calculations  are  performed  by  the  method  of  distribu- 
1-3 

ted  singularities  which  assumes  a  stagnation  line  at  the  trailing  edge, 
satisfying  the  Kutta  condition.  In  reality  viscous  interaction  of  the 
propeller  slipstream  with  the  external  flow  alters  this  simple  flow  pattern 
so  that,  in  addition  to  the  stagnation  line  at  the  trailing  edge,  a  stagna¬ 
tion  line  exists  on  the  outer  surface  of  the  shroud.  For  instance,  in  the 
static  condition  the  inviscid  flow  pattern  shown  in  Figure  17  is  altered  by 
slipstream  interaction  so  that  the  actual  flow  pattern  is  that  shown  in 
Figure  18.  Inviscid  flow  theory  would  have  the  two  stagnation  lines  coinci¬ 
ding  at  the  trailing  edge.  This  forward  movement  of  one  stagnation  line 
reduces  the  circulation  below  that  which  would  be  calculated  by  the  distribu¬ 
ted  singularities  methods.  Thus  the  shrcud  thrust  will  be  less  than  the 
theoretical  value,  the  reduction  being  a  function  of  the  distance  between 
the  outer  surface  stagnation  line  and  the  trailing  edge. 

Therefore,  an  extensive  investigation  was  undertaken  to  determine 
the  flow  pattern  and  pressure  distribution  about  the  shroud  of  a  shrouded 
propeller  in  the  static  condition,  in  rearward  motion,  and  in  forward  motion. 
The  inclusion  of  the  case  of  rearward  motion  was  prompted  by  the  recent 
application  of  shrouded  propellers  on  VTOL  aircraft.  For  static  testing  the 
propeller  was  mounted  on  a  propeller  test  stand  as  shown  in  Figure  19.  The 
entire  propeller  test  stand  used  in  the  static  testing  was  mounted  on  an 
automobile  chassis,  facing  aft  as  shown  in  Figure  20,  for  rearward  mother, 
studies  and  facing  forward  (Figure  21)  for  forward  motion  studies. 
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Typical  stagnation  lines  are  shown  in  the  shroud  tuft  pictures  in 
Figure  22.  Since  fluctuations  in  the  flow  made  the  stagnation  lines  some¬ 
what  irregular  in  some  cases,  the  locations  used  in  data  correlation  were 
averages  over  all  tufts  on  the  four  pictures  taken  in  each  case.  Figure 
23  shows  the  stagnation  line  location  as  a  function  of  the  shroud  thrust 
coefficient.  At  5°  pitch  angle  the  stagnation  line  is  located  fairly 
far  forward,  around  the  three-quarter  chord  position.  As  the  shroud  circu¬ 
lation  increases,  either  because  of  pitch  angle  or  rotational  speed  increase, 
the  stagnation  line  moves  toward  the  trailing  edge,  rapidly  at  first  and 
then  more  slowly,  approaching  a  chord  fraction  between  0.90  and  0.92. 

Although  there  is  some  irregularity  due  to  the  flow  fluctuations,  the  curves 
are  of  the  same  general  shape  for  all  rotational  speedp. 

Rearward  Motion 

In  rearward  motion  essentially  three  types  of  shroud  pressure  dis¬ 
tribution  were  observed  as  sketched  in  Figure  K  .  Type  A,  occurring  at  the 
higher  circulations  and  lower  speeds,  is  the  same  type  of  pressure  distribu¬ 
tion  as  observed  in  the  static  case.  Type  B,  occurring  at  intermediate 
circulations  and  speeds,  is  characterized  by  the  occuri  2nce  of  positive 
values  of  the  pressure  coefficient  on  the  inner  shroud  surface  and  one 
intersection  in  the  curve.  Type  C,  at  lower  circulations  and  higher  speeds, 
exhibits  two  intersections  with  lower  pressure  on  the  outer  surface  over 
most  of  the  shroud.  The  transition  from  Type  A  to  Type  B  is  much  more 
sudden  than  that  from  B  to  C.  In  some  cases  two  types  of  pressure  distribu¬ 
tion  were  observed  at  the  same  pitch  angle,  rotational  speed,  and  linear 
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speed.  Thus  the  occurrence  of  a  particular  type  of  pressure  distribution 
seems  to  be  a  question  of  stability. 

Typical  stagnation  lines  and  conditions  of  the  wire  tufts  above, 
in  front  of,  and  behind  the  shroud  are  shown  in  Figure  25.  The  flow  patterns 
about  the  shroud  as  indicated  by  the  wire  tufts  are  drawn  in  Figure  26.  At 
all  speeds  the  shroud  is  surrounded  by  a  ring  vortex  which  encloses  another 

ring  vortex  on  the  outer  surface  of  the  shroud.  For  the  sake  of  brevity 

the  0.0  chord  position  will  be  referred  to  as  the  leading  edge  and  the  1.0 
chord  position  as  the  trailing  edge  in  both  rearward  and  forward  motion, 

even  though  in  the  rearward  case  these  terms  are  misnomers.  The  transition 

of  flow  patterns  proceeds  as  follows: 

(1)  (Figure  26a)  At  low  speed  the  enclosed  vortex  on  tne  outer 
surface  is  located  near  the  trailing  edge  and  is  very  small.  This  condition 
results  in  a  stagnation  line  located  at  high  chord  position.  The  flow 
turns  sharply  around  the  leading  edge  of  the  shroud  causing  a  low  pressure 
peak  near  the  leading  edge  and  resulting  in  a  pressure  distribution  of  Type 
A,  similar  to  that  of  the  static  condition. 

(2)  (Figure  26a)  As  speed  increases,  the  enclosed  vortex  grows 
and  the  stagnation  line  moves  toward  lower  chord  positions.  The  pressure 
distribution  remains  of  Type  A,  but  the  low  pressure  peak  decreases  since 
the  turn  about  the  leading  edge  is  more  gradual. 

(3)  (Figure  26b)  When  the  stagnation  line  nears  the  leading  edge, 
the  enclosed  vortexhas  become  quite  large.  The  turn  around  the  leading 
edge  is  very  gradual  so  that  the  low  pressure  peak  has  disappeared.  The 
increase  in  pressure  across  the  propeller  now  results  in  positive  pressure 
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coefficients  on  the  inner  surface  near  the  trailing  edge, since  tho  pressure 
upstream  of  the  propeller  has  increased.  This  results  in  a  pressure  distri¬ 
bution  of  Type  B. 

(A)  (Figure  26b)  Further  increase  in  speed  causes  the  enclosed 
vortex  to  be  inclined  rearward  with  a  small  stagnation  line  movement  toward 
larger  chord  positions.  The  increased  velocity  on  the  outer  surface  causes 
the  pressure  on  the  outer  surface  to  decrease.  However,  the  lowest  pressures 
still  .ccur  near  the  leading  edge  on  the  inner  surface  so  that  the  pressure 
distribution  is  still  of  Type  B. 

(5)  (Figure  26c)  Still  further  speed  increase  causes  the  inclina¬ 
tion  of  the  enclosed  vortex  to  increase  so  that  it  becomes  pressed  down  on 
the  outer  surface  of  the  shroud.  This  increased  velocity  on  the  outer 
surface  of  the  shroud  results  in  lower  pressures  on  the  outer  surface  than 
on  the  inner.  The  stagnation  line  reverses  directions  again  and  moves 
toward  the  leading  edge  as  the  vortex  is  pressed  down  on  the  surface. 

In  Figure  27  the  shroud  thrust  coefficient  is  plotted  against  the 
advance  ratio,  with  the  pressure  distribution  type  given  for  each  point. 

The  points  are  well  correlated  by  one  curve  for  each  pitch  angle.  Although 
some  several  points  fall  off  these  curves,  most  of  these  points  represent 
pressure  distributions  of  types  not  consistent  with  those  of  neighboring 
points  on  the  curves.  This  results  from  the  above  mentioned  temporary 
occurrence  of  pressure  distributions  of  a  type  not  the  most  stable  at 
the  particular  values  of  pitch  angle,  rotational  speed,  and  linear  speed. 

The  shape  of  the  lines  at  high  advance  ratio  is  questionable  because  of 
the  scarcity  of  points  in  that  region. 
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The  stagnation  line  location  is  plotted  in  Figure  28  as  a  function 
of  advance  ratio.  As  mentioned  above,  the  trend  is  for  the  stagnation  line 
to  move  toward  the  leading  edge  as  J  increases.  Upon  reaching  the  vicinity 
of  the  leading  edge,  the  motion  of  the  stagnation  line  is  reversed,  and  it 
moves  toward  higher  chord  positions  as  J  increases.  Still  further  increase 
in  J  causes  another  reversal  and  movement  toward  the  leading  edge  again. 

t igure  29  shows  a  correlation  of  shroud  thrust  coefficient  against 
stagnation  line  location.  Again  ic  is  possible  to  draw  one  line  for  each 
pitch  angle.  However,  in  this  case,  errors  in  measurement  compound  each 
other  to  produce  more  scatter. 

Forward  Motion 

Typical  stagnation  lines  and  conditions  of  the  wire  tufts  in  for¬ 
ward  motion  are  shown  in  Figure  30.  Flow  patterns  are  drawn  in  Figure  31. 
Even  at  low  speeds  the  stagnation  line  is  located  forward  of  the  quarter 
chord  point, yielding  a  flow  pattern  as  shown  in  Figure  31a. 

Large  increases  in  speed  cause  the  stagnation  line  to  move  to  the 
leading  edge,  with  a  flow  pattern  as  in  Figure  31b.  As  speed  is  further 
increased,  leading  edge  separation  and  reattachment  occurs  on  the  outer 
surface.  The  speed  at  which  separation  occurs  is  lower  at  the  lower  circu¬ 
lations  . 

The  stagnation  line  location  is  plotted  in  Figure  32  as  a  function 
of  advance  ratio.  It  should  be  noted  that  here  the  initial  forward  movement 
of  the  stagn  tion  line  is  the  same  for  all  propeller  rotational  speeds, 
while  in  the  case  of  rearward  motion  this  initial  forward  movement  decreased 
as  rotational  speed  increased.  The  stagnation  line  moves  rapidly  forward 
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toward  the  leading  edge  with  the  initial  increase  in  advance  ratio.  It  then 
continues  to  approach  the  leading  edge,  but  more  slowly.  After  the  stagna¬ 
tion  line  reaches  the  leading  edge,  leading  edge  separation  occurs  with 
further  increase  in  advance  ratio.  Thus  the  low  points  at  high  advance 
ratio  and  low  circulation  do  not  represent  the  location  of  the  above- 
mentioned  stagnation  line,  which  is  still  at  the  leading  edge,  but  rather 
the  location  of  the  reattachment  after  the  separation  bubble. 

CONCLUSION 

In  conclusion,  a  shrouded  propeller  has  been  successfully  employed 
on  the  XAZ-1  research  aircraft  and  on  the  U.  S.  Army  XV-11A  STOL  research 
aircraft.  Significant  increases  in  static  and  low  speed  thrust  have  been 
obtained  by  using  the  shrouded  propeller  on  these  aircraft,  with  little 
degradation  of  high  speed  performance.  Adaptation  of  the  propeller  to  the 
shroud  by  twisting  the  propeller  to  match  the  inflow  velocity  profile 
significantly  reduces  shroud  leading  edge  separation  and  improves  perfor¬ 
mance.  The  shroud  introduces  no  significant  control  and  handling  problems. 

Viscous  interaction  of  the  propeller  slipstream  with  the  external 
flow  does  indeed  cause  a  stagnation  line  to  be  located  on  the  outer  surface 
of  the  shroud  of  a  shrouded  propeller,  as  well  as  at  the  trailing  edge.  The 
forward  movement  of  this  stagnation  line  causes  a  reduction  in  circulation 
about  the  shroud  and,  consequently,  in  shroud  thrust.  Furthermore,  the 
existence  of  this  stagnation  line  on  the  outer  surface  of  the  shroud  means 
that  the  actual  thrust  will  be  less  than  that  calculated  by  the  distributed 
singularities  methods. 
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The  thrust  reduction  in  the  static  condition  and  in  rearward 


motion  can  be  correlated  with  the  distance  between  the  stagnation  line  and 
the  trailing  edge.  In  forward  motion  the  suddenness  of  the  forward  movement 
of  the  stagnation  line  precludes  such  a  correlation,  although  the  shroud 
thrust  does  decrease  as  the  stagnation  line  moves  forward. 
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FLOW  A  HOUND  SHROUD 
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STANDARD  PROPELLER 


Figure  1.  Shroud  Leading  Edge  Separation. 
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Figure  2.  Experimental  Shrouded  Propeller  on  AG-14 
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MOD/F/ED  TH/ST  PROPEL LEA  /N  SHROUD 

Figure  3.  Shroud  Leading  Edge  Separation  on  AG-14  Shroud. 
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Figure  5.  XAZ-1  MARVELETTE  Research  Aircraft. 
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Figure  10.  Propeller  Angle  of  Attack  Distribution,  XAZ-l  Shroud, 
Twisted  Propeller. 


s 


Figure  11.  Determination  of  Break-Even  Velocity  of  XAZ- 


Figure  12.  Comparison  of  Experimental  Total  Thrust  of  XAZ- 
with  Prediction  of  Simple  Momentum  Theory. 
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Figure  16.  XV-11A  in  Flight. 


Figure  17.  Inviscid  Flow  Pattern  -  Static  Condition 
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Figure  23.  Stagnation  Line  Location  -  Static  Condition. 
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Figure  24.  Pressure  Distribution  Types  -  Rearward  Motion 
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Figure  26b.  I  Low  Pattern  -  Rearward  Motion. 
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Figure  28.  Stagnation  Line  Location  -  Rearward  Motion. 
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Figure  29.  Relation  of  Shroud  Thrust  to  Stagnation  Line  Location. 
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Figure  31.  Flow  Pattern  -  Forward  Motion 
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Figure  32.  Stagnation 


Location  -  Forward  Motion. 
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ABSTRACT 

An  analytical  method  is  developed  to  determine  the 
effect  of  propeller  slipstream  on  wing  ljft,  drag  and  stability. 
The  method  is  used  to  predict  the  onset  of  stall  of  slipstream- 
immersed  wings.  The  use  of  propeller  slipstream  to  augment  the 
stability  characteristics  of  a  tilt-wing  VTOL  aircraft  during 
transition  is  also  investigated.  Finally,  a  simplified  analy¬ 
sis  is  made  of  the  effect  of  a  nonuniform  jet  velocity. 

INTRODUCTION 

The  problem  of  the  effect  of  propeller  slipstream 
on  wing  lift  was  first  investigated  by  Kcning  some  thirty 
years  ago  (Reference  1).  His  analysis  utilizes  the  assumption 
that  the  slipstream  velocity  is  only  slightly  higher  than  the 
free  stream  velocity  Other  investigators  have  extended 
Koning's  work  for  applications  to  more  practical  operating 
conditions.  In  general,  the  methods  are  relatively  complex 
and  solutions  to  specific  problems  often  require  the  aid  of 
automatic  computing  equipment. 


In  Reference  2  a  simple  analysis  is  presented  for  the 
lift  of  a  wing  spanning  a  slipstream  of  uniform  velocity.  In 
this  paper  the  analysis  is  extended  and  applied  to  practical 
V/STOL  aircraft  wing-propeller  configurations.  Correlation 
with  available  test  data  is  found  to  be  satisfactory. 

The  method  is  then  used  to  predict  the  onset  of  stall 
of  wings  immersed  in  slipstreams.  Also,  the  use  of  propeller 
slipstream  to  augment  the  stability  characteristics  of  a  tilt¬ 
wing  VTOL  aircraft  during  transition  is  investigated. 

Finally,  an  analysis  is  made  of  the  effect  on  the 
wing  lift  of  the  nonuniform  velocity  of  the  slipstream. 


PREDICTION  OF  LIFT  AND  DRAG  COEFFICIENTS 

An  expression  for  the  lift  of  a  wing  spanning  a  slip¬ 
stream  of  uniform  velocity  is  derived  in  Reference  2.  For  a 
rectangular  wing  the  total  lift  coefficient  is 


1  + 


(1) 


The  incremental  lift  coefficient  due  to  the  propeller  slip¬ 
stream  is 
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(2) 


The  spanwise  incremental  lift  distribution  assumes  the  shape  of 
the  symmetrical  curve 
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For  a  wing  extended  through  one  or  more  slipstreams 
the  tocal  lift  coefficient  is  obtained  as  the  sum  of  the  basic 
lift  coefficient  of  the  wing,  which  is  generated  by  the  wing 
when  no  slipstream  exists,  and  the  incremental  lift  coefficient 
due  to  the  effect  of  the  slipstream.  Thus, 

Q-Qw  +  aCl(|)  (3) 

In  this  equation  CL(%,  can  be  determined  by  conventional 
methods  (e.g.,  Reference  3).  The  incremental  part  aCL(-|pj 
can  be  rewritten  from  Equation  (2)  in  terms  of  thrust  coeffi¬ 
cient  as  follows: 


40,(1*)-.  0.468 - ■ — - - cos  (a,-*)  ‘os°1’  ,1  ND‘a 
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where  the  effective  angle  of  attack  of  the  slipstream  as  is 

“s“  V*-$-«u  +  H  (5) 

In  Equation  (5)  the  slipstream  deflection  angle  4>  ,  as  shown 
in  Figure  1,  is 
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and  the  empirical  factor  A  for  the  flap  deflection  angle  is 
dependent  on  the  flap  type  and  flap  chord. 

The  total  drag  coefficient  consists  of  the  following: 

(a)  Induced  drag  coefficient  of  the  basic  wing 

(b)  Induced  drag  coefficient  from  the  incremental 
lift  coefficient  due  to  the  slipstreams 

(c)  Profile  drag  coefficient  of  the  wing  areas 
immersed  in  the  slipstreams 

(d)  Profile  drag  coefficient  of  the  wing  areas  not 
immersed  in  the  slipstreams 

Thus,  the  total  drag  coefficient  is 

Ce=  CD.+ ACji+Cpos+C^  (7) 

The  induced  drag  coefficients  C0l  and  aCDi  can  be 
determined  from  the  respective  lift  coefficients  Cuw  and  aCu 
by  usual  methods.  The  profile  drag  coefficients  CDoS  and 
COOW  are  based  on  available  section  airfoil  data.  In  obtain¬ 
ing  these  data  the  angle  iw+ 4>  represents  the  angle  of  attack 
within  the  slipstream-immersed  area,  and  the  angle  iw+cfT=«w  is 
used  for  the  other  portions  of  the  wing. 

CORRELATION  WITH  EXPERIMENTAL  DATA 

As  indicated  in  Reference  2,  lift  and  drag  coeffi¬ 
cients  determined  by  the  analytical  method  presented  above  are 
in  good  agreement  with  available  experimental  data  (Reference 
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4,  5  and  6).  Additional  correlation  with  the  test  data  from 
Reference  7  is  made  in  this  section. 

In  Reference  7  test  data  are  presented  of  che  lift 
coefficient  of  a  propeller-wing  configuration  with  single- 
slotted  flaps.  By  subtracting  the  corresponding  zero  thrust 
values,  the  incremental  lift  coefficients  due  to  the  effect 
of  slipstreams  plus  the  contribution  from  the  propeller  thrust, 
aQ,5(P)  »  are  obtained.  These  are  shown  in  Figures  2  and  3  for 
flap  deflection  angles  of  0°  and  40°,  respectively.  The  solid 
lines  on  Figures  2  and  3  represent  the  calculated  values,  and 
are  seen  to  correlate  well  with  the  test  data. 

In  Reference  7  longitudinal  force  coefficient  data 
are  presented  instead  of  drag  coefficient.  Correlation  is  thus 
only  possible  for  the  longitudinal  force  coefficients.  The 
results,  as  shown  in  Figures  4  and  5,  are  seen  to  be  satisfac¬ 
tory. 

PREDICTION  OF  STALL 

The  usual  method  of  predicting  the  onset  of  wing 
stall  is  to  find  the  geometric  angle  of  attack  at  which  the 
spanwise  lift  coefficient  distribution  curve  first  reaches 
the  maximum  airfoil  section  lift  coefficient  curve  (Reference 
8). 


The  same  method  is  utilized  here  for  a  wing  immersed 
in  propeller  slipstreams.  The  wing  areas  immersed  in  the 
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slipstreams  and  those  which  are  not  are  treated  separately. 

In  this  context  it  should  be  noted  that  the  velocity  of  the 
slipstream  is  higher  than  that  of  the  free  stream.  Hence,  the 
Reynolds  number  of  the  slipstream,  and  consequently  the  max¬ 
imum  airfoil  section  lift  coefficient  Cima>(  of  the  wing 
area  inside  the  slipstream  are  higher.  Data  on  this  subject 
can  be  obtained  from  Reference  9.  Furthermore,  in  the  slip¬ 
stream-immersed  areas,  the  local  lift  coefficient  C|  should 
be  based  on  the  slipstream  dynamic  pressure,  which  is  higher 
than  the  free  stream  dynamic  pressure.  The  combined  effect  of 
the  increase  of  C\nax  and  the  decrease  of  C,  results  in  a 
delay  of  wing  stall  in  the  slipstream- immersed  areas. 

In  determining  the  onset  of  stall  in  the  wing  areas 
immersed  in  the  slipstreams,  the  direction  of  propeller  rota¬ 
tion  has  to  be  taken  into  consideration.  As  pointed  out  before, 
for  a  nonrotating  slipstream  the  incremental  lift  coefficient 
distribution  curve  is  symmetrical  with  respect  to  the  slip¬ 
stream  centerline.  However,  as  a  result  of  the  slipstream 
rotation,  the  local  angle  of  attack  is  increased  on  that 
portion  of  the  wing  in  the  propeller  upwash  and  vice  versa, 
thereby  shifting  the  peak  of  the  incremental  lift  coefficient 
distribution  curve  away  from  the  centerline.  This  phenomenon 
is  noticed  in  the  experiments  as  reported  in  Reference  10. 

A  numerical  analysis  was  performed  on  the  lift  coeffi¬ 
cient  distribution  of  the  wing  model  used  in  Reference  5. 

Shown  in  Figure  7  is  the  lift  coefficient  distribution  curve 
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with  the  propeller  rotating  up  at  wing  tip.  In  Figure  7  the 
corresponding  result  is  shown  for  the  propeller  rotating  down 
at  wing  tip  In  both  cases,  due  to  a  comparatively  low  angle 
of  attack  no  stall  occurs  within  the  slipstream-immersed 
areas . 


When  the  wing  geometric  angle  of  attack  is  higher, 
the  situation  becomes  somewhat  different.  In  Figure  8,  where 
the  propeller  blades  rotate  up  at  the  vTing  tip,  stall  still 
does  not  occur  within  the  slipstreams.  On  the  other  hand, 
with  the  propeller  blades  rotating  down  at  the  wing  tip,  the 
peak  of  the  C,  curve,  which  is  now  inboard  of  the  propeller 
centerline,  exceeds  the  local  Cimax  .  As  a  result  the  slip¬ 
stream-  immersed  areas  are  partially  stalled,  as  shown  in 
Figure  9. 


Test  data  of  Reference  6  are  used  to  verify  the  above 
mentioned  method  of  predicting  the  onset  of  wing  stall.  The 
results  are  presented  in  Figure  10.  The  experimental  data 
shown  correspond  to  the  wing  angles  of  attack  where  the  lift 
curve  slope  begins  to  change  noticeably.  It  is  seen  that  for 
constant  values  of  CTS  ,  up  to  CT.S  -0.5  the  calculated  values 
of  the  wing  angle  of  attack  at  initial  stall  within  the  slip¬ 
streams  correlate  well  with  test  data. 

EFFECT  ON  V/STOL  AIRCRAFT  STABILITY 

The  feasibility  of  utilizing  the  slipstream  for 
stability  augmentation  was  investigated  for  a  4-propeller 
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tilt-wing  VTOL  aircraft  having  a  variable-angle  flap. 

An  analog  computer  analysis  was  conducted  of  the 
longitudinal  stability  characteristics  for  transition  flight 
speeds  of  30  and  70  knots.  The  results  are  as  follows: 

Computer  Solution  At  30  Knots 

The  response  of  the  unstabilized  aircraft  to  a  flap 
pulse  input  is  shown  in  Figure  11.  It  is  noted  that  the 
response  consists  of  an  unstable  aperiodic  mode  with  a  time 
to  double  of  0.9  seconds. 


The  variable-angle  flap  stabilizer  system  investigated 
is  described  by  the  following  equation: 

b  +  k,  9  +  k3w  =  0  (8) 


where 


6  is  flap  stabilizer  deflection 
6  is  fuselage  attitude 
w  is  normal  perturbation  velocity 
k,  is  stabilizer  attitude  gain 
k3  is  stabilizer  normal  velocity  gain 


The  effect  of  the  constants  and  k3  were  investiga¬ 
ted  and  it  was  found  that  the  aircraft  becomes  very  stable  with 
the  stabilizer  setting  as  described  by 

<5-ise- o.32w=  o  (9) 
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The  transient  response  of  the  aircraft  with  this  stabilizer 
setting  is  shown  in  Figure  12. 

Computer  Solution  At  70  Knots 

The  basic  unstabilized  aircraft's  response  to  a  flap 
pulse  input  is  shown  in  Figure  13.  As  expected,  the  response 
for  this  flight  condition  is  similar  to  that  of  a  typical 
fixed-wing  aircraft.  The  aircraft  response  consists  of  a  long 
period  oscillation  (phugoid  mode)  with  a  period  of  10  seconds 
and  a  time  to  half  amplitude  of  9  seconds. 

The  use  of  a  flap  stabilizer  corresponding  to 

<5- 15  0-0.169  w«  0  (10) 

results  in  the  response  shown  in  Figure  14.  As  noted  from 
this  figure,  the  response  is  very  rapid  and  stable. 

EFFECT  OF  NONUNIFORM  VELOCITY  PROFILE 

In  the  analysis  presented  in  Reference  2,  the  slip¬ 
stream  is  assumed  to  be  of  uniform  velocity.  This  theory  is 
now  extended  in  order  to  determine  the  gross  effect  of  the 
nonuniform  velocity  profile  on  the  wing  lift  characteristics. 
For  this  purpose  a  slipstream  with  two  concentric  velocity 
zones  is  considered.  The  velocity  pattern  is  either  as  in  case 
(a)  or  case  (b)  of  Figure  15.  In  the  outer  ring  there  is  a 
perturbation  velocity  potential  §,  ,  and  in  the  central  region 
the  potential  is  •  The  potential  in  the  free  stream  is 

defined  as  .  These  potentials  are  solutions  of  the  two- 
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dimensional  Laplace's  Equation  with  the  following  boundary 
conditions : 


rs  =  r  =  rs 


o  r  <  r 
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The  potentials  at  the  wing  trailing  edge  are  found  to 
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The  derivation  of  these  expressions  is  presented  in  full  detail 
in  Reference  11. 


The  lift  coefficient  of  the  wing  spanning  a  slipstream 
with  two  velocity  zones  based  on  a  free  stream  dynamic  pressure 
becomes 

Cl” +  (13) 

Tt  Jo 

A  numerical  investigation  was  made  to  compare  the 
total  lift  coefficient  produced  by  slipstreams  of  the  same 
thrust  but  with  different  velocity  profiles  as  follows: 

(a)  Velocity  in  the  outer  .  _ng  is  lower  than  that  in 
the  central  region 

(b)  Velocity  in  the  outer  ring  is  higher  than  that  in 
the  central  region 
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(c)  Uniform  velocity 

It  was  found  that  the  total  lift  coefficient  in  case  (a)  is 
the  highest.  The  result  of  this  simplified  analysis  is  quali 
tatively  substantiated  by  the  test  data  of  Reference  12. 
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FIGURE  I:  DEFLECTION  OF  PROPELLER  SLIPSTREAM  DUE  TO  FLIGHT 
SPEED 
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O  TEST  DATA 
-  CALCULATED  VALUES 


°*  w,  DEG.  “w,  DEG. 

(a)  CjS=  1.0  (b)CTS=0.95 


FIGURE  2:  CORRELATION  OF  PREDICTED  INCREASE  OF  LIFT  COEFFICIENT 
DUE  TO  PROPELLER  THRUST  AND  SLIPSTREAM  WITH  TEST 
DATA  («5f  =  0°) 
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(a)  Cj^3=  1.0 


(b)  Cj g=  0.95 


FIGURE  3:  CORRELATION  OF  PREDICTED  INCREASE  OF  LIFT  COEFFICIENT 
DUE  TO  PROPELLER  THRUST  AND  SLIPSTREAM  WITH  TEST 
DATA  (if  =  40°) 
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FIGURE  6:  PREDICTED  SPANWISE  LIFT  COEFFICIENT  DISTRIBUTION  AND 
EXTENT  OF  WING  STALL  (PROPELLER  ROTATING  UP  AT 
WING  TIP,  =15  DEGREES,  CT(S=0.58,  /R=  5.00,  NACA  0015 
AIRFOIL  SECTION). 
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STALLED  AREA 

FIGURE  7:  PREDICTED  SPANWISE  LIFT  COEFFICIENT  DISTRIBUTION  AND 

EXTENT  OF  WING  STALL  (PROPELLER  ROTATING  DOWN  AT  WING 
TIP,  =15  DEGREES,  CT>S=0.58,  /*?=5.00,  NACA  0015 
AIRFOIL  SECTION). 


FIGURE  8:  PREDICTED  SPANWISE  LIFT  COEFFICIENT  DISTRIBUTION  AND 
EXTENT  OF  WING  STALL  (PROPELLER  ROTATING  UP  AT  WING 
TIP,  ~w  =50  DEGREES,  CTS=0.88,  >6=5.0,  NACA  0015 
AIRFOIL  SECTION). 
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STALLED  AREA 


FIGURE  9:  PREDICTED  SPANWISE  LIFT  COEFFICIENT  DISTRIBUTION  AND 
EXTENT  OF  WING  STALL  (PROPELLER  ROTATING  DOWN  AT 
WING  TIP,~w  =  50  DEGREES,  CTjS=0.88,  /f?=  5.0,  NACA  0015 
AIRFOIL  SECTION) 


WING  INITIAL  STALL  ANGLE 


THRUST  COEFFICIENT 


FIGURE  10:  CORRELATION  OF  PREDICTED  WING  INITIAL  STALL  WITHIN 
SLIPSTREAM  WITH  TEST  DATA  (PROPELLER  ROTATING 
DOWN  AT  WING  TIP,4?=5.0,  NACA  0015  AIRFOIL  SECTION), 
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FIGURE  13:  TRANSIENT  h  TSPONSE  OF  THE  UNSTABILIZED  AIRCRAFT 
TO  A  FLAP  CONTROL  INPUT  (V0  =  70  KNOTS). 
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FIGURE  14:  TRANSIENT  RESPONSE  WITH  NORMAL  VELOCITY  AND 
ATTITUDE  FEEDBACK  STABILIZATION  (Vn  =  70  KNOTS). 


FIGURE  15:  VELOCITY  ZONES  AND  PERTURBATION  POTENTIALS  IN 
SLIPSTREAM 
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velocity 

yu,  ratio  of  free  stream  velocity  to  slipstream 
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Mz  ratio  of  free  stream  velocity  to  slipstream 

velocity  in  central  region 
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forward  speed 

free  stream  perturbation  velocity  potential 
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$2  slipstream  perturbatio:  velocity  potential  in 
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4>Si  slipstream  perturbation  velocity  potential  at  wing 
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$S2  slipstream  perturbation  velocity  potential  at  wing 
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AERODYNAMIC  PROPERTIES  OF  AIRFOILS 
IN  NONUNIFORMLY  SHEARED  FLOWS* 

by 

W.  G.  Brady  and  G.  R.  Ludwig 

Cornell  Aeronautical  Laboratory,  Inc. 
Buffalo,  New  York,  U.S.A. 


SUMMARY 

A  theoretical  method,  which  requires  the  use  of  a  digital  computer, 
was  developed  to  predict  pressure  distributions  on  airfoils  in  two-dimensional 
nonuniformly  sheared  flows.  The  theory  is  applicable  to  airfoils  of  arbitrary 
profile,  including  those  with  flaps  and  other  auxiliary  lift  devices.  To  verify 
this  theory,  aerodynamic  characteristics  of  an  airfoil  were  investigated  both 
theoretically  and  experimentally  in  a  simple  nonuniformly  sheared  flow 
consisting  of  two  segments  with  shears  of  equal  magnitude  but  of  opposite 
sign.  Agreement  between  computed  and  experimental  pressure  distributions 
was  at  least  as  good  as  the  agreement  obtained  between  exact  potential  flow 
theory  and  experimental  data  for  an  unsheared  flow. 

Previous  investigations  have  demonstrated  experimentally  that  the 
lift  characteristics  of  an  airfoil  in  a  two-dimensional  nonuniformly  sheared 
jet,  similar  to  a  section  through  a  propeller  slipstream,  and  in  an  axially 
symmetric  nonuniformly  sheared  jet  similar  to  a  propeller  slipstream  are 
strongly  dependent  on  the  airfoil's  location  in  the  flow.  On  the  basis  of  the 
present  theoretical  and  experimental  results,  a  mechanism  is  postulated  by 
which  the  large  changes  in  lift  characteristics  can  occur.  This  mechanism 
suggests  that  the  presence  of  nonuniform  shear  in  a  propeller  slipstream 
must  be  considered  if  accurate  pressure  distributions  and,  hence,  lift  and 

This  paper  is  based  on  work  performed  for  the  U.  S„  Army  Aviation  Materiel 
Laboratories,  Fort  Eustis,  Virginia. 


stalling  characteristics  of  an  airfoil  immersed  in  a  propeller  slipstream 
with  significant  shear  are  to  be  obtained. 
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Function  defined  by  Equations  (1). 

Displacement  of  dividing  streamline  (see  Figure  2). 

Vertical  height  of  airfoil  midchord  above  tunnel  centerline, 
feet. 


Smoothing  parameter  used  in  computer  program. 

Vertical  distance  from  tunnel  centerline  to  wall  (see  Figure  2). 
Static  pressure  at  (x,  y),  pounds  per  square  foot. 

Static  pressure  in  free  stream,  pounds  per  square  foot. 
Uniform  flow  free-stream  velocity,  feet  per  second. 

Sheared  flow  free-stream  velocity,  feet  per  second. 
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Sheared  flow  free-stream  velocity  at  y  -  h  ,  feet  per  second. 

Sheared  flow  free -stream  velocity  of  stagnation  streamline, 
feet  per  second. 

Free-stream  velocity  at  lower  wind-tunnel  wall,  feet  per 
second. 

Free-stream  velocity  at  upper  wind-tunnel  wall,  feet  per 
second. 

Free-stream  velocity  along  vertical  centerline  in  wind  tunnel, 
feet  per  second. 

Average  free-stream  velocity  in  wind-tunnel  test  section, 
feet  per  second. 

Velocity  component  in  x  -direction,  feet  per  second. 

Velocity  component  in  y -direction,  feet  per  second. 

Coordinate  along  axis  parallel  to  wind-tunnel  walls,  positive 
downstream;  x  =  0  at  airfoil  midchord. 

Coordinate  perpendicular  to  X  -axis,  positive  up;  y  =  0  at 
wind-tunnel  centerline. 

Airfoil  angle  of  attack,  positive  leading  edge  up,  degrees. 

Airfoil  angle  of  attack  for  zero  lift. 

Airfoil  angle  of  attack  measured  from  zero  lift  angle  of 
attack;  os  =  cc-Oc.jz0  . 

Vortex  sheet  singularity  strength  of  airfoil  aft  lower  surface. 
Vortex  sheet  singularity  strength  of  airfoil  aft  upper  surface. 
Vortex  sheet  singularity  strength  of  airfoil  leading  edge  surface. 
Vortex  sheet  singularity  strength  of  lower  wind-tunnel  wall. 
Vortex  sheet  singularity  strength  of  upper  wind-tunnel  wall. 

Air  density,  assumed  constant,  slugs  per  cubic  foot. 

Stream  function;  see  Equation  (1). 
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INTRODUCTION 


Many  of  the  analytical  problems  of  wing  aerodynamics  in  flows  at 
low  subsonic  speeds  are  treated  successfully  within  the  framework  of  potential 
flow  theory.  However,  potential  flow  theory  is  inadequate  for  predicting 
analytically  the  interaction  between  a  wing  and  a  propeller  slipstream  when 
the  slipstream  intersects  the  wing  and  the  flow  in  the  propeller  slipstream 
is  not  irrotational,  i.e.  ,  when  the  slipstream  flow  has  significant  shear.  An 
additional  difficulty  which  gives  rise  to  nonlinearities  in  the  analytical  treat¬ 
ment  results  from  the  interaction  of  the  free  boundaries  of  the  slipstream 
flow  and  the  wing. 

A  number  of  analytical  treatments  of  the  aerodynamics  of  sheared 
flows  have  been  published;  typical  are  References  1  through  4.  Exact  aero¬ 
dynamic  theories  for  uniform  two-dimensional  shear  (constant  vorticity)  are 
available  (References  1  and  4),  but  the  available  solutions  for  nonuniform 
shear  (Reference  3,  for  example)  are  not  particularly  useful  to  the  aero- 
dynamicist  because  of  limitations  of  the  theoretical  assumptions.  These  are 
primarily  restrictions  to  small  shear  and  shear  gradient.  Three-dimensional 
sheared  flow  theory  (Reference  2,  for  example)  has  hardly  advanced  beyond 
attempts  to  arrive  at  a  formulation  which  is  mathematically  tractable. 

This  paper  reports  the  latest  results  of  an  experimental  and 
theoretical  program  which  has  been  concerned  with  the  aerodynamic  char¬ 
acteristics  of  airfoils  in  sheared  flows.  This  work  and  the  earlier  work 
(References  4  through  8)  were  conducted  for  the  U.  S.  Army  Aviation 
Materiel  Laboratories  (formerly  U,  S.  Army  Transportation  Research 
Command).  A  point  of  particular  interest  from  the  earlier  studies  was  the 
discovery  that  nonuniform  shear  in  a  two-dimensional  flow  could  have  a 
marked  effect  on  airfoil  stall  characteristics  (Reference  6).  It  is  noted  in 
Reference  6  that  in  the  specific  two-dimensional  nonuniformity  sheared  flow 
treated,  large  variations  in  maximum  lift  were  obtained  near  the  flow 
centerline  when  the  airfoil  vertical  position  was  varied  by  distances  of  the 
order  of  magnitude  of  the  airfoil  thickness.  The  current  program  was 
designed  as  the  initial  step  toward  investigating  the  reasons  for  this  unusual 
behavior. 
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Since  the  maximum  lift  and  stall  behavior  are  dependent  on  the 
properties  of  the  boundary  layer  on  the  airfoil,  and  these  properties,  in 
turn,  are  primarily  influenced  by  the  pressure  distribution,  it  was  decided 
that  any  theoretical  investigation  must  be  one  which  will  produce  accurate 
pressure  distributions  on  an  airfoil  in  a  nonuniformly  sheared  flow.  Of 
course,  the  unusual  behavior  near  stall  may  have  been  caused  by  other  effects 
such  as  a  dir  :nteraction  of  the  boundary  layer  with  the  inviscid  sheared 
flow  in  such  a  ..  nner  as  to  promote  or  delay  boundary-layer  separation,  or 
by  subtle  wind-tunnel  wall  effects  associated  with  the  sheared  flow.  The 
direct  interaction  possibility  has  been  considered  previously  with  essentially 
negative  results  (Reference  6).  The  effect  of  wind-tunnel  walls  could  be 
determined  by  including  them  in  the  analysis,  which,  if  successfully  completed, 
could  be  used  to  investigate  wind-tunnel  wall  effects  by  varying  the  wall 
positions. 

In  keeping  with  the  above  observations,  a  digital  computer  program 
was  developed  to  compute  pressure  distributions  on  a  two-dimensional  airfoil 
in  a  particular  inviscid,  nonuniformly  sheared  flow.  The  sheared  flow 
velocity  profile  was  kept  as  simple  as  possible  for  this  initial  development, 
while  still  incorporating  what  was  believed  to  be  all  essential  features  of  a 
nonuniformly  sheared  flow.  The  free -stream  velocity  profile  was  such  that 
shears  of  equal  magnitude  but  opposite  sign  occurred  above  and  below  the 
centerline  of  the  flow.  Wind-tunnel  walls  were  also  included  in  the  theoretical 
program  so  that  the  flow  model  would  duplicate  the  experiment  as  nearly  as 
possible.  At  the  same  time,  experimental  airfoil  pressure  distribution  data 
were  obtained  in  a  wind  tunnel  for  a  free -stream  velocity  profile  nearly  the 
same  as  that  used  in.  the  theory.  These  data  were  obtained  to  enable  the 
accuracy  of  the  theoretical  program  to  be  verified. 

Once  the  analytical  technique  is  established  as  to  its  accuracy  on 
the  basis  of  the  nonuniformly  sheared  flow  treated  in  the  present  program, 
its  extension  to  the  more  complicated  nonuniformly  sheared  flow  of 
Reference  6  is  justified.  The  analytical  technique  would  then  provide  a 
powerful  tool  for  the  investigation  of  flow  mechanisms  responsible  for  the 
unusual  lift  behavior  near  stall  shown  in  Reference  6. 
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THEORETICAL  PROGRAM 


The  main  characteristic  of  a  sheared  flow  which  distinguishes  it 
from  a  nonshea  red  flow  is  the  fact  that  it  is  rotational.  The  powerful  tech¬ 
niques  of  potential  theory,  including  conformal  mapping,  that  can  be  used 
for  two-dimensional  irrotational  flow  problems,  are  generally  not  applicable 
to  rotational  flows. 

Rather  than  Laplace's  equation,  the  stream  function  must  satisfy 

V V  =  /(>/  (Eq.  1) 

i.e.  ,  Vz*p  is  a  constant  along  a  streamline  (along  which  the  stream  function 
is  also  a  constant).  It  can  be  shown  that,  in  two-dimensional  flow,  this 
constant  is  the  vorticity;  that  is,  the  vorticity  is  a  constant  along  a  stream¬ 
line  in  an  inviscid  two-dimensional  rotational  flow.  Except  for  the  case  of 
uniform  shear  (  -f(P)  =  constant),  or  the  exceptional  case  ■fl'f)  -  cup  , 

O >  =  constant,  one  is  faced  with  solving  a  nonlinear  partial  differential 
equation.  However,  if  the  free-stream  velocity  profile  can  be  approximated 
by  piecewise  linear  segments,  as  in  Figure  1,  then  the  streamlines  passing 
through  those  points  at  which  the  velocity  gradient  changes  (points  (a),  (b), 
and  (c)  in  Figure  1,  for  example)  separate  regions  of  constant  vorticity 
throughout  the  flow.  This  approximation  can  be  used  to  change  the  problem 
of  solving  a  nonlinear  partial  differential  equation  with  known  boundary  con¬ 
ditions  to  one  with  a  linear  equation  {  Vz<p  =  constant)  within  regions  with 
free  (and,  hence,  unknown)  boundaries. 

The  above  approach  was  used  to  formulate  a  theory  for  finding  the 
pressure  distribution  on  any  airfoil  at  angle  of  attack  in  an  inviscid,  incom¬ 
pressible,  two-dimensional  nonuniforrnly  sheared  flow.  A  digital  computer 
program  was  written  for  a  two-dimensional  velocity  profile  like  that  of 
Figure  1(b).  This  velocity  profile  has  the  advantage  of  being  one  of  the 
simplest  nonuniforrnly  sheared  flows  to  which  the  theory  can  be  applied. 
Application  of  the  theory  in  its  development  stage  to  the  more  complicated 
sheared  flow  used  in  the  experimental  work  of  Reference  6  (which  included 
the  added  complexity  of  free  jet  boundaries)  was  felt  to  be  premature. 
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Wind-tunnel  walls  were  incorporated  into  the  program  so  that  the  computed 
results  would  be  directly  comparable  to  the  results  of  the  experimental 
program,  and  to  allow  the  possibility  of  investigating  wind-tunnel  wall  effects 
at  a  later  date.  The  theoretical  program  is  not  limited  to  small  magnitudes 
of  the  shear  or  to  small  flow  disturbances.  This  latter  feature  is  essential 
if  inviscid  pressure  distributions  are  to  be  obtained  for  airfoil  angles  of 
attack  near  stall. 

Details  of  the  theory  and  its  implementation  on  an  IBM  7044  digital 
computer  are  presented  in  Reference  9.  A  brief  outline  of  the  basic  ideas 
incorporated  and  certain  numerical  techniques  of  importance  to  the  program 
are  discussed  in  the  following  paragraphs. 

Description  of  Computer  Program 

The  two-dimen:: ional  flow  model  assumed  for  the  computer  program 
is  shown  in  Figure  2.  A  Cartesian  coordinate  system  is  used.  The  x.  -axis 
is  parallel  to  the  wind-tunnel  walls  and  positive  in  the  direction  of  the  flow. 
The  origin  is  centered  vertically  between  the  walls  and  x  =  0  is  located  at 
the  airfoil  midchord.  The  wind-tunnel  walls  and  the  airfoil  surface  are 
represented  by  bound  vortex  sheets  of  varying  strength.  The  free-stream 
velocity  variation  in  the  y  -direction  is  piecewise  linear  and  symmetric 
about  y  -  0  (point  (a)  in  Figure  2).  On  either  side  of  the  streamline  which 
passes  through  point  (a)  (denoted  henceforth  as  the  "dividing  streamline") 
the  distributed  vorticity  is  constant.  If  the  assumed  free-stream  velocity 
profile,  V(<j)  ,  is  symmetric  (  U,  -  U2  in  Figure  2),  as  is  the  case  in  the 
work  presented  here,  the  vorticity  on  either  side  of  the  dividing  streamline 
is  equal  in  magnitude  but  of  opposite  sign.  If  the  airfoil  is  of  nonsymmetric 
profile,  or  is  nonsymmetrically  disposed  in  the  wind  tunnel,  the  resulting 
disturbance  flow  will  distort  the  dividing  streamline  in  some  fashion.  This 
distorted  shape  is  denoted  by  y(Z)  in  Figure  2.  Equation  (1)  states  that 
the  vorticity  in  a  rotational  flow  is  constant  along  streamlines;  therefore, 
the  dividing  streamline  still  separates  regions  of  constant  vorticity.  If  the 
shape  of  the  dividing  streamline  g(x)  were  known,  then  it  would  be  possible 
to  compute  the  variation  in  strength  of  the  bound  vortex  sheets  representing 
the  wind-tunnel  walls  and  airfoil  surface. 
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If,  for  the  present,  it  is  assumed  that  g(z)  is  known,  then  the 
procedure  used  in  calculating  the  singularity  strength  distributions  on  the 
airfoil  surface  and  wind-tunnel  walls  is  as  follows.  The  surfaces  are 
represented  by  small  segments  of  vortex  sheets  each  of  which  has  a  constant 
singularity  strength.  The  integrals  over  segments  of  the  vortex  sheets  for 
use  in  the  equations  for  velocity  components  can  then  be  evaluated  in  closed 
form.  Hence,  the  velocity  components  are  expressible  as  finite  series  in 
terms  of  the  vortex  sheet  singularity  strengths.  The  vortex  sheets 
representing  the  wind-tunnel  walls  are  assumed  to  have  known  constant 
strengths  far  upstream  and  downstream  of  the  airfoil.  The  boundary 
conditions  for  no  normal  flow  across  the  wind-tunnel  walls  and  airfoil 
surface  are  next  applied  at  the  center  of  each  of  the  vortex  sheet  segments. 

If  these  boundary  conditions  were  the  only  conditions  to  be  satisfied,  there 
would  result  a  determinant  set  of  linear  simultaneous  algebraic  equations. 
However,  in  addition,  the  Kutta- Joukowski  condition  at  the  airfoil  trailing 
edge  must  be  satisfied.  This  constraint  is  implemented  in  the  computer 
program  by  requiring  that  the  sum  of  the  singularities  representing  the 
upper  and  lower  airfoil  surfaces  vanish  in  the  limit  as  the  trailing  edge  is 
approached.  This  results  in  an  additional  equation,  and  the  problem  is  now 
overdetermined.  The  overdeterminancy  results  from  the  choice  of  constant 
singularity  strength  for  each  vortex  sheet  segment.  If  it  had  been  assumed 
that  the  singularity  strength  of  each  segment  varied  linearly,  there  would 
have  been  as  many  equations  as  unknowns.  The  constant  strength  assumption 
was  made  to  limit  the  number  of  linear  equations  requiring  solution. 

The  method  used  to  solve  the  overdetermined  set  of  equations 
followed  that  of  Reference  10  which  presents  a  numerical  technique  for 
obtaining  smoothly  varying  solutions  to  such  sets  when  it  is  known  beforehand 
that  the  solutions  should  be  smooth.  The  smoothing  technique  is  required 
to  suppress  oscillations  arising  from  small  errors  in  the  numerical  computa¬ 
tions.  These  oscillations  can  become  very  large  as  is  illustrated  in  Figure  3 
where  computed  pressure  distributions  on  a  symmetric  Joukowski  airfoil  in 
a  uniform  stream  are  compared  to  the  exact  solution  obtained  through 
conformal  mapping.  It  is  evident  that  smoothing  is  required  and  that  the 
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process  works  quite  well,  at  least  for  the  present  model.  The  maximum 
deviation  Detween  the  smoothed  numerical  solution  and  the  exact  solution 
occvirred  at  the  suction  peak  and  amounted  to  a  3  percent  overshoot.  The 
size  of  the  vortex  sheet  segments  chosen  to  represent  the  airfoil  surface 
in  this  region  may  have  contributed  to  this  small  error.  Hence,  it  would 
appear  that  the  process  used  to  calculate  the  singularity  strength  dis¬ 
tribution  on  the  surfaces  provides  satisfactory  accuracy. 

An  important  part  of  the  computer  program  is  an  iterative  technique 
for  determining  the  function  q(z)  numerically.  Mass-flow  continuity  is 
the  basis  for  this  iterative  technique  in  the  developed  computer  program.  If 
an  initial  shape  q0(z)  is  assumed,  the  strengths  of  the  singularity  distribu¬ 
tions  (vortex  sheets  representing  the  wind-tunnel  walls  and  the  airfoil 
surface)  can  be  determined  as  outlined  above.  Once  these  distributions  are 
known,  the  first  iterated  streamline  shape,  qt(z)  ,  is  determined  by 
computing,  at  each  value  of  x  considered,  the  vertical  point  (  g* '  (x)  ,  say) 

in  the  flow  where  the  mass  flow  is  equal  to  that  on  either  side  of  the  dividing 
streamline  in  the  undisturbed  free  stream,  and  then  by  taking  the  average 
of  q0(%)  and  q*(x )  .  Thus 

*  {[ 9o(*)  +  gH*)  • 

Given  q,  (x)  ,  again  the  singularity  strength  distributions  for  the  bound 

vortex  sheets  representing  the  wind-tunnel  walls  and  the  airfoil  can  be 
obtained,  from  which  the  second  iterated  streamline  shape,  g2(Z)  ,  can  be 
determined.  If  the  iterative  process  converges,  Um  a  (x)-q(z).  Hopefully, 

L  ^  OO 

the  rate  of  convergence  is  such  that  this  limit  is  very  nearly  attained  after 
only  a  few  iterations.  For  the  configurations  run  thus  far,  five  to  six 
iterations  appeared  to  be  sufficient  to  provide  an  acceptable  solution. 

The  singularity  strength  distributions  required  to  satisfy  the  boundary 
conditions  at  the  wind-tunnel  wails  and  the  airfoil  surface  are  those  which 
correspond  to  the  g(x)  evaluated  in  the  final  iteration.  The  flow  field  is 
thus  completely  determined.  Velocities  on  the  airfoil  surface  can  then  be 
converted  to  static  pressures  by  means  of  the  Bernoulli  equation  for 
rotational  flow, 
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(Eq.  2) 


YP(u.Z+vZ)+p  -  |  p[u(c/)Y  +  p0  , 

applied  along  the  airfoil  stagnation  streamline.  The  proper  ’’alue  of  U(y) 
in  the  freestream  to  be  used  in  Equation  (2)  is  determined  by  computing  the 
mass  flow  between  the  wind-tunnel  wall  and  the  airfoil.  The  value  of  U(^) 
which  corresponds  to  this  mass  flow  in  the  undisturbed  freestream  is  the 
correct  value.  It  should  be  noted  that  the  use  of  mass  flow  continuity  in  the 
iterative  process  and  in  determining  U(y)  requires  a  knowr.  reference 
streamline  of  the  flow.  In  the  present  case,  a  wind-tunnel  wall  provides 
this  reference  streamline. 

Using  the  program  outlined  above,  pressure  distributions  on  a 
symmetric  Joukowski  airfoil  at  angle  of  attack  were  computed  for  several 
locations  in  the  simple  nonuniformly  sheared  flow.  The  results  are  compared 
to  experiment  in  a  later  section. 

EXPERIMENTAL  PROGRAM 

The  experimental  program  was  designed  to  provide  airfoil  pressure 
distribution  data  by  means  of  which  the  accuracy  of  the  theoretical  results 
could  be  determined.  At  the  same  time,  additional  experimental  evidence 
regarding  airfoil  maximum  lift  behavior  in  sheared  flows  would  be  obtained. 

The  experiments  were  made  in  the  subsonic  leg  of  the  CAL  One- 
Foot  High-Speed  Wind  Tunnel.  This  leg  of  the  wind  tunnel  has  a  test  section 
with  a  cross  section  of  17  inches  by  24  inches  and  is  operated  as  a  closed- 
throat  nonreturn  tunnel.  The  two-dimensional  airfoil  used  in  this  research 
has  a  symmetric  Joukowski  profile  with  a  thicknes s -chord  ratio  of  17  percent 
and  a  chord  of  6  inches.  The  wing  was  mounted  in  the  wind-tunnel  test 
section  by  a  supporting  structure  which  permitted  independent  change  in 
geometric  angle  of  attack  and  vertical  location  of  the  wing  in  the  wind  tunnel 
by  means  of  manual  adjustment. 

Instrumentation  was  relatively  simple.  Flow  calibrations  in  the 
wind-tunnel  test  section  were  performed  with  a  conventional  3/ 16-inch- 
diameter  pitot-static  probe  which  included  static  pressure  taps  from  which 
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flow  angularity  data  were  obtained.  Static  and  dynamic  pressures  were 
measured  with  manometers.  Thirty-two  static  pressure  taps  were  distri¬ 
buted  over  the  top  and  bottom  surfaces  of  the  wing  at  a  section  0.7  5  inch  to 
one  side  (left  looking  downstream)  of  the  wind-tunnel  centerline.  These  taps 
were  connected  to  an  inclined  manometer  bank  and  the  resulting  pressures 
were  recorded  photographically.  Free- stream  static  pressure  was  obtained 
from  a  pressure  tap  on  the  ceiling  of  the  test  section.  This  tap  was  approxi¬ 
mately  2-1/2  chord  lengths  upstream  of  the  model  midchord. 

The  sheared  flow  in  the  wind-tunnel  test  section  was  produced  by 
a  screen  placed  slightly  more  than  3  feet  (between  6  and  7  wing  chords) 
upstream  of  the  wing.  This  screen  consisted  of  an  array  of  circular  rods 
spanning  the  wind-tunnel  section  horizontally  and  secured  by  a  frame  which 
was  clamped  between  two  sections  of  the  wind-tunnel  circuit.  The  spacings 
between  rods  and  the  rod  diameters  were  varied  so  as  to  introduce  variable 
losses  across  the  flow.  By  proper  spacing  and  rod  size  variation,  the 
vertical  distribution  of  losses  in  the  flow  at  the  screen  can  be  such  that  the 
desired  sheared  flow  velocity  contour  is  obtained  in  the  test  section. 

The  theory  presented  in  Appendix  II  of  Reference  5  was  used  in  the 
design  of  the  screen.  Application  of  this  theory  is  somewhat  laborious  and 
the  results  are  dependent  on  the  accuracy  of  the  empirical  relationship  between 
local  screen  solidity  and  loss  coefficient.  Apparently  the  empirical  curve 
presented  in  Reference  5  is  not  reliable  at  high  values  of  solidity  because 
the  use  of  this  curve  provided  a  velocity  distribution  which  was  considerably 
different  from  the  design  velocity  distribution  near  the  centerline  of  the 
test  section.  The  final  configuration  used  for  the  shear  screen  was  tailored 
by  a  trial  and  error  process  of  small  modifications  to  the  theoretical  design 
until  it  was  judged  that  further  improvement  could  not  be  obtained  without 
undue  expenditure  of  time. 

A  number  of  velocity  profiles  were  measured  in  the  test  section  in 
a  vertical  plane  at  the  model  axis  of  rotation  (midchord  of  the  wing)  using 
the  3/ 1  6-inch-diameter  pitot-static  probe.  Figure  4  presents  velocity  dis¬ 
tributions  for  the  final  shear  screen  configuration  used  during  pressure 
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data  tests.  It  can  be  seen  that  the  velocity  distribution  approached  the 
design  distribution  with  reasonable  accuracy  but  that  seme  irregularities 
were  present  in  the  experimental  data  in  the  vicinity  of  5  inches  above  and 
below  the  tunnel  centerline.  These  ir regularities  are  not  believed  to  be  of 
importance  because  they  were  not  large  and  because  the  test  region  only 
extended  to  3  inches  above  and  below  the  centerline. 

In  addition  to  the  velocity  calibrations  of  the  sheared  flow,  cali¬ 
brations  of  flow  angularity  in  the  vertical  plane  were  obtained  using  the 
static  pressure  angularity  taps  on  the  3/ 16-inch  probe.  The  variation  of 
flow  angularity  over  the  test  section  height  was  small.  A  total  variation  of 
0.65  degrees  was  observed.  The  geometric  angle  of  attack  of  the  airfoil 
was  corrected  for  this  small  variation  in  flow  angularity. 

The  wind-tunnel  airfoil  tests  were  of  three  types:  uniform  flow, 
low-turbulence  tests;  uniform  flow,  high-turbulence  tests;  and  sheared 
flow  tests.  The  uniform-flow  tests,  with  and  without  high  turbulence,  were 
designed  as  an  experimental  control.  The  shear  screens  used  to  generate 
the  sheared  flow  also  generate  higher  turbulence  levels  in  the  flow.  It  was 
possible  that  this  higher  turbulence  level  would,  in  itself,  have  some 
marked  effect  on  airfoil  separation  characteristics.  The  low  and  high 
turbulence  tests  in  uniform  flow  were  to  resolve  this  possibility.  The 
uniform  flow,  low  turbulence  tests  were  performed  in  the  clear  wind  tunnel 
and  the  uniform  flow,  high  turbulence  tests  were  performed  with  screens 
of  uniform  solidity  mounted  upstream  from  the  test  section.  These  screens 
were  similar  to  the  shear  screen  except  that  the  horizontal  rods  used  were 
uniformly  spaced.  Two  such  turbulence  screens  were  used  in  separate 
tests,  one  with  rod  diameter  and  spacing  corresponding  to  that  portion  of  the 
shear  screen  with  low  solidity  (relatively  large  open  area  ratio)  and  one 
screen  corresponding  to  the  high  solidity  portion  of  the  shear  screen. 
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EXPERIMENTAL  RESULTS 


Uniform  Flow  Tests 

In  Figure  5,  pressure  distributions  measured  on  the  airfoil  for 
several  conditions  of  free-stream  velocity  and  turbulence  are  shown.  The 
angle  of  attack  for  these  measurements  was  10  degrees.  All  three  sets  of 
data  agree  satisfactorily  with  the  pressure  distributions  calculated  from 
potential  flow  theory  except  in  the  region  near  the  suction  peak.  The 
divergence  between  theory  and  experiment  and  between  the  various  experi¬ 
mental  conditions  decreased  with  decreasing  angle  of  attack.  At  angles  of 
attack  higher  than  10  degrees,  the  divergence  between  theory  and  experiment 
increased  but  differences  between  the  various  experimental  configurations 
remained  about  the  same.  The  maximum  sectional  lift  coefficient,  obtained 
by  integration  of  the  pressure  coefficient  distributions,  occurred  at  an  angle 
of  attack  of  13  degrees  for  all  cases  and  the  numerical  values  ranged  between 
1. 24  and  1. 29. 

In  general,  the  uniform  flow  experiments  showed  that  the  addition 
of  turbulence  by  means  of  screens  influenced  the  data  only  very  slightly 
(see  Reference  9  for  further  details).  However,  since  there  was  a  small 
effect  of  turbulence  on  the  airfoil  pressure  distributions  in  uniform  flow,  it 
was  concluded  that  the  uniform  flow  data  used  for  comparison  with  the 
sheared-flow  data  should  be  those  obtained  with  the  low-solidity  turbulence 
screen  in  place.  These  data  were  chosen  over  those  with  the  high  solidity 
turbulence  screen  in  place  because  the  velocity  and  angularity  calibrations 
in  the  wind-tunnel  test  section  showed  that  a  more  nearly  uniform  flow  was 
obtained  with  the  low  solidity  screen.  The  free-stream  velocity  used  in  the 
uniform  flow  control  experiments  was  66  fps  which  is  approximately  midway 
between  the  minimum  and  maximum  velocities  encountered  in  the  sheared- 
flow  tests. 

Sheared-F’ow  Tests 

A  major  problem  in  discussing  the  aerodynamics  of  airfoils  in 
sheared  flows  is  the  question  of  what  constitutes  a  proper  reference  velocity. 


From  a  practical  viewpoint,  the  lift  generated  is  of  most  interest.  Hence, 
some  average  velocity,  taken  a3  a  constant  independent  of  airfoil  location, 
should  be  used.  Another  reference  velocity  which  has  been  used  is  the  velocity 
which  would  occur  in  the  undisturbed  sheared  flow  (with  no  airfoil  in  the  flow) 
at  the  location  of  the  midchord  of  the  airfoil.  This  velocity,  which  will  be 
called  the  midchord  velocity,  is  convenient  for  data  reduction  since  it  takes 
approximate  account  of  variations  in  velocity  with  variations  in  airfoil  loca¬ 
tion  in  the  sheared  flow.  Finally,  a  velocity  can  be  defined  which  takes  into 
complete  account  all  variations  in  the  effective  velocity  of  oncoming  flow, 
either  due  to  changes  in  airfoil  location  in  the  sheared  flow,  or  duo  to 
changes  in  airfoil  angle  of  attack  in  the  sheared  flow.  This  velocity  is  the 
one  which  occurs  in  the  undisturbed  sheared  Cow  far  upstream  on  the 
stagnation  streamline  (the  streamline  which  encloses  the  airfoil).  It  is  this 
multiplicity  of  reference  velocities  in  sheared  Cows  which  tends  to  cause 
confusion  and  may  encourage  misleading  conclusions  regarding  experimental 
aerodynamic  data  obtained  in  sheared  flows. 

The  dynamic  pressure  which  corresponds  to  the  stagnation  stream¬ 
line  reference  velocity  is  equal  to  the  difference  between  the  pressure  at  the 
stagnation  point  on  the  airfoil  and  the  free-stream  static  pressure.  If 
pressure  taps  are  placed  sufficiently  close  together  on  the  surface  of  the 
airfoil,  it  is  possible  to  obtain  this  dynamic  pressure  from  the  experimental 
airfoil  pressure  distributions.  Bernoulli's  equation  along  the  stagnation 
streamline  is 

P  +  J  p(**+  v2)  =  p0  i-  j p(U5)z 

where  pg  is  the  free-stream  static  pressure  and  U5  is  the  velocity  in  the 
undisturbed  sheared  flow  far  upstream  which  corresponds  to  the  stagnation 
streamline.  Using  this  relation,  the  pressure  coefficient,  (Cp)s  >  based  on 
the  dynamic  pressure  of  the  stagnation  streamline  becomes 


Note  that,  since  the  maximum  positive  value  of  (Cf,)s  is  unity,  the  appropriate 
dynamic  pressure  is  simply  the  maximum  positive  value  of  (  p- p0  )  which  is 
reached  on  the  airfoil.  In  addition,  the  remaining  two  reference  velocities 
which  were  mentioned  above,  provide  pressure  coefficients  Cp  (based  on  an 
average  velocity,  U  )  and  (Cp)R  (based  on  midchord  velocity,  0 ^  )  which  are 
defined  by 


^  ,  _  P  '  fa.  _  /J^s\Z_ 

=  ±P(DR)Z  "  (  uj  '  (Ut)* 


•Section  force  coefficients  can  be  referred  to  any  of  the  reference  velocities 
to  obtain,  for  example,  ,  (Cj  ) p  or  (Cj  )s  . 

The  pressure  distributions  obtained  on  the  Joukowski  airfoil  over 
a  wide  range  of  angles  of  attack  for  various  locations  in  the  sheared  flow 
were  reduced  to  the  coefficient  form  (Cp  )p  .  These  distributions  were  then 
integrated  numerically  to  obtain  the  section  lift  coefficient  (Cj  based 
on  midchord  velocity.  In  addition,  the  section  lift  coefficient,  ,  based 
on  the  average  velocity  of  the  undisturbed  sheared  flow,  U  ,  was  computed 
by  using  the  relation  Cj  -  \cj)  •  The  results  obtained  for  are 

shown  in  Figure  6  and  for  ( C £  )#  in  Figure  7,  both  plotted  versus  angle  of 
attack.  In  both  figures,  the  vertical  location  of  the  midchord  of  the  airfoil 
is  designated  by  the  value  of  h /c  in  the  legend,  where  h  is  the  height  of 
the  midchord  above  the  shcared-flow  centerline  and  c  is  the  chord  of  the 
airfoil. 

It  is  obvious  from  the  curves  of  versus  Ct  (Figure  6)  that  the 
magnitude  of  the  section  lift  is  strongly  influenced  by  the  location  of  the  air¬ 
foil  in  the  flow.  It  is  equally  obvious  from  the  curves  of  (C^  )  %  versus 
(Figure  7)  that  a  simple  linear  correction  for  the  variation  of  local  (midchord) 
dync.mic  pressure  with  airfoil  position  in  the  sheared  flow  does  not  completely 
account  for  the  large  variations  in  sectional  lift.  The  trends  shown  on  both 
of  these  figures  are  similar,  on  a  much  reduced  scale,  to  those  found  in 
Reference  6.  In  Reference  6,  however,  the  shear  in  the  flow  was  greater  by 
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a  factor  of  4  and  the  velocity  profile  contained  free  boundaries  as  in 
Figure  1(c).  Either  of  these  factors,  the  higher  shear  or  the  free  boundaries, 
might  account  for  the  larger  effects  of  airfoil  location  on  the  section  lift 
coefficient  curves  found  in  that  work. 

The  curves  of  Figures  6  and  7  show  that  at  zero  angle  of  attack 
there  is  a  positive  section  lift  coefficient  for  h  >  O  and  a  negative  section 
lift  coefficient  for  <  O  .  This  phenomenon  is  similar  to  that  predicted 
by  Tsien's  theory  (Reference  1)  for  a  symmetric  Joukowski  airfoil  in  a 
uniformly  sheared  flow.  Tsien's  theory  has  been  applied  to  the  current 
test  configuration  by  assuming  that  the  local  value  of  shear  at  the  airfoil 
midchord  could  be  used  to  predict  the  angle  of  attack  for  zero  lift  (  0cy:O  ). 

For  this  purpose,  the  experimental  velocity  profile  of  Figure  4  was  used  to 
estimate  the  mean  shear  at  the  airfoil  midchord  position.  Comparison 
between  values  of  Ocy.Q  calculated  in  this  way  and  experimental  values  of 
CK/.0  provided  no  conclusion  on  the  validity  of  the  assumption  that  local 
shear  may  be  used  to  compute  oc.j,0  in  a  nonuniformly  sheared  flow, 
primarily  because  of  the  uncertainties  in  finding  the  value  of  local  shear  from 
the  experimental  velocity  profiles. 

In  addition  to  the  change  in  angle  of  attack  for  zero  lift  with  airfoil 
location  in  the  nonuniformly  sheared  flow,  the  behavior  of  (Cy  )%  near  maximum 
lift  appears  to  be  dependent  on  the  magnitude  of  the  shear  (Figure  7).  At 
airfoil  positions  well  below  the  centerline  of  the  velocity  distribution,  the  stall 
is  abrupt.  At  positions  near  the  centerline  and  above  the  centerline  the 
stall  is  much  more  gradual.  Some  cases  which  exhibit  abrupt  stall  (  —  =  -1/6 
and  -1/3  for  example)  provide  higher  values  of  ( Cy  )#  just  before  stall  than 
those  obtained  when  the  airfoil  stalls  gradually  (  h.  -  +1/2,  +1/6,  and  +1/12). 

Some  degree  of  correlation  can  be  brought  to  the  ty  data  by 
referring  the  section  lift  to  the  stagnation  streamline  reference  velocity  as 
discussed  previously.  The  results  of  this  method  of  data  reduction  are 
shown  in  Figure  8  as  (  Cy  )$  versus  oc  =  Oe.  -Qy.0-  The  angle  of  attack  for 
zero  lift  has  been  subtracted  from  the  usual  angle  of  attack  in  this  figure 
to  facilitate  direct  comparison  of  the  shapes  of  the  curves  formed  by  the 
( Cy  )s  data.  Note  that  the  data  shown  in  Figure  8  lie  on  or  below  the 
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inviscid  uniform  flow  theoretical  curve  and  fall  into  two  approximate 
groupings;  one  for  positions  below  the  sheared  flow  centerline,  and  one  for 
positions  near  and  above  the  sheared  flow  centerline.  The  data  group 
corresponding  to  the  airfoil  near  and  above  the  sheared-flow  centerline 
exhibit  a  gradual  stall  and,  generally,  higher  values  of  (Cj  )s  at  all  positive 
values  of  &  than  the  data  group  corresponding  to  the  airfoil  below  the 
centerline.  This  behavior,  high  values  of  (Cj  )s  in  combination  with  gradual 
stall  and  low  values  of  [Cj  )5  in  combination  with  abrupt  stall,  is  in  contrast 
to  that  exhibited  by  the  ( C i  )f  data,  for  which  no  correlation  of  stall 
characteristics  with  magnitude  of  ( C 4  )*  just  prior  to  stall  was  found. 

Examination  of  the  pressure  distributions  measured  on  the  airfoil 
at  the  various  locations  in  the  sheared  flow  provides  at  least  a  partial  ex¬ 
planation  for  the  apparent  correlation  of  the  data  when  reduced  to  the 
form  (Cj)s  .  The  pressure  distributions  measured  on  the  airfoil  located  at 
=  +1/6  and  -1/6  and  at  an  angle  of  attack  of  12.9  degrees  for  each  case 
are  shown  in  Figure  9.  These  particular  pressure  distributions  are  presented 
because  they  resulted  in  approximately  the  iargest  differences  in  (  Cj  )*  data 
just  below  stall.  A  major  point  of  interest  in  the  data  of  Figure  9  is  that  the 
largest  contribution  to  the  difference  in  (  Cj  )*  occurs  on  the  airfoil's  lower 
surface,  with  the  data  for  —  =  -1/6  providing  the  highest  contribution 

to  (Cj,  )*  .  Apparently,  reducing  the  lift  data  to  the  form  (  Cj  )s  by  accounting 
for  the  differences  in  ( Cp  )#  at  the  stagnation  point  on  the  airfoil,  partially 
corrects  the  data  for  the  differences  in  the  lower  surface  pressure  distribu¬ 
tion.  However,  viscous  effects,  which  are  of  more  importance  on  the  upper 
surface,  and  the  differences  due  to  shear  in  angle  of  attack  for  zero  lift  are 
not  properly  accounted  for  by  reducing  the  data  to  the  form  (  Cj  }s  . 

The  upper  surface  pressure  distribution  shown  in  Figure  9  for 
h.  =  +1/6  has  a  substantially  lower  suction  peak  and  less  adverse  pressure 
gradient  just  downstream  of  the  peak  than  that  shown  for  h  =  -1/6.  These 
results  are  typical  of  all  the  pressure  distributions  measured  below  stall  but 
at  relatively  high  angles  of  attack,  except  for  some  cases  where  a  leading- 
edge  bubble  formed,  reducing  the  magnitude  of  the  suction  peak  which  was 
attained.  The  existence  of  such  bubbles  was  detected  through  lampblack  and 
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kerosene  flow  visualization  studies.  The  presence  was  also  detectable  in 
some  cases  by  irregularities  in  the  pressure  distribution  just  downstream 
of  the  suction  peak. 

It  was  found,  in  general,  that  those  configurations  which  exhibited  an 
abrupt  stall  were  also  configurations  which  resulted  in  the  sharpest  suction 
peaks  and  most  adverse  pressure  gradients  just  downstream  of  the  peak 
(as  for  Jl  =  -1/6  in  Figure  9)  up  to  angles  of  attack  where  the  formation 
of  a  leading-edge  bubble  did  not  severely  distort  the  upper  surface  pressure 
distribution.  The  abrupt  stall  which  was  obtained  resulted  from  the  bursting 
of  the  leading- edge  bubble.  In  contrast,  the  configurations  which  exhibited 
gradual  stall  were  configurations  in  which  bubble  formation  was  delayed  to 
much  higher  angles  of  attack  and  which  resulted  in  the  least  sharp  suction 
peaks  and  least  adverse  pressure  gradients  before  the  formation  of  the  bubble 
(as  for  h.  -  +1/6  in  Figure  9).  In  these  latter  cases,  stall  occurred  by  a 
progressive  forward  movement  of  the  upper  surface  separation  point. 

Both  of  the  above  cases  appear  to  be  consistent  with  the  hypothesis 
that  it  is  the  actioii  of  the  pressure  distribution  on  the  upper  surface 
boundary  layer  which  det>  rmines  the  stalling  characteristics  of  an  airfoil 
in  a  sheared  flow.  It  is  probably  not  necessary  to  consider  the  direct  inter¬ 
action  between  shear  in  the  flow  external  to  the  boundary  layer  and  the 
behavior  of  the  boundary  layer  when  investigating  the  stall  behavior  of  an 
airfoil  in  a  sheared  flow.  It  is,  however,  necessary  to  account  for  the 
manner  in  which  the  shear  affects  the  pressure  distribution  on  the  airfoil. 

In  summary,  the  above  observations  suggest  that  the  unusual  lift 
behavior  of  an  airfoil  in  a  two-dimensional  nonumiformly  sheared  flow  is  a 
consequence  of  at  least  two  more  or  less  independent  effects.  These  are: 
changes  in  stall  behavior  as  a  result  of  differences  in  the  upper  surface 
pressure  gradients,  and  changes  in  lift  behavior  associated  with  changes  in 
stagnation  pressure  of  the  streamline  which  intersects  the  airfoil  as  angle 
of  attack  is  varied.  It  is  worth  noting  that  the  stagnation  streamline  dynamic 
pressure  differed  from  the  dynamic  pressure  referenced  to  the  airfoil 
midchord  by  more  than  50  percent  in  some  instances. 
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COMPUTED  RESULTS  AND  COMPARISON  WITH  EXPERIMENT 


The  computer  program  has  been  used  to  compute  pressure  distribu¬ 
tions  on  the  Joukowski  airfoil  at  two  vertical  locations  in  the  nonuniformly 
sheared  flow;  ~  -  -1/3  and  —  =  +1/2.  At  —  =  -1/3,  pressure  distribu- 
tions  were  computed  for  (X  -  7 .  1  degrees  and  Ct  -  9.  5  degrees.  At  p-  =  +1/2, 
the  pressure  distribution  at  Ct  =8.2  degrees  was  computed. 

Figures  10  and  11  show  comparisons  between  computed  and 
experimental  values  of  ( Cp  )g  versus  for  =  -1/3  and  Cc  =  7.1  and 
Ct  -  9.  5  degrees.  Also  shown  on  these  figures  are  the  computed  pressure 
distributions  for  the  airfoil  in  uniform  flow  at  the  same  location  in  the  wind- 
tunnel  test  section.  These  latter  pressure  distributions  differ  from  the 
potential  flow  solutions  for  an  airfoil  in  an  inviscid  uniform  stream  of 
infinite  extent  because  of  the  presence  of  the  wind-tunnel  walls.  For  the 
particular  cases  computed,  however,  the  differences  were  very  small.  In 
both  figures,  the  computed  sheared  flow  pressure  distributions  agree  satis¬ 
factorily  with  the  experimental  data  on  the  lower  surface  and  both  experi¬ 
mental  and  computed  results  differ  appreciably  from  the  uniform  flow 
pressure  distribution  in  the  vicinity  of  the  stagnation  point.  The  moderate 
differences  between  theory  and  experiment  near  the  midchord  on  the  lower 
surface  are  about  the  same  as  those  found  in  uniform  flow,  an  example  of 
which  is  given  in  Figure  5.  At  Ct  =  7.  1  degrees,  the  computed  upper  surface 
pressure  distribution  agrees  well  with  the  experimental  data  except  near  the 
suction  peak.  At  the  higher  angle  of  attack,  Ct  -  9.5  degrees,  the  experi¬ 
mental  data  differ  from  the  computed  pressure  distribution  over  the  majority 
of  the  upper  surface.  In  both  cases,  the  difference  between  the  computed 
pressure  distributions  and  the  experimental  data  for  the  sheared  flow  were 
comparable  to  the  differences  observed  between  theory  and  experiment  in 
the  uniform  flow  control  tests  at  comparable  angles  of  attack.  These 
differences  are,  therefore,  probably  attributable  to  viscous  effects  which 
become  larger  as  angle  of  attack  is  increased.  Note  that  at  both  angles  of 
attack,  the  computed  pressure  gradient  immediately  downstream  of  the 
suction  peak  is  more  adverse  than  that  computed  for  uniform  flow.  This 
result  is  consistent  with  the  observed  abrupt  stall  which  occurs  with  the 

airfoil  at  fl  =  —  1  / 3 . 
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Figure  12  presents  a  comparison  of  computed  and  experimental 
pressure  distributions  for  h.  -  +1/2  and  (X  -  8.2  degrees.  The  agreement 
between  theory  and  experiment  for  this  case  is  very  good  over  most  of  the 
airfoil  surface.  The  experimental  data  show  a  small  loss  in  the  suction  near 
the  peak  and  a  small  increase  in  stagnation  pressure.  The  reason  for  the 
latter  difference  is  not  known  at  present  but  may  be  attributable  to  experi¬ 
mental  error.  The  loss  near  the  suction  peak  in  this  case,  again  probably 
due  to  viscous  effects,  is  quite  small  because  of  the  mild  pressure  gradient 
on  the  upper  surface.  Note  that  the  computed  and  experimental  pressure 
distributions  lie  above  the  computed  uniform  flow  curve.  Such  behavior 
cannot  be  attributed  to  viscous  effects.  These  would  reduce  the  experimental 
suction  pressures  on  the  upper  surface.  The  computed  pressure  gradient 
immediately  downstream  from  the  suction  peak  is  less  adverse  for  the 
sheared  flow  than  for  the  uniform  flow  in  this  case.  This  is  consistent 
with  the  gradual  stall  observed  with  the  airfoil  in  this  location. 

The  computed  pressure  distributions  were  integrated  to  get  lift, 
moment,  and  pressure  drag  coefficients  referenced  to  midchord  dynamic 
pressure.  The  results  are  compared  to  appropriate  experimental  results 
in  Reference  9.  Agreement  between  computed  and  experimental  values  of 
(Cj  anc*  ^Cmc/^R  was  satisfactory,  except,  perhaps,  for  (  Cmc/,  )r  at 
_  =  -1/3,  Oc  =  9.  5  degrees  where  the  computed  value  was  -0.  01  1  and  the 

experimental  value  -0.028.  The  pressure  drag  coefficients  computed  from 
the  theory  for  the  sheared  flow  and  the  uniform  flow  c  onfigurations  varied 
between  -0.  01  and  -0.  016.  The  negative  values  may  reflec  ;  r 

prediction  of  leading-edge  suction  peak  for  uniform  flow  which  undoubtedly 
carries  over  into  the  nonuniformly  sheared  flow  calculations.  Jn  contrast, 
the  experimental  pressure  drag  coefficients  had  minimum  values  of 
approximately  0.005. 

CONCLUDING  REMARKS 

The  good  agreement  between  the  computed  and  experimental 
pressure  distributions  (allowing  for  viscous  effects  of  similar  magnitude  to 
those  ooserved  in  uniform  flow)  and  the  fact  that  the  computed  distributions 
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differ  in  important  respects  from  computed  uniform  flow  pressure  distribu¬ 
tions  indicate  that  the  theoretical  model  successfully  predicts  the  important 
features  of  the  airfoil  aerodynamic  characteristics  in  the  unseparated  non- 
uniformly  sheared  flow.  A  particular  notable  feature  of  the  computed  results 
is  that  they  appear  to  account  for  the  stall  behavior  solely  on  the  basis  of  the 
inviscid  flow  pressure  distribution. 

For  the  particular  airfoil  and  nonuniformly  sheared  flow  used  in  the 
experimental  portion  of  this  program,  a  marked  degree  of  correlation  of  the 
section  lift  data  is  obtained  for  differing  vertical  positions  in  the  flow  if 
section  lift  coefficients  are  referenced  to  the  airfoil  stagnation  streamline 
dynamic  pressure  (  ),  and  angle  of  attack  is  measured  from  the  zero- 

lift  angle  of  attack.  The  two  primary  factors  in  determining  the  observed 
two-dimensional  aerodynamic  behavior  of  airfoils  in  nonuniformly  sheared 
flows  then  appear  to  be  the  variation  of  stagnation  streamline  dynamic 
pressure  with  angle  of  attack  for  a  given  wing  position  in  the  sheared  flow, 
and  the  inviscid  influence  of  the  nonuniformly  sheared  flow  in  determining  the 
pressure  distribution  on  the  airfoil  upper  surface.  The  former  factor  deter¬ 
mines  the  apparent  lift-curve  slope  well  below  stall  and  the  latter  determines 
the  stalling  behavior.  Both  factors  appear  to  depend  on  the  overall  distribu¬ 
tion  of  shear  in  the  flow. 

The  experimental  results  presented  in  this  paper  and  in  Reference  6 
have  indicated  that  airfoil  location  in  a  nonuniformly  sheared  flow  plays  a 
substantial  role  in  determining  the  lift  and  stalling  characteristics  of  the 
airfoil.  The  results  presented  in  Reference  6,  which  were  measured  with 
the  same  airfoil  in  a  two-dimensional  nonuniformly  sheared  flow  more 
representative  of  a  section  through  a  propeller  slipstream,  have  shown  that 
the  differences  in  lift  characteristics  can  become  very  large  indeed.  For 
example,  relocation  of  the  airfoil  midchord  in  a  region  of  large  shear 
gradient  from  a  position  approximately  one-half  of  its  thickness  above  the 
sheared  flow  centerline  to  a  position  equally  distant  from  but  below  the 
centerline  increased  the  maximum  sectional  lift  available  before  stall  by  a 
factor  of  more  than  2.  Such  large  differences  in  the  sectional  lift  available 
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before  stall  plus  the  change  in  stalling  characteristics  in  the  two-dimensional 
nonuniformly  sheared  flow  suggest  that  there  may  be  appreciable  benefits 
available  for  a  wing  immersed  in  a  three-dimensional  propeller  slipstream 
with  significant  shear  if  the  overall  configuration  is  properly  designed. 

The  theoretical  model  outlined  briefly  in  this  paper  would,  with 
extensions  to  account  for  additional  complexities  of  the  flow  such  as  free 
boundaries,  provide  a  powerful  tool  for  investigation  of  the  effec  of  non- 
uniform  shear  on  aerodynamic  characteristics  of  arbitrary  airfoils  in  two- 
dimensional  nonuniformly  sheared  flows.  The  results  of  these  investigations 
could  be  used  as  a  starting  point  to  attack  the  more  complex  three- 
dimensional  problem  of  a  wing  immersed  in  a  propeller  slipstream  with 
significant  nonuniform  shear. 
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igure  2.  SCHEMATIC  OF  FLOW  MODEL  FOR  IBM  COMPUTER  PROGRAM. 


AIRFOIL  COORDINATE,  x/c  ,  FRACTION  OF  CHORD 


PRESSURE  COEFFICIENT,^,  VERSUS  FRACTION  OF  CHORD, X/t, 

FOR  UNIFORM  FLOW  FROM  POTENTIAL  FLOW  THEORY,  AND  FROM 
COMPUTER  PROGRAM  WITH  AND  WITHOUT  SMOOTH  I NG ;«C= 1 0 . 0  DEGREES. 
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Figure  4.  MEASURED  FREE-STREAM  VELOCITY  DISTRIBUTION  IN  NONUNIFORMLY 
SHEARED  FLOW  AT  MODEL  MIDCHORD  STATION. 
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PRESSURE  COEFFICIENT 


SECTION  LIFT  COEFFICIENT  REFERENCED  TO  DYNAMIC  PRESSURE  CORRESPONDING  TO  AVERAGE  V-LOCIT''. 


ANGLE  OF  ATTACK,  or, degrees 


Figure  6.  NOND I  MENS  I ONAL  SECTION  LIFT,?,  ,  VERSUS  ANGLE  OF  ATTACK,  cc, 
NONUNI FORMLY  SHEARED  FLOW. 
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SECTION  LIFT  COEFFICIENT  REFERENCED  TO  MIDCHORD  FREE-STREAM  DYNAMIC  PRESSURE,  (c 


St  T ION  LIFT  COEFFICIENT  REFERENCED  TO  DYNAMIC  PRESSURE  OF  STAGNATION  STREAMLINE,  (Cj) 


Figure  8.  SECTION  LIFT  COEFFICIENT  REFERENCED  TO  DYNAMIC  PRESSURE  OF 
STAGNATION  STREAMLINE,Cct)s,  VERSUS  ANGLE  OF  ATTACK  MEASURED 
FROM  ZERO  LIFT, a  ,  NONUNIFORMLY  SHEARED  FLOW. 
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PRESSURE  COEFFICIENT  REFERENCED  TO 
MIDCHORD  DYNAMIC  PRESSURE, (Cp)* 


PRESSURE  COEFFICIENT  REFERENCED  TO 
MIOCHORD  DYNAMIC  PRESSURE,  (C p)/z 


AIRFOIL  COORDINATE,  t  c  .  FRACTION  OF  CHORD 


Figure  10.  COMPARISON  OF  COMPUTED  AND  MEASURED  PRESSURE  COEFFICIENTS, 
Up)g,  VERSUS  FRACTION  OF  CHORD, i/c,  NONUNIFORMLY 
SHEARED  FLOW,  h/c  =-i/3,cx  =  7.  I  DEGREES. 


PRESSURE  COEFFICIENT  REFERENCED  TO 
MIDCHORD  DYNAMIC  PRESSURE, CCp)g 


PRESSURE  COEFFICIENT  REFERENCED  TO 
MIDCHORD  DYNAMIC  PRESSURE, 


Figure  12.  COMPARISON  OF  COMPUTED  AND  MEASURED  PRESSURE  COEFFICIENTS, 
VERSUS  FRACTION  OF  CHORD, x/c,  NONUNIFORMLY 
SHEARED  FLOW,h/c  =  + 1/2  ,<X.  =8. 2  DEGREES. 


35 


NOTES 


EXPERIMENTAL  INVESTIGATION  OF  COMPOUND 
HELICOPTER  AERODYNAMIC  INTERFERENCE  EFFECTS 

by 

RICHARD  M.  SEGEL 

Supervisor,  Advanced  Rotorcraft  Research  Projects  Group 
Aircraft  Advanced  Research  Section 

LAWRENCE  J.  BAIN 
Supervisor,  Dynamics  Research  Group 
Aircraft  Advanced  Research  Section 

(6)  Sikorsky  Aircraft,  Division  of  United  Aircraft  Corporation 
Stratford,  Connecticut,  U.  S.A.,  1966.  All  Rights  Reserved. 


ABSTRACT 

The  presence  of  numerous  lifting  components  on  a  compound 
helicopter  results  in  mutual  aerodynamic  interferences  between  them. 

A  generalized  model  of  a  compound  helicopter  was  tested  during  March, 

1966  in  the  United  Aircraft  18  Foot  Wind  Tunnel  to  investigate  these 
effects  at  speeds  up  to  300  knots.  Three  strain  gage  balances  measured 
rotor,  wing  and  fuselage  forces  and  moments  independently.  In  addition, 
wing  static  pressure  distribution,  oscillatory  airloads,  tail  downwash 
angle,  and  rotor  blade  flapping  ard  stresses  were  measured.  Configuration 
variables  included  three  wing  sizes,  each  tested  at  three  vertical  locations 
on  the  fuselage. 

A  preliminary  analysis  of  :he  data  is  presented  which  shows  the 
decrease  in  wing  loading  induced  by  rotor  downwash,  and  the  influence 
of  wing  size  and  position  on  this  effect.  The  wing  lift/drag  ratio  is  shown 
to  decrease  with  increasing  rotor  lift.  In  the  practical  operating  range  of 
compound  helicopters,  the  influence  of  the  wing  upon  rotor  characteristics 
was  found  to  be  small.  The  measured  downwash  at  the  tail  conforms  to 
fully  developed  momentum  values. 


INTRODUCTION 


The  successful  marriage  of  fixed  wing  and  helicopter  technologies 
in  the  development  of  the  compound  helicopter  requires  an  understanding  of 
the  aerodynamic  interactions  among  the  vehicle's  lifting  components.  The 
relative  effects  of  a  number  of  design  and  operating  parameters  such  as 
wing  size  and  vertical  position,  and  rotor/wing  lift  sharing  upon  performance 
and  handling  characteristics  must  be  determined  before  the  concept  can  be 
fully  developed.  This  paper  describes  an  AVLABS  sponsored  experimental 
investigation  of  the  interference  effects  occurring  in  compound  helicopters 
between  rotor,  wing,  fuselage  and  horizontal  stabilizer  over  a  broad  range 
of  configuration  variables,  Reference  1. 

The  existence  of  these  interference  effects  was  revealed  in  small 
scale  wind  tunnel  tests  which  were  conducted  as  early  as  1955  and  more 
recently  in  the  UAC  4x6  foot  wind  tunnel.  In  addition,  flight  research  in 
this  area  began  with  the  addition  of  wings  to  an  S-55  helicopter  in  1957-8  and 
is  continuing  at  the  present  with  the  Sikorsky  S-61F  research  aircraft, 
shown  in  Figure  1.  This  aircraft  is  instrumented  to  measure  rotor  and 
wing  loads  independently  and  has  attained  speeds  in  excess  of  200  knots. 

The  S-61F  program  and  other  AVLABS  sponsored  flight  research  investi¬ 
gations  of  compound  helicopters  have  helped  to  define  the  effects  of  rotor/ 
wing  lift  sharing,  but  none  has  provided  the  flexibility  or  precision  which 
may  be  obtained  in  the  wind  tunnel.  Consequently,  AVLABS  awarded  a 
contract  to  Sikorsky  Aircraft  in  February  1965,  to  perform  a  comprehensive 
wind  tunnel  investigation  to  study  the  aerodynamic  characteristics  and 
interference  effects  of  rotor -wing -fuselage  combinations. 

WIND  TUNNEL  MODEL 

The  wind  tunnel  model  used  in  this  investigation  is  shpwn  in 
Figure  2.  Representing  a  generalized  compound  helicopter,  it  is  approx¬ 
imately  nine  feet  long  and  is  equipped  with  a  four  blade  dynamically  scaled 
rotor  nine  feet  in  diameter.  Rotor  power  is  supplied  by  a  19  horsepower 
electric  motor  through  a  reduction  gearbox.  The  rotor  controls  are 
hydraulically  actuated  and  are  capable  of  duplicating  all  motions  of  a  full 
scale  control  system.  The  model  is  equipped  with  three  internal  six- 
component  strain  gage  balances  which  measure  forces  and  moments  on  the 
rotor,  wing  and  fuselage  independently.  A  floating  horizontal  tail  surface 
allows  measurement  of  the  mean  tail  downwash  angle.  In  order  to  measure 
the  various  effects  of  wing  size  and  vertical  location,  the  model  was 
constructed  with  a  rectangular  fuselage  cross-section  in  the  vicinity  of  the 
wing  mounting  points,  thus  allowing  the  mounting  of  any  of  the  three  wings 
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at  any  vertical  location.  The  wings  are  mounted  from  side  plates  which  are 
attached  to  a  box  structure  suspended  from  the  wing  balance.  The  wing  loads 
are  thus  measured  independently  of  fuselage  loads.  The  three  wings  have 
similar  planforms  with  an  aspect  ratio  of  6  and  taper  ratio  of  1.5.  The 
wing  spans  are  equal  to  30,  50  and  75  percent  of  rotor  diameter  and  cover 
the  range  of  interest  for  most  compound  configurations.  Three  wing  vertical 
positions,  determined  by  model  geometry  are  10.9,  15.7  and  20.4  percent 
of  the  rotor  diameter  below  the  rotor.  Thus  nine  combinations  of  wing  size 
and  location,  in  addition  to  the  wingless  configuration  were  tested.  The 
general  arrangement  and  dimensions  of  the  model  are  shown  in  Figure  3. 

The  fiberglass  rotor  blades  were  dynamically  scaled  about  t-y 
three  axes  to  allow  simulation  of  forward  speeds  and  rotor  speeds  which  are 
twice  the  actual  values.  These  blades  were  fully  strain  gaged  to  measure 
the  affects  of  the  v/ing  upon  vibratory  stress.  The  "flight”  conditions  which 
were  investigated  are  shown  in  Figure  4,  a  plot  of  rotor  tip  speed  versus 
forward  speed.  The  five  test  speeds  simulate  operation  from  120  knots  to 
high  speed  compound  flight  at  300  knots. 

ANALYSIS  OF  TEST  DATA 

Effect  of  Rotor  On  Wings 

The  effect  of  rotor  downwash  upon  wing  loading  and  performance 
was  investigated  for  the  nine  configurations  throughout  the  speed  range.  As 
was  anticipated,  the  most  significant  measured  effect  was  a  decrease  in  wing 
lift  which  was  approximately  proportional  to  the  rotor  momentum  downwash. 
An  example  of  this  effect  is  shown  for  the  large  wing  at  120  knots  in  Figure 
5.  The  wing  lift  coefficient,  CL\y  is  plotted  as  a  function  of  rotor  lift 
coefficient -solidity  ratio,  Cl r/O-’  >  for  several  values  of  fuselage  angle 
of  attack,  a  f.  In  addition  to  the  linear  relationship  between  wing  and  rotor 
lift,  it  is  seen  that  rotor  trim  condition  affects  the  magnitude  of  the  lift 
decrease.  For  example,  at  Clr  /&  =  -04  and  Of  =  4  degrees,  an  increase 
in  collective  pitch  irom  0  degrees,  where  the  rotor  is  tilted  rearward  and 
near  autorotation;  to  8  degrees,  where  the  rotor  is  tilted  forward;  results 
in  an  increase  in  Cl^  of  .06.  This  is  believed  due  to  the  change  in  rotor 
loading  distribution  and  the  corresponding  change  in  orientation  between 
the  rotor  disk  and  wing. 

The  linear  decrease  in  wing  lift  may  be  expressed  as  a  change  in 
wing  effective  angle  of  attack.  This  effective  rotor  downwash  angle  is 
compared,  for  the  large  wing,  with  the  value  of  the  momentum  downwash 
in  the  rotor  disk  in  Figure  6.  The  average  experimental  rotor  downwash 
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angle  per  unit  Cjn  /  O'  is  presented  as  a  function  to  forward  velocity  for  the 
large  wing  at  low,  mid,  and  high  positions.  The  downwash  angle  for  the 
low  wing  is  approximately  equal  to  the  momentum  value  and  increases  as 
the  wing  is  placed  nearer  to  the  rotor.  Expressing  this  angle  as  a  linear 
function  of  the  momentum  downwash  gives  the  expression 

Aa“  =k^2^"  [Clr/ct]  <Eq.  D 

The  average  experimental  value  of  the  constant,  k,  was  determined 
for  the  medium  and  large  wings  and  is  shown  in  Figure  7  as  a  function  of 
the  nondimensional  vertical  separation  between  wing  and  rotor,  Z/D.  The 
effect  of  the  downwash  on  the  large  wing  is  significantly  larger  than  on  the 
medium  wing,  and  increases  for  both  wings  approximately  linearly  as  the 
wing/rotor  spacing  is  reduced  in  the  practical  range  of  interest  shown. 

Although  the  effect  of  rotor  downwash  on  the  total  wing  lift  was 
expressed  as  a  simple  angle  of  attack  change,  the  effect  on  wing  drag  was 
found  to  be  more  complex.  An  increase  in  rotor  downwash  caused  an 
increase  in  wing  drag  due  to  the  rearward  inclination  of  the  wing  resultant 
force  vector.  When  combined  with  the  loss  of  lift,  this  resulted  in  a 
signficiant  deterioration  of  the  wing  lift/drag  ratio.  An  example  of  this 
effect  is  shown  in  Figure  8  for  the  large  wing  at  120  knots  the  lowest  speed 
investigated.  To  obtain  this  plot,  wing  lift  and  drag  coefficients  were 
plotted  against  Clr  ,(j  »  and  the  results  were  crossplotted  as  shown. 

It  may  be  seen  that  the  decrease  in  L/D  is  due  primarily  to  a  large  decrease 
in  lift  at  nearly  constant  drag.  This  effect  decreases  rapidly  with  forward 
speed  as  the  rotor  downwash  angle  decreases.  Consequently  at  high  speeds, 
where  a  decrease  in  wing  L/D  would  be  critical,  the  effect  is  small. 

In  order  to  obtain  a  more  complete  picture  of  the  rotor  downwash 
effects,  the  medium  wing  was  constructed  with  static  pressure  taps  at  40 
and  70  percent  semi-span  on  top  and  bottom  surfaces.  The  resulting 
static  pressures  were  examined  and  the  chordwise  distributions  integrated 
to  obtain  wing  local  normal  force  coefficients. 

An  interesting  effect  was  obtained  at  120  knots  where  the  spanwise 
normal  force  distribution  is  shown  in  Figure  9  for  the  high  and  raid  positions 
and  several  values  of  rotor  lift.  In  the  upper  plot,  three  curves  are  shown 
for  the  high  wing  with  no  rotor,  and  for  values  of  C Lr/C  ^qual  to  .033 
and  .063.  The  effect  of  the  downwash  is  to  decrease  tne  lift  on  the  star¬ 
board  wing,  beneath  the  advancing  blade,  without  a  large  effect  on  the  port 
wing.  However,  with  the  wing  in  the  mid  position,  the  lift  is  more  sym¬ 
metrical  as  shown  in  the  lower  plot.  The  downwash  effect  upon  the  low 
wing  was  also  found  to  be  symmetrical. 
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Effect  of  Wings  on  Rotor 

A  major  portion  of  this  investigation  was  directed  toward 
establishing  the  effects  of  wing  geometry  and  loading  upon  rotor  perform¬ 
ance  and  stresses.  The  model  was  initially  tested  without  wings  through 
the  speed  range  and  over  a  wide  range  of  rotor  operating  conditions  to 
obtain  basic  rotor  characteristics.  Comparison  of  the  basic  data  with 
values  measured  in  the  presence  of  the  wings  revealed  that  the  effect  of 
any  of  the  wings  upon  the  rotor  is  small  at  speeds  below  200  knots. 

Figure  10  illustrates  the  effect  of  the  large  wing  upon  rotor  lift 
and  angle  of  attack  at  120  knots.  These  data  were  obtained  by  trimming 
out  first  harmonic  flapping  with  respect  to  the  shaft.  Thus  the  constant 
values  of  Of  are  equivalent  to  equal  tip -path  plane  angles  of  attack.  The 
wing  free  stream  angle  of  attack  is  equal  to  Of  +  4  degrees  because  of  the 
four  degrees  wing  incidence.  It  may  be  seen  that  the  wing,  even  at  high 
angles  of  attack,  causes  only  modest  decreases  in  rotor  lift.  Similar  data 
are  shown  for  the  large  wing  in  Figure  11  and  12  for  velocities  of  200  and 
300  knots  respectively.  At  200  knots,  the  decrease  in  Clr /cr  is  more 
pronounced,  while  at  300  knots,  the  interference  effect  is  substantial. 
However,  it  should  be  noted  that  the  large  wing  is  oversized  for  cruise  at 
300  knots.  The  condition  at  df  =  4  degrees  with  four  degrees  incidence  is 
equivalent  to  a  wing  lift  more  than  twice  the  hovering  lift  capacity  of  the 
rotor. 


The  effect  of  the  medium  size  wing  upon  the  rotor  is  much  smaller 
than  that  of  the  large  wing.  For  example  at  200  knots,  Figure  13,  there  is 
no  measurable  loss  of  rotor  lift  at  any  condition.  The  effect  of  the  medium 
and  large  wings  upon  rotor  lift  is  summarized  in  Figure  14,  a  plot  of  the 
change  in  Cl p /<t  Per  degree  change  in  wing  angle  of  attack  as  a  function 
of  forward  speed. 

It  is  seen  that  the  interference  effect  increases  rapidly  with 
forward  speed,  but  that  it  is  relatively  large  only  for  the  large  wing. 
Furthermore,  since  the  rotor  lift  sensitivity  with  respect  to  changes  in 
angle  of  attack  also  increases  rapidly  with  forward  speed,  the  wing  induced 
loss  in  rotor  lift  is  easily  offset  by  a  small  cyclic  pitch  change. 

Rotor  blade  vibratory  stresses  were  measured  at  five  flapwise, 
three  chordwise  and  two  torsion  stations  on  one  blade.  Time  histories 
and  peak-to-peak  values  were  recorded  at  all  test  conditions  and  analysis 
of  a  considerable  portion  of  these  data,  including  the  large  wing  at  the  high 
position,  has  not  revealed  any  significant  effect  of  the  wing  on  blade  stresses. 
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Tail  Downwash  Angle 

To  investigate  the  downwash  effects  on  a  tail,  a  0012  airfoil 
section  with  a  span  equal  to  36  percent  of  the  rotor  diameter  was  mounted 
as  a  floating  tail  on  the  compound  helicopter  model.  The  integrated  down- 
wash  angle  at  the  tail  was  measured  by  means  of  a  potentiometer  which 
sensed  the  tail  angle  of  attack.  This  angle  was  measured  over  simulated 
forward  speeds  from  120  to  300  knots,  fuselage  angles  of  attack,  df 
from  -8  to  +  12  degrees  and  rotor  lift  coefficient -solidity  ratios,  Clr/<t  , 
from  zero  to  0. 1. 

The  downwash  angle  resulting  from  the  wing  at  the  mid  position 
and  the  fuselage  is  shown  in  Figure  15  as  a  function  of  a  *  Both  theoretical 
*  (Reference  2)  and  experimental  values  are  included  ana  the  correlation  for 
all  three  size  wings  is  good.  The  measured  angle  for  the  fuselage  alone  is 
also  shown.  In  Figure  16  the  downwash  angle  resulting  from  the  rotor  is 
shown  for  forward  speeds  of  120,  200  and  300  knots  as  a  function  of  Clr/q-  • 
The  derivative  of  the  downwash  angle  with  respect  to  Clr  j&  varies 
inversely  with  the  forward  speed  as  would  be  expected  from  momentum 
considerations.  The  net  angle  predicted  by  the  fully  developed  momentum 
flow  can  be  checked,  for  example,  at  Clr/  a  -  0.04  and  200  knots.  The 
calculated  downwash  angle  is  2.6  degrees  which  compares  with  the  measured 
value  of  3.5  degrees  including  an  effect  of  approximately  1  degree  from  the 
fuselage.  Also  shown  in  Figure  16  is  the  variation  of  the  downwash  angle 
with  including  the  effects  of  the  small  and  medium  wings  at  the  mid 

position.  At  the  zero  rotor  lift  condition  the  angle  for  the  medium  wing  is 
on  the  order  of  4.5  degrees.  The  addition  of  the  angle  for  the  wing  alone 
from  Figure  15  and  the  angle  for  the  rotor  alone  also  gives  approximately 
4.5  degrees,  when  one  degree  is  deducted  to  avoid  double  inclusion  of  the 
fuselage  effect.  For  larger  values  of  Clr/ct  this  type  of  addition  is  no 
longer  possible  since  the  contribution  of  tne  wing  to  the  net  downwash  angle 
decreases  with  increasing  rotor  loading.  This  effect  is  believed  to  be 
caused  by  the  rotor  downwash  displacing  the  wing  trailing  vortices  further 
below  the  tail,  thus  reducing  their  influence  on  the  downwash  angle.  The 
curves  show  that  with  the  small  wing  the  downwash  angle  at  the  tail  is 
essentially  unaffected  by  the  wing  when  Cf- /&  reaches  a  value  of  0.06. 
Although  operation  of  the  rotor  much  beyoncr  this  point  at  200  knots  is  not 
generally  practical  it  can  be  seen  that  the  medium  wing  curve  is  also  tending 
toward  a  union  with  the  rotor  only  curve. 


*  Because  of  the  unusually  large  nondimensional  distance  of  the  tail  behind 
the  wing  it  was  necessary  to  extrapolate  the  published  data. 
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The  results  shown  in  Figures  15  and  16  are  representative  of  the 
remainder  of  the  tail  downwash  data,  but  may  not  be  applicable  to  wing-tail 
configurations  other  than  those  used  in  this  program. 

Fuselage  and  wing  balance  measurements  taken  over  a  range  of 
conditions  with  the  three  wings  positioned  at  the  high,  mid  and  low  points 
on  the  fuselage  were  reduced  and  compared  with  various  published  NACA 
test  data  on  fuselage -wing  combinations  and  no  significant  departure  from 
those  results  was  found.  The  complete  analysis  of  these  measurements 
will  be  given  in  the  final  program  report. 

ADDITIONAL  RESEARCH 

Currently  being  implemented  under  an  AVLABS  contract  is  a 
rotor -propeller  interference  investigation  through  the  forward  speed  range 
to  300  knots.  For  this  program  dynamically  scaled  propellers  will  be 
mounted  on  the  model  in  various  configurations  to  obtain  experimental 
stress  and  aerodynamic  data  necessary  for  correlation  with  theory  and  for 
the  design  of  high  performance  rotary-fixed  wing  aircraft  using  propellers 
for  auxiliary  forward  thrust. 

Future  studies  with  the  compound  helicopter  model  which  would 
provide  valuable  data  include: 

a.  A  repetition  of  a  selected  number  of  the  test  conditions 
reported  in  the  present  paper  utilizing  metal  dynamic  model  rotor  blade 
scaled  for  full  scale  Mach  number.  The  results  of  this  program  could  be 
compared  with  the  present  data  to  define  the  effects  of  compressibility 
on  rotor -wing  interference. 

b.  An  exploratory  investigation  of  rotor  instabilities  including 
stall  and  classical  flutter  and  torsional  divergence.  In  addition  transient 
rotor  response  due  to  step  control  inputs  and  angle  of  attack  change  could 
be  studied. 

CONCLUSIONS 

Based  on  the  preliminary  analysis  of  data,  the  following  conclusions 
were  reached: 

1.  The  rotor  induced  downwash  at  the  wing  is  approximately 
equal  to  the  momentum  downwash  in  the  rotor  disk.  However  wing  size 
and  vertical  distance  from  the  rotor  are  important  parameters  in  accurately 
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determining  the  actual  mean  downwash  value. 

2.  In  addition  to  the  loss  in  wing  lift  caused  by  the  rotor  down- 
wash,  a  small  increase  in  wing  drag  occurs.  The  combined  result  is  a 
decrease  in  the  wing  lift /drag  ratio,  at  low  speeds.  However  this  effect 
decreases  with  forward  speed  and  is  small  when  wing  performance  becomes 
important. 

3.  At  pure  helicopter  speeds,  the  rotor  downwash  induces  an 
unsymmetrical  spanwise  decrease  in  wing  loading.  The  larger  reduction  in 
wing  lift  occurs  beneath  the  advancing  blade.  This  effect  is  largest  for  the 
high  wing  configuration. 

4.  The  effec'  of  the  wing  is  to  decrease  rotor  lift.  However 
the  loss  in  lift  is  small  ai  low  speeds  for  all  conditions.  It  becomes 
large  only  at  high  speeds  with  large  wings  at  loading  corresponding  to 
extreme  operating  conditions  and  is  not  considered  to  be  significant  in  the 
design  of  a  compound  helicopter. 

5.  The  presence  of  a  wing  has  no  significant  effect  on  rotor 
blade  stresses. 

6.  The  rotor  induced  downwash  at  rhe  tail  is  approximately 
equal  to  the  full  developed  momentum  downwash  in  the  absence  of  the  wing. 
The  contribution  to  the  downwash  at  the  tail  induced  by  the  wing,  in  the 
configurations  tested,  decreases  rapidly  as  rotor  lift  is  increased. 
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FIGURE  I.  S-6IF  RESEARCH  AIRCRAFT. 


FIGURE  2.  SIKORSKY  COMPOUND  HELICOPTER  MODEL  IN  18  FOOT  WIND  TUNNEL 


FIGURE  3.  GENERAL  ARRANGEMENT  OF  SIKORSKY 
COMPOUND  HELICOPTER  MODEL. 
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FIGURE  4.  ROTOR  TEST  CONDITIONS, 
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ROTOR  LIFT  COEFFICIENT- SOLIDITY  RATIO  ,  CLr /«■ 

FIGURE  5.  EFFECT  OF  ROTOR  LIFT 


ADVANCE  RATIO,  fi 

FIGURE  6.  EFFECT  OF  FORWARD  SPEED  ON  EFFECTIVE 
ROTOR  OOWNWASH  ANGLE. 


WING  DRAG  COEFFICIENT, 


FIGURE  7  EFFECT  OF  WING  POSITION 
ON  DOWNWASH  FACTOR. 


FIGURE  8.  EFFECT  OF  ROTOR  LIFT  ON  WING  CL  AND  Cq. 


14 


LOCAL  NORMAL  FORCE  COEFFIENT 


V  ■  120  KNOTS 


PERCENT  SEMISPAN 

FIGURE  9.  SPANWISE  NORMAL  FORCE  DISTRIBUTION 
ON  THE  MEOIUM  WING  AT  120  KNOTS. 
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FIGURE  10.  ROTOR  LIFT  COEFFICIENT  -  SOLIOITY  RATIO 
V.S.  ROTOR  ANGLE  OF  ATTACK  AT  120 KNOTS. 


ROTOR  LIFT  COEFFICIENT -SOLIDITY  RATIO, 


ROTOR  ANGLE  OF  ATTACK  ,  «c  ,  DEG. 


FIGURE  II.  ROTOR  LIFT  COEFFICIENT- SOLIDITY  RATIO 
V.S.  ROTOR  ANGLE  OF  ATTACK  AT  200  KNOTS. 


LIFT  COEFFICIENT -SOLIDITY  RATIO,  CLR/, 


FIGURE  12.  ROTOR  LIFT  COEFFICIENT- SOLIDITY  RATIO 
VS.  ROTOR  ANGLE  OF  ATTACK  AT  300  KNOTS. 


ROTOR  LIFT 


ROTOR  ANGLE  OF  ATTACK, «e  , DEG. 


FIGURE  13.  ROTOR  LIFT  COEFFICIENT- SOLIDITY  RATIO 
V.S.  ROTOR  ANGLE  OF  ATTACK  AT  200  KNOTS. 
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FIGURE  14.  EFFECT  OF  FORWARD  SPEEO  ON  CHANGE  IN 
ROTOR  LIFT  PER  UNIT  CHANGE  IN  WING  ANGLE  OF  ATTACK. 


20 


DOWNWASH  AT  THE  TAIL,  €t,  DEG. 


FIGURE  15.  EFFECT  OF  FUSELAGE  ANGLE  OF  ATTACK 
ON  DOWNWASH  ANGLE  AT  THE  TAIL. 


DOWN  WASH  AT  THE  TAIL, «t, DEG 


FIGURE  16.  EFFECT  OF  ROTOR  LIFT  ON  DOWNWASH 


ANGLE  AT  THE  TAIL. 


MAXIMUM  LIFT  COEFFICIENT  FOR  STOL  AIRCRAFT: 

A  CRITICAL  REVIEW 
by 

D.  C.  Whittley 

Chief  Research  Engineer  (Aerodynamics) 

The  De  Havilland  Aircraft  of  Canada,  Limited, 

Downsview,  Ontario,  Canada. 

INTRODUCTION 

Various  investigators  have  studied  wing  theory  in  an  attempt  to 
establish  the  absolute  maximum  lift  coefficient  which  can  be  generated 
from  a  given  wing  of  finite  aspect  ratio  This  problem  is  an  important 
one  to  those  engaged  in  STOL  technology  because  take-off  and  landing  per¬ 
formance  depends  vitally  upon  it.  Also,  it  is  important  to  establish  the 
"absolute"  performance  potential,  so  that  reasonable  objectives  can  be  set 
for  research  programs  and  proper  conclusions  drawn  from  design  studies. 

The  present  paper  reviews  the  previous  work  on  the  subject  and 
draws  different  conclusions  in  some  respects.  Many  of  the  ideas  expressed 
are  not  regarded  as  absolute  but  are  presented  to  promote  discussion. 
STATEMENT  OF  THE  PROBLEM 

It  is  well  known  that  CLmax  is  a  parameter  of  prime  importance  in 
the  calculation  of  take-off  performance.  It  could  be  said  that  wing  stall 
usually  occurs  prematurely  because  of  boundary  layer  and  viscous  effects 


so  that  the  full  potential  of  the  wing  is  seldom  attained  in  practice.  How¬ 
ever,  various  boundary  layer  control  systems  (suction  and  blowing)  provide 
a  means  of  largely  overcoming  the  adverse  viscous  effects  so  that  the  full 
potential  of  the  wing  can  be  more  nearly  attained.  The  questions  which 
arise  are  these:  what  is  the  full  potential  lift  of  a  wing  and  how  nearly  can 
it  be  attained  (even  with  the  help  of  boundary  layer  control  devices)  ?  And, 
again,  what  is  the  drag  (particularly  induced  drag)  associated  with  this  lift? 
Further  to  this,  the  question  of  super -circulation  is  considered  briefly. 

Various  investigators  have  attempted  to  predict  the  limiting  value  of 

CLmax  (*  "  ®)  while  others  have  used  these  predictions  to  obtain  absolute 

(6) 

potential  performance  of  STOL  aircraft  (for  example,  Sutcliffe,  and 
Johnston  &  Saint,  All  the  investigators  ^  ”  ^who  have  developed  "large 
angle"  theories  suggest  that  the  absolute  value  of  CLmax  on  a  wing  depends 
linearly  on  aspect  ratio.  That  is,  CLmax  =  (constant)  A.  The  value  of  the 
constant  arrived  at  varies  from  0.  855  to  1.  94  and  so  gives  very  wide  limits 
for  the  possible  absolute  value  of  CLmax* 

Also,  it  is  quite  apparent  that  formulae  of  this  kind  tend  to  quite 
astounding  limits,  as  aspect  ratio  becomes  large.  For  instance,  the  form 
used  by  Sutcliffe  states  that  for  an  aspect  ratio  of  20,  the  maximum  poten¬ 
tial  aerodynamic  lift  coefficient  is  31.  4! 

The  present  note  suggests  that  expressions  of  this  kind  do  not  give  a 
sensible  guide  to  the  potential  value  of  CL.max  (except  at  quite  low  values  of 
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aspect  ratio  outside  the  range  of  our  particular  interest)  and  suggests 
another  limit  which  may  be  more  realistic. 

BRIEF  REVIEW  OF  LARGE  ANGLE  WING  THEORY 


Various  investigators  have  made  theoretical  predictions  of  the  lift 
generated  by  finite  wings  at  high  incidence  ^  “  ®).  These  theories  are 
unanimous  in  that  they  predict  a  limit  to  the  lift  which  can  be  achieved. 

This  limit  arises  simply  because,  at  high  a  ,  the  local  force  vector  is 
inclined  so  far  backward  by  the  local  downwash  that  the  ability  of  the  wing 
to  produce  lift  becomes  seriously  impaired,  (the  direction  of  the  lift  force 
being,  by  definition,  at  right  angles  to  the  distance  up-stream  relative  wind). 

For  wings  with  elliptic  spanwise  loading  some  limiting  values  of  lift 
have  been  given  as  follows: 


CL  =  .  855  A 
CL  =  1.  21  A 
Cl  =  1.  90  A 


(2)  (3) 
(4) 

(1)(5) 


Analysis  shows  that  the  differences  between  these  various  forms 
depend  upon  the  assumptions  which  are  made  concerning  the  final  inclination 
of  the  vortex  wake  at  infinity  down  stream. 

The  first  expression  is  based  on  the  assumption  that  the  vortex  sheet 
remains  flat  (i.  e.  it  does  not  roll  up)  and  has  a  transversal  velocity  equal  to 
ro/b  ;  in  this  case  the  vortex  sheet  follows  a  convex  path,  the  downwash  at 
infinity  down  stream  being  greater  than  that  at  the  wing. 
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In  the  second  expression,  it  is  assumed  that  the  vortex  sheet 
remains  flat  and  follows  a  straight  path  such  that  the  downwash  at  infinity 
down  stream  is  the  same  as  that  at  the  wing.  This  case  arose  because  it 
happened  to  represent  a  convenient  mathematical  approximation  but  it  has 
only  weak  physical  justification. 

The  third  equation  results  from  the  assumption  that  the  vortex 
sheet  rolls  up  and  finally  attains  a  stable  configuration  in  the  form  of  two 
symmetrical  vortex  cores  of  opposite  and  equal  circulation  r0  ;  the  dis¬ 
tance  between  their  centers  being  given  by  *-  b/  4  where  ’b’  is  the  wing 
span.  The  transversal  velocity  at  infinity  down  stream  is  less  than  that  at 
the  wing  and  the  vortex  sheet,  therefore,  follows  a  concave  path  as  it  moves 

down  stream.  The  transversal  velocity  in  this  case  is  equal  to  \  —  - 

L  b 

This  latter  case  has  ample  physical  justification  and  is  equally  good 
in  terms  of  mathematical  accuracy.  It  would,  therefore,  appear  to  give 
the  preferred  equation,  namely,  limiting  Cl  4  1.  90  A. 

(5) 

The  derivation  of  this  equation  was  first  given  by  Helmbold  in  1956  . 
In  retrospect,  it  would  appear  to  be  a  brilliant  piece  of  manipulation.  The 
data  of  (3) (4)  whiCh  came  four  or  five  years  later  appear  very  much  as  an 
anti-climax,  particularly  since  the  first  and  second  equations  can  also  be 
derived  directly  from  Helmbold' s  analysis  (See  Appendix  "A"). 

The  work  of  Cone  (1961)  elaborates  on  the  "rolling  up"  process  and 
describes  it  in  much  more  mathematical  detail  so  as  to  provide  general 
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information  about  the  downwash  field  in  the  vicinity  of  the  empennage.  He 
arrives  at  substantially  the  same  result  as  Helmbold  with  regard  to  the 
limitation  of  maximum  lift  coefficients;  namely,  Cl  =  1.  94  A. 

Appendices  "A”  and  "B”  deal  with  the  large  angle  wing  theory  in 
more  detail.  Appendix  "C"  investigates  the  implications  of  large  angle 
theory  on  the  interpretation  of  wind  tunnel  test  data,  while  Appendix  ’D” 
speculates  about  the  possible  ratio  between  translational  and  rotational 
energy  in  the  wake. 

THE  WORK  OF  NICHOLS  AND  SUTCLIFFE 

/g\ 

Sutcliffe  '  '(page  191)  reviews  the  simple  theory  put  forward  by 
Nichols  It  is  a  momentum  concept  in  which  the  lift  and  drag  of  a  wing 
of  finite  aspect  ratio  is  related  to  the  downwash,  £w  ,  at  the  wing.  This 
would  appear  to  be  incorrect  because  lift  and  drag  depend  upon  the  final 
changes  experienced  by  the  associated  mass  and  not  on  the  local  changes 
which  occur  at  the  wing.  Consideration  should  be  given  to  the  final  con¬ 
ditions  of  the  flow  at  infinity  down  stream  to  obtain  a  correct  result  of 
simple  momentum  theory,  (See  for  example,  Nichols’  theory,  there¬ 

fore,  cannot  be  considered  very  seriously.  Comments  to  this  effect  were 
made  by  the  writer  on  page  231  of  ^  as  part  of  the  discussion  to  Sutcliffe's 
paper.  However,  errors  in  drag  are  not  very  great,  provided  Cl/a  is  kept 
small;  for  example,  Sutcliffe  limited  Cl/a  ^  0-  6. 

Nichol’s  equations  are  given  for  reference,  as  follows: 
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jr  A  „  •  a 
=—  sin 

(1) 

cDi 

C0S  fw] 

(2) 

whereas  more  correct  results  have  been  derived  from  simple  momentum 
theory  as  follows:  (See 

CL  =-y  sin  £„  (3) 

CDi  (4) 

One  obvious  objection  to  these  simplified  forms  of  momentum  wing 
theory  is  that  the  energy  dissipated  in  the  wake  is  assumed  to  be  purely 
translational.  In  fact,  we  know  that  the  wake  takes  the  form  of  a  strong 
contra-rotating  vortex  pair  having  some  translational  motion,  so  that  a 
large  part  of  the  energy  is  bound  up  in  this  rotational  motion. 

A  PROPOSED  LIMIT  TO  MAXIMUM  LIFT 

It  has  been  pointed  out  that  Helmbold's  limit  to  maximum  lift 
(Cl  =  1.  90  A)  yields  exceedingly  high  .alues  of  lift  coefficient  for  moderate 
values  of  aspect  ratio  (say  A  =  10). 

A  more  practical  limit  to  the  lift  of  a  wing  might  well  be  obtained  by 
assuming  that  the  local  flow  at  the  wing  is  deflected  to  90°;  values  greater 
than  this  being  arbitrarily  considered  extreme.  In  effect,  this  corresponds 
to  an  unflapped  wing  at  90°  angle  of  attack  in  potential  flow.  In  practice, 
similar  flow  conditions  would  more  probably  be  achieved  with  a  flap  angle 
of  75°  with  the  aerofoil  at  30°  incidence. 

This  suggested  limit  also  corresponds  to  the  condition  when  the  two 
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stagnation  points  approach  coincidence  (if  A  =  oo  ),  and  any  further  increase 
in  circulation  would  cause  the  stagnation  points  to  leave  the  surface.  Under 
these  circumstances,  fluid  must  circulate  right  around  the  aerofoil,  which 
it  is  not  likely  to  do  in  practice  if  the  trailing  edge  is  sharp.  Thus,  a  second 
limitation  to  lift  could  be  considered;  namely,  the  value  of  lift  achieved  when 
the  stagnation  points  coincide.  Such  a  limit  would  permit  local  flow  angles 
in  excess  of  90°  (for  A<oo).  On  the  basis  of  these  suggestions,  it  is  possible 
to  obtain  expressions  for  limiting  Cl-  The  ’lift’  of  an  elliptic  wing  of  finite 
aspect  ratio  may  be  written: 


Cl  =  2t  (1  +t/c)  sin  (  a  -  £w  )  (5) 

where  (1  +t/c)  is  an  approximate  correction  factor  to  the  sectional  lift -curve 
slope  to  account  for  wing  thickness  chord  ratio  (t/c).  The  ’lift’,  (to  which 
corresponds)  acts  at  right  angles  to  the  local  free  stream  velocity).  Let 
a  =  90°  represent  a  practical  limit  to  lift,  then 

limiting  Cl  =  2»  (1  +  t/c)  cos  (6) 

Now  following  an  analysis  similar  to  that  of  Appendix  "B": 

Local  velocity  =  V  cos 

Also  lift  L  -  L*  cos  £*w  (8) 


Therefore,  Cl 
and  equation  (6)  becomes: 


h' 

-^pV2cos2£wS 


cl 

cos3£w 


giving: 


limiting  Cl  =  2*  (1  +  t/c)  cos 

limiting  (^L/A)  _  2  n  (1  +  t/c) 

cos*£w 


(9) 

(10) 
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But  <?w  is  a  known  function  of  ^L/A  (See  Figure  A  2,  Appendix  ’’A”),  and 
therefore,  t-JZ?)  plots  uniquely  against  ^/A  according  to  equation  (U) 

(See  Figure  1).  The  variation  of  limiting  Cl  with  aspect  ratio  and  thickness 
chord  ratio  is  shown  in  Figure  2,  this  being  based  on  a  lift-curve  slope  of 
2  jt  (1  +  t/c)  according  to  equation  (5). 

H  coincidence  of  the  stagnation  points  is  used  as  a  basis  for  limita¬ 
tion  of  lift,  then  the  average  local  flow  direction  close  to  the  wing  exceeds 
90°  and  an  alternative  expression  for  limiting  Cl  is  given  as  follows: 

limiting  Cl  =  2 r  (1  +  t/c)  cos3£w  (11) 

The  implications  of  such  an  assumption  are  also  shown  in  Figure  1. 

It  should  be  noted  that  the  previous  derivations  have  dealt  with  the 
limiting  lift  without  physical  chord  extension  devices. 

Leading  edge  devices  are  an  essential  part  of  high  lift  technology 
(nose  blowing,  slats,  etc. ).  Leading  edge  slats  may  appear  to  give  a  chord 
extension  but  seldom  give  an  effective  increase  in  wing  area.  They  do, 
however,  represent  a  very  powerful  means  of  extending  the  lift  curve  to 
higher  angles  of  attack. 

Flap  extension  schemes,  such  as  the  Fowler  flap,  have  been  used  at 
moderate  values  of  lift  coefficient  in  spite  of  the  large  nose  down  pitching 
moment  which  they  generate.  However,  since  tail  lift  coefficient  required 
to  trim  increases  linearly  with  wing  lift  coefficient,  it  seems  unlikely  that 
high  lift  flaps  using  targe  flap  extensions  will  find  a  place  in  STOL  technology. 
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THE  WORK  OF  CLARENCE  D.  CONE 

Cone  M  makes  a  lengthy  analysis  of  the  vorticity  in  the  wake  and 
obtains  a  limit  to  lift  coefficient  given  by: 

CLmax  =1*94  A  (12) 

Once  again,  we  note  that  other  investigators,  using  similar  concepts,  have 
arrived  at  formulae  of  similar  form  but  with  a  different  constant.  For 
example,  both  Ribner  and  Hancock  give  Cl  =  0.  855  A  The  relative  merits 
of  these  various  expressions  are  discussed  in  the  section  on  Large  Angle 
Wing  Theory,  and  in  Appendix  "A". 

However,  Cone  does  recognize  that,  in  the  limit,  expressions  such 
as  equation  (12)  yield  unrealistic  answers  in  attempts  to  provide  a  more 
practical  limit  for  wings  of  aspect  ratio  greater  than  6.  5.  This  he  gives  in 
Page  31,  as  CLmax  =  4*  .  Cone  chose,  as  a  limiting  factor,  the 
conditions  which  cause  front  and  rear  stagnation  points  to  coincide,  and 
considers  a  Joukowski  profile  as  obtained  by  conformal  transformation  of  a 
circle  of  radius  'a'  with  the  chord  of  the  profile  being  given  as  *c\  He  then 
gives  an  expression  for  lift  as: 

Cl  =  8?r  a/c  sin  a  (13) 

This  equation  happens  to  be  true  for  many  (if  not  all)  profiles  which  have 
been  transformed  from  a  circle  of  radius  'a'  to  a  profile  of  chord  *c\ 

Cone  then  goes  on  to  state  that  the  maximum  value  of  a/c  is  1/2,  and 
so  for  this  case:  Cl  =  4*  sin  a  (I4i 


9 


which  in  the  limit  yields  CLmax  =  4*  =  12.  6. 

The  first  thing  we  note  here  concerning  equation  (14)  is  that  for 
small  angles  the  lift-curve  slope  tends  to  4*  (rather  than  the  more  familiar 
value  of  2*  ).  On  further  investigation  it  is  found  that  the  symmetrical 
Joukowski  aerofoil  can  be  obtained  in  a  full  range  of  thickness/chord  ratios. 

At  one  extreme  it  becomes  a  flat  plate  for  which  a/c  =  1/4  (a  well  known 
result);  at  the  other  end  of  the  range  the  transformation  degenerates  to  a 
one  to  one  correspondence,  and  the  flow  round  the  circle  of  radius  ’a’  remains 
unchanged.  It  is  this  case  for  which  a/c  =  1/2. 

It  would  appear  that  Cone  unwittingly  chose  this  rather  uninteresting 
case  to  obtain  equation  (14)  -  it  simply  expresses  the  lift  on  a  circular 
cylinder  as  the  stagnation  points  move  toward  one  another.  The  lift 
coefficient  corresponding  to  coincidence  of  the  stagnation  points  on  a 
circular  cylinder  is  well  known  to  be  equal  to  4  *  . 

Cone  limits  the  aerofoils  lift  coefficient  of  aspect  ratios  greater  than 
6.  5  to  CLmax  =  4*  .  It  would  seem  that  a  value  of  2*  would  be  closer  to 
reality  for  those  engaged  in  the  design  of  conventional  types  of  aircraft.  If 
limiting  Cl  =  2‘n'  is  accepted  as  a  more  realistic  limit,  then  vortex 
considerations  would  limit  the  lift  attainable  at  values  of  aspect  ratio  less 
than  3.  3.  (In  other  words,  the  formula  C^max  =  1.  94  A  would  not  apply 
above  an  aspect  ratio  of  3.  3. )  The  use  of  such  relationships  should,  there¬ 
fore,  be  avoided  in  take-off  calculations  for  conventional  STOL  transport 
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aircraft.  Similarly,  the  ’’complete"  induced  drag  polar  shown  by  Cone 
Figure  23,  and  discussed  by  Schairer,  has  little  meaning  in  this  context, 
although  it  might  have  a  bearing  on  the  Concorde,  for  example. 

Figure  25  of  Cone's  report  further  substantiates  the  verity  of  the 

o 

standard  form  of  the  induced  drag  equation  Cdi  =  Cl  /it A,  up  to  the 
maximum  values  of  lift  that  are  likely  to  be  achieved  in  practice,  provided, 
once  again,  that  we  confine  attention  to  a  moderate  or  high  aspect  ratio  wing. 
OTHER  LIMITING  FACTORS 

It  has  been  suggested,  that  a  practical  limitation  to  lift  on  a  wing  of 
moderate  aspect  ratio  would  occur  when  the  local  flow  deflection  reachs  a 
value  of  90°  or  that,  in  any  case,  it  is  not  likely  to  exceed  a  value  corres¬ 
ponding  to  coincidence  of  the  stagnation  points.  Either  of  these  conditions 
would  appear  to  be  quite  extreme  when  one  tries  to  visualize  the  corres¬ 
ponding  flow  patterns  round  the  wing,  and  so,  the  useable  value  of  lift  is 
likely  to  be  somewhat  less.  Also,  lift  limitations  would  be  imposed  for 
other  quite  different  reasons;  these  are  briefly: 

Effect  of  Fuselage 

Lift  coefficients  are  normally  based  on  gross  wing  area  and  there 
will  be  a  'lift  loss"  due  to  the  blanking  effect  of  the  fuselage.  Lift  will  not 
be  completely  lost  in  this  region  but  will  fall  far  below  the  peak  potential  as 
compared  to  other  sections  of  the  wing. 
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Effect  of  Trim 


High  lift  coefficients  usually  create  large  nose-down  moments,  and  a 
significant  down  load  on  the  tail  is  then  required  to  trim  the  aircraft.  This 
down  load  detracts  from  the  overall  lift  coefficient  which  would  otherwise  be 
achieved. 

Part-span  Flaps 

It  is  not  always  convenient  to  employ  full-span  flaps.  Although 
ailerons  may  be  drooped,  they  seldom  are  as  effective  as  the  flap  for  two 
reasons:  first,  in  order  to  retain  roll  effectiveness,  droop  angles  are  often 
limited  to  approximatley  40°;  secondly,  it  is  not  always  possible  to  incor¬ 
porate  quite  such  exotic  high  lift  devices  in  the  aileron  (as  compared  to  the 
flap)  because  the  aileron  is  required  to  serve  a  double  function.  Once  again, 
the  lift  tends  to  fall  below  the  maximum  potential  value. 

Effect  of  Flap  Drag 

A  very  large  flap  angle  often  produces  a  large  drag:  this  hinders 
take-off  and  can  also  cause  embarrassment  in  the  event  of  a  single  engine 
failure.  Consequently,  it  often  turns  out  that  the  'optimum'  flap  angle  for 
take-off  is  considerably  less  than  maximum. 

Use  of  B.  L.  C. 

It  becomes  increasingly  costly  in  terms  of  air  supply  to  maintain 
attached  flow  at  large  flap  angles  when  using  conventional  blowing  for 
boundary  layer  control.  Once  again,  it  may  be  desirable  to  restrict  the 
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flap  angle  or  settle  for  conditions  in  which  the  flow  over  the  flap  is  fully 
attached. 

No  attempt  has  been  made  to  estimate  the  loss  of  lift  due  to  such 
effects  since  each  design  will  differ  and  must  be  considered  on  its  own 
merits.  It  is,  nevertheless,  interesting  to  note  that  large  scale  experiments 
have  been  carried  out  on  a  transport  type  aircraft  which  employed  quite 
powerful  B.  L.  C.  flaps.  The  maximum  lift  coefficient  reported  with  a  flap 
angle  of  70°  at  a  =  13°  was  approximately  3.0 
SUPERCIRCULATION  AND  THE  JET  FLAP 

In  the  previous  sections,  consideration  has  been  given  to  the 
maximum  lift  coefficient  which  can  be  achieved  on  a  wing  with  ’potential" 
flow.  It  was  assumed  that  potential  flow  conditions  could  be  closely  attained 
by  means  of  conventional  flap  blowing  schemes,  together  with  some  form  of 
leading  edge  device.  The  amount  of  flap  blowing  was  assumed  to  be  just 
sufficient  to  meet  the  requirement  (CM  critical^  so  *^at  super  circulation 
effects  were  not  significant. 

Wing  Lift  with  Jet  Flap 

If  a  jet  flap  (or  augmentor  wing)  scheme  is  used,  then  the  flap  is 
said  to  be  ’’over  blown".  This  is  the  condition  known  as  supercirculation 
and  large  values  of  lift  coefficient  can  now  be  achieved  that  are  in  excess  of 
the  limit  given  by  equation  (9).  The  reasons  for  this  can  be  described  in 
various  ways;  one  is  to  state  that  the  lift  of  a  two-dimensional  wing  depends 
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upon  the  circulation,  which,  from  a  mathematical  point  of  view,  is 
distributed  along  the  wing  chord  in  such  a  manner  as  to  satisfy  the  boundary 
condition  (zero  flow  normal  to  the  wing  surface).  In  the  case  of  the  jet  flap, 
a  jet  sheet  extends  behind  the  wing  and  sustains  a  pressure  change  across 
it.  Circulation  in  the  field  now  extends  along  the  wing  chord  and  along  the 
jet  sheet  so  that. the  total  integrated  lift  on  the  wing  (Spur)  can  now  be 
considerably  greater  than  for  the  aerofoil  without  jet  sheet. 

The  variation  of  lift  coefficient  with  jet  strength  and  jet  angle  has 
been  predicted  by  Spence  for  a  two-dimensional  jet  augmented  flap.  This 
theory  has  been  reasonably  well  substantiated  by  experiment  and  if 
anything,  the  lift  tends  to  exceed  the  prediction.  The  jet  flap  in  three- 
dimensions  has  been  considered  theoretically  by  various  investigators, 
including  Maskell  &  Spence  The  limitations  to  lift  which  might  apply 

in  this  case  are  not  easy  to  predict  because  the  problem  is  now  further 
complicated  by  the  presence  of  the  jet  sheet.  The  extension  of  the  work  by 
Cone  to  include  the  effect  of  a  jet  sheet  would  be  of  great  interest  in  this 
regard. 

In  the  absence  of  an  adequate  theory,  one  might  be  permitted  to 
speculate  as  follows: 

One  of  the  previous  limits  for  the  case  without  jet  sheet  was  set  by 
coincidence  of  the  two  stagnation  points.  This  gave  a  limit  to  lift  coefficient 
of  2ir  (approximately)  for  a  thin  aerofoil  of  infinite  aspect  ratio,  but  as 
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already  indicated.  thi»  will  no  longer  apply. 

The  equation  6-33  (c)^^  quotes  Strand  to  give  the  following 
expression  for  lift  of  a  jet  flapped  wing: 

Clt  =  2*  sin  «<  +(3.  9  Cj*1  +  Cj)  sin  (  0  +  a  ) 
where  e  is  the  jet  angle  relative  to  the  chord  line  and  the  theory  has  been 
extended  to  apply  for  large  angles  of  deflection. 

Lift  coefficient  would  reach  a  maximum  for  a  =  90°,  0  =  0°. 
Physically  this  is  an  unrealistic  condition  and  a  more  realistic  case  is 
arbitrarily  suggested  as  a  =  30°,  0  =  60°,  giving  (  0  +  a  )  =  90°. 

Thus,  iCLT)max  =  «  +(3.9  Cj*  +  Cj)  having  an  aerodynamic  lift 
component  of  (  ir  +3.9  Cj* )  as  compared  to  2  K  ,  approximately,  without  a 
jet  sheet.  This  suggests  an  increase  of  about  40°o  in  the  aerodynamic  lift  of 
a  two-dimensional  jet  flap  with  Cj  =  2  as  compared  to  the  potential  lift  of  a 
wing  without  jet  sheet. 

Wing  Drag  With  Jet  Flap 

Consideration  of  the  downwash  at  infinity  down  stream  and  the 
corresponding  magnitude  of  the  associated  mass  (according  to  momentum 
wing  theory)  leads  to  a  reconsideration  of  the  induced  drag  of  a  wing  with  jet 
flap. 

(13) 

The  well  known  result  given  by  Maskell  7  can  be  derived  as  follows: 

Assuming  that  the  jet  itself  is  subject  to  the  influence  of  the  vortex 
system,  it  will  acquire  a  downward  velocity  at  infinity  down  stream  equal  to 
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that  of  the  associated  mass. 


Then  lift  =  V0*  sin  £«,  +  M j  Vj  sin  £a 

giving  CLt  =(|irA  +  2Cj^-^^D 


sin  £r 


Also,  equating  work  done  to  energy  in  the  wake, 

°i  vo  :  )  vo(vo  sin  4>)  +  aT  vo 

where  AT  is  the  thrust  loss  due  to  deflection  of  the  jet  sheet  far  below  the 

wing. 

ACt  =  (1  -  cos  )  Cj  =  Cj  for  Sa  small 

Thus,  CDi  =.§10^0  ^A  +  2Cj) 
giving  CDi  = 

VirA  +  2  Cj/ 

Alternatively,  and  in  a  similar  manner,  the  results  of  Appendix  "D" 

can  be  used  to  obtain  a  different  relationship  between  lift  and  drag  of  a  jet 

flapped  wing  with  elliptic  spanwise  loading.  The  associated  mass  is  now 

o 

larger  by  a  factor  _j(_  ,  the  downwash  at  infinity  down  stream  less  by  a 

2 

factor  of  2  ,  and  the  ratio  of  total  energy  to  translational  energy  carries 


a  factor  (7f  ). 

Thus,  in  a  similar  manner  as  before: 


Lift  =  p 


(J2")(JT )  vo sin  4 +  M,  vj sin  4 


giving  CLt  =  „A  +  2 

Aiso,  D,  V0  (£)’($)  V0(V0  sin  <?„)“  ♦  .T  V0 
TTius,  CDi  -  2  Cj] 
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The  modifying  factor  in  the  square  brackets  is  very  close  to  unitv 


for  combinations  of  Cj  and  A  of  interest. 
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APPENDIX  "A" 

(51 

The  work  of  Helmbold  v  '  is  repeated  here  for  ease  of  reference  and 
because  of  its  fundamental  value. 

"With  a  spanwise  elliptic  distribution  of  circulation,  T =  rJiHiY/b)1 
the  trailing- vortex  system  far  downstream  of  the  wing  will  move  like  a  rigid 
flat  strip  with  the  proper  velocity  normal  to  itself,  provided  that 

its  cross  section  is  hypothetically  endowed  with  the  ability  to  resist  the 
deforming  stress  of  the  pressure  fieid.  The  transversal  momentum  carried 
with  the  vortex  sheet  increases  per  unit  length  by  , 

and  the  kinetic  energy  of  the  absolute  motion  increases  per  unit  length  by 

a/2). 

If  the  cross-sectional  rigidity  is  instantaneously  destroyed  and  the 


sheet  is  left  free  to  yield  to  the  deforming  stresses  of  the  pressure  field, 
the  sheet  will  start  rolling  itself  in  from  its  side  edges  and  finally  attain  a 
stable  configuration  in  the  form  of  two  symmetrical  vortex  cores  of  opposite 
and  equal  circulation  H,  .  The  distance  between  their  centers  of  gravity  is 
“Is  ftb/A  ,  and  their  proper  velocity  is 

w  =  Ti/2*  t*2  n/rr’b  *  (2/r*)wo  (i) 

at  least  approximately.  Hence  W0  =  ( 7(  /2)  W  .  Since  no  outer  forces 
are  acting  on  the  vortex  system  during  the  ’’deformation,  the  transversal 
momentum  and  the  total  energy  per  unit  length  remain  unchanged  by  the 
deformation; 

d  Jx/6  ac  *  6  Jzm/6x.  and  d E/dac  *  dE0/dx0 
(With  the  deformed  vortex  system,  the  total  energy  consists  of  kinetic 
energy  and  work  done  by  pressure  forces;  with  the  undeformed  vortex  sheet, 
the  kinetic  energy  is  the  total  energy  because  then  there  are  no  perturbation 
velocities  along  the  vortex  sheet. ) 

’’The  completely  rolled- in  vortex  system  is  inclined  by  the  angle  6 
toward  the  undisturbed  relative  velocity  V.  This  angle  is  determined  by 
sin  6  =  w'/V.  The  length  of  the  final  vortex  system  increases  per  unit  time 
by  dx/dt  =  V  cos  6  .  The  components  of  the  aerodynamic  force  on  the  wing 
are  the  lift. 

L  =  V  •  (dJz/dx)  -  Djtg  4  (2) 

and  the  induced  drag  -  Dj  =  (dE/dx)  cos  &  (3) 
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Since  DjV  =  dE/dt  =  (dE/dx)  V  cos  6 


‘According  to  preceding  statements, 

L  =  /o(*bj4)Vw0  yO  (Tb^Xw.’/aJsiN  8 

=  /3  0r/fc)bVl$jN  6  [l-Or/'-O  sim*^) 

(4) 

D[  -  yO  COS  S 

=  /3(x9/3z)  bV’sii/s  cos$ 

(5) 

a  a 

or,  after  division  by  b/2 

Cl/A  =  {**/*)  SIN  s] 

(6) 

CDi/A  =  (n/'l6)  SIN *8  COS  & 

(7) 

A 

Here  A  =  tr/S  denotes  the  aspect  ratio.  The  reduced  lift  coefficient  Cl/A 
attains  its  maximum  value,  Cl  lim/^  =  2/(3  J~Z)  =  1.  90, 

for  sin  8  =  2/(iT  J~3)  =  0.  368  with  a  corresponding  value  of  the  reduced 
drag  coefficient 

CDl/A  =  (n2/l 2)  J^r^A/3  =  2.  40 

It  will  now  be  shown  how  the  various  assumptions  concerning  the 
transversal  velocity  of  the  jet  sheet  at  infinity  down  stream  can  be  used  to 
obtain  the  three  equations  described  in  the  section  on  Large  Angle  Wing 
Theory.  Helmbold's  notation  is  retained  so  that  the  downwash  angle  at 
infinity  downstream  is  written  as  8  ,  rather  than  as  in  other  parts  of 
the  text. 

Helmbold  gives: 

L  =  /ofab/^Vwo  -/o(7rb/4Kw«J/2)siN  S  (A- 1) 
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(A-2) 


Di  -yO('irb/<)(w'0/i)cos  S 


Three  cases  are  now  considered: 
(i)  SIN  8  *  Wo/ V 
UO  $iN&*W0/2V 
(iii)  SIN  8  *  (Wo/vKZ/X*) 


corresponding  to  Cl=  0.  855A 

(A-3) 

corresponding  to  Cl  =  1.  21 A 

(A- 4) 

corresponding  to  Helmbold 
and  to  equation  Cl  =  1.  90  A 

(A- 5) 


Substituting  these  values  for  W0  in  equation  (A  -  1),  we  obtain: 


(i)  CJ->  =  SIN  *•  °'815  (A- 6) 

yielding  a  limiting  Cl  =  0.  855  A 

=  7t[sIN  8  -  SIN3S]  SIN  6*  0.576  (A- 7) 

yielding  a  limiting  Cl  =  1.  205  A 
(iii)i^L=  ^-3£siN  8  -(f  )siN36j  Sin  8  =  0.368  (A-8) 

yielding  a  limiting  Cl  =  1.  90  A 

Corresponding  equations  for  induced  drag  may  also  be  obtained  thus: 
Helmbold  gives  Dj  =  /°M(  ¥)  COS  8  (A- 9) 

Using  equations  (A-3),  (A-4),  and  (A- 5),  we  have: 

=(jt)  sin's  COS  S 

<“>  £5L  =  ^)siN'Scos6 

A 

(fjsiN’ScoSS 

Each  equation  for  drag,  with  its  corresponding  equation  for  lift  gives 
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a  lift  to  drag  relationship  shown  in  Figure  A-l. 

The  variation  of  £w  with  Cl can  be  found  using  tan  Sw  =  Di/L 

(from  equations  A-6  to  A-12),  and  this  is  shown  as  a  plot  of  sin  £wvs  Cl/A 

in  Figure  A-2  for  each  of  the  three  cases.  Prandtl’s  linear  theory  gives 
CL 

sin<?w  =  -  which  is  also  shown  in  Figure  A-2  for  comparison. 

'n'A 

APPENDIX  "B" 

Appendix  "B"  demonstrates  how  equation  A-7,  Appendix  "A"  can  be 
derived  using  well  known  results  from  Prandtl’s  lifting  theory.  This  serves 
to  give  a  more  physical  insight  into  the  mechanism  of  lift  vector  rotation 
which  causes  the  limitation  of  lift.  It  also  provides  a  basis  for  equations  (6) 
to  (11)  of  the  main  text. 


It  is  assumed  here  that  the  wing  has  elliptic  spanwise  loading  and 
that  the  trailing  vorticity  lies  in  a  plane  inclined  to  the  free  stream  and  does 
not  roll  up  -  the  inclination  being  equal  to  the  induced  downwash  angle 
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calculated  at  the  wing  lifting  line.  In  practice,  the  vortex  sheet  leaves  the 
wing  at  this  angle,  and  the  basis  of  the  present  assumption  is  that  the 
deviation  from  it  is  insignificantly  close  to  the  wing  where  the  effect  of  the 
sheet  on  the  wing  is  predominant.  This  is  simply  a  convenience  for  the  sake 
of  mathematical  simplicity  but  nevertheless  it  sounds  reasonably  plausible. 
(This  assumption  leads  to  limiting  Cl  =  1.  21  A,  as  previously  shown). 

Let  the  local  'lift'  L  be  inclined  backward  through  the  (downwash) 
angle  Sw 


Then  lift  L  =  h'  •  cos  £w 

cl  - 

(B-l) 

and  drag  D  =  if  •  sin  <?w 

(B-2) 

and,  therefore  Cd  =  Cl  tan 

(B-3) 

cLv 


The  transversal  velocity  of  the  vortex  sheet  is  usually  given  as 
.  Now,  in  large  angle  theory,  these  values  of  Cl  and  V  must  assume 


local  values  which  account  for  inclination  and  for  local  changes  in  the 
velocity. 

Thus,  (CLV)  LOCAL  =  s  •  ^ 

therefore,  sin  £w  =  %  *  (E 


•  WJ  V.  y  _  v-i.  v0 

V.*  cos  *<fw  S  J  cos*£w 

COS 

(B-4) 

cos  ^ 

(B-5) 

-  sin^  £w  ) 

(B-6) 

This  being  identical  to  equation  (A-7)  which  gives  rise  to  a  limiting 
value  of  Cl  =  1.  20  5 A. 
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APPENDIX  "C" 


This  appendix  investigates  the  implications  of  large  angle  theory 
with  regard  to  experimental  measurements  of  lift  made  on  wings  of  finite 
aspect  ratio  when  comparisons  are  required  with  two-dimensional  theory. 
In  such  experiments,  the  measured  lift  for  a  given  vorticity  is  less  than  it 
would  be  using  linear  theory  because  of  the  inclination  of  the  force  vector 
(the  effect  of  which  is  usually  small,  and  therefore  ignored). 

Helmbold  gives: 

l  =  /°(&)  sm  S 

where  W0  can  be  considered  as  a  measure  of  vorticity 

since  T.  »  w0  b 

Let  L  =  /a  (ZJL)  Vw„ 

where  h'  is  the  lift  which  would  result  using  small  angle  theory. 


Now  W°  i  Z  <V 

V  T* 


•  _  cL'r 

Evaluating: 

CL/A 

.5 

1.0 

1.  5 

(Cl  /a)2 

.  25 

1.0 

2.  25 

.01 

.02 

.09 

Cl"  /cL 

1.01 

1.04 

1.  10 

24 

.495 


.96 


1.36 


Cl/A 

-  thus,  giving  an  approximate  correction  factor  for  the  purpose 


described. 


APPENDIX  ’D" 

Ratio  of  Translation  to  Rotational  Energy  in  the  Wake 

It  has  been  noted  that  one  obvious  objection  to  simple  momentum 

wing  theory  is  that  all  of  the  residual  energy  in  the  wake  is  assumed  to  be 

translational.  Making  this  assumption,  then  the  associated  mass  is  given  by: 

,  this  value  is  consistent  with  equation  (4)  of  the  main  text. 

For  &  small,  Helmboid  gives  Cl  =  sin  S  (equation  A- 8)  as 

compared  to  Cl  =  ^  sin  6  for  Prandtl  (equation  A-6),  where  6  is 

the  downwash  angle  at  infinity  down  stream.  This  implies  that,  for  a  given 

wing  lift,  the  associated  mass  given  by  Halmbold  is  greater  by  a  factor  of 

_.a 

-j—  ,  and,  that  the  downwash  is  correspondingly  smaller  by  a  factor  of 

approximately. 

The  total  residual  energy  in  the  wake  gives  rise  to  induced  drag  and 
it  is  known  that  the  Prandtl  result  is  substantially  correct  in  this  regard. 

It  is,  therefore,  possible  to  deduce  the  ratio  of  translational  to  rotational 
enerty  in  the  wake  according  to  Halmbold’ s  model. 

E  (translational)  +  E  (rotational)  =  E  (total) 

Let  -m  represent  the  associated  mass  per  unit  time  which  is 
deflected  according  to  Prandtl,  and  let  W0  represent  its  corresponding 
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translational  velocity: 

Then  -J-  s-yvncj/ 

giving  ♦  k„J 

when  Hr  =  k,  m  U0* 

Thus,  k„  the  fraction  of  rotational  energy  is  given  by  (1  -  ), 

that  is  nearly  80%  of  the  wake  energy  is  rotational. 

*  The  total  energy  in  the  wake,  shed  from  an  elliptically  loaded  wing,  is 
known  from  wing  theory,  (for  example,  Sears,  1955,  Page  227).  This  is 
derived  from  consideration  of  an  undeformed  vortex  sheet  at  infinity  down 
stream.  The  fact  that  the  vortex  sheet  might  roll  up  into  a  closed  pair  does 
not  alter  the  wake  energy.  (This  was  pointed  out  by  Halmbold).  Thus,  when 
considering  momentum  wing  theory,  the  wake  energy  is  known  and  does  not 
depend  upon  assumed  changes  in  the  value  of  downwash  (  )  or  correspond¬ 

ing  changes  in  the  associated  mass. 
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